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I. Introduction 



The stopping power of a material is defined as the space-rate of 
loss of energy suffered by a particle of nuclear or atomic dimensions 
in its passage through the material. Since the time of Thomson and 
Rutherford the subject has been investigated both theoretically and 
experimentally in its many complexities, and it early provided a 
source of tests of atomic structure and mechanics. In recent years 
the interest of theoretical nuclear physicists in the various facets 
of the subject has waned, since it is largely a property of the extra- 
nuclear part of the atom. Nevertheless the subject is still a vital 
one in the experimental and applied side of nuclear physics. The 
need for good stopping power formulas and/or data for all types of 
nuclear particles, over the whole range of energies, and throxigh all 
types of matter is almost too obvioiis to mention - e.g, , in cosmic 
ray studies, range determinations, shielding requirements, and so on. 
Unfortunately for those who need this information, stopping 
poir/er formulas provided by theoreticians are in most cases unable to 
give accurate predictions without the use of semi-empirical constants 
adjusted to conform to previous experiment. This is due to the 
physical and mathematical complexities of the subject. By this date 
theory and experiment are sufficiently advanced that for heavy par- 
ticles, such as protons or alpha particles, the stopping power can 
be predicted with accxiracy for most types of matter usually encountered, 
provided the energy of the particles exceeds a few million electron- 
volts. However if one passes to the energy region lower than this, the 
predictions becone less accurate and even empirical adjustments are 



difficult. In such regions experimental data must largely be relied 
on at present. Such data are still far from complete, even for 
stopping by basic elements. 

For direct measurement of stopping power within this energy 
region the Van de Graaff generator is well suited, because of the 
range of accelerating potentials available and because of its ability 
to produce accelerated particles of precisely defined and accurately 
measured energies. For this dissertation the ?n:iter has elected to 
report on the experimental results of stopping power of various ele- 
ments for protons from the Van de Graaff generator at the Ohio State 
University. Specifically, measurements are herein reported of the 
stopping pov/er for protons of the metals, copper and nickel, and of 
the gases, nitrogen, neon, argon, krypton, and xenon. The energy 
range over which these meastirements are taken is from 400 to 1200 kev. 

Brief mention is made of the associated problem of "straggling" 
in stopping power, or the mean square deviation in energy loss per 
unit thickness of matter for Individual protons. This phenomenon is 
a result of the fact that slowing down of particles is accomplished 
throxigh random collisions vjith atcmis of the stopping medium and is 
therefore subject to statistical variations. This factor is less 
important than that of average stopping power, and both theoretically 
and experimentally it has been studied less fully. Accurate analysis 
of straggling is just as difficult in theory as analysis of stopping 
power itself, and experimentally it is even more difficult, requiring 
generally a precision of measurement better than stopping power 
measurement by an order of magnitude. The technioue employed by the 



writer for determination of stopping does make it possible to deter- 
mine variations in sane, under favorable experimental conditions. It 
has been found possible to report stopping power straggling for the 
metallic elements, copper and nickel, although precise accuracy is 
lacking. Such results are considered as a "by-product" of the major 
investigation into stopping power itself. 

A, Definitions and Units, 

Stopping pov;er has been defined above; and unless otherwise 
specified or implied it refers to the average energy loss of a large 
number of similar particles of the same energy, per unit length of 
path through the material traversed. It is dependent upon the type 
of particle, the energy (or velocity) of the particle, and the stop- 
ping medixm. The concept of stopping pov/er involves a continuous 
change in particle energy, and yet a specified value of stopping 
poiver is valid only for a fixed value of particle energy. This in- 
dicates that stopping poxver must be recognized as a derivative, 

- dE/dx, itself a function of energy. 

If the stopping medium is composed of a single element, one can 
use the concept of atomic stopping pov/er, defined as stopping power 
divided by the number of atoms per unit cube of the stopping material, 
(This is sonetimes called the" "stopping cross-section" of the atom.) 

It may be considered as the average contribution of each atom of the 
elemental material toward slov/ing the incident particles. In symbols 



Theoretical formulas have been derived in the Gaussian (c.g.s.) 
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systCTi of units, so that stopping power is given in ergs /cm. For 
experimental work - both in determination and use of stopping power - 
energy changes are more usually expressed in electron-volts or 
similar units, while thicknesses are conveniently expressed in mg/cm , 
Thus stopping power data from experiment is usually reported in units 
of kev-cn^/mg. The conversion factor between the two system of units 
is easily seen to be: 



Stopping power (ergs/cm) 



1.602 X 10“ X Stopping power 

(kev-cm^/mg) x^o , (I-A-1) 



where 



Z’ = density of medium (g/cra^) . 

For atomic stopping power the proper c.g.s. \mits are erg-cm^ • 

p 

Experimental conventions make preferable the units ev-cm , The 
conversion factor betv/een the ltwo systems is obviously: 

p Ip 

At* stopping poiver (erg-cn^) ® 1.602 x 10***"^ x At. stopping 

po\7er ( ev-c m^) . (I-A-2) 

The conversion factor from stopping power to atomic stopping 

power in the c.g.s. system is: 

(At. stopping power) s 

cgs 

A 



(Stopping power) ggg / N 



.23 



X (Stopping 



6.023 X 10 x^ 

power) , (I-A-3) 

cgs 

where N is the number of atoms per cubic centimeter of stopping 
material, and A is the atomic weight of the stopping element. 



From the above equations it is easily deduced that in experimen- 
tal units the conversion factor from stooping power to atonic stopping 



pov/er is: 

(At, stopning pov/er) 



X (stopping 



exp. 



6,023 X 

power ) 



exp. 



(I-A-4) 

The concept of stopping pov;er "straggling" is understood most 
easily by visualizing a statistical distribution of relative fre- 
quencies of stopping T)ov;er losses for individual particles, under 
identical conditions. If the standard deviation about the mean of 
this distribution is designated as - 0 ., then beccnes the mean 

square fluctuation , or the variance; and one best describes the 

straggling by ~ ( JX^) . (In finite increments we may consider 

dx 

straggling as , v/hich is used rather than since, 

as we shall see later, the former is Independent of thickness of 

medium v<hile the latter is not.) This quantity may also be vTitten, 

as sha-m by statistical theory (Hoel, 1947), as ^ (p;^) - ('E)^J , 

where E in this expression indicates the energy of an Individual 

particle among a large population of similar narticles incident on 

the stopping medium with the sane energy. 

The units for 4(il" ) are evidently ergs^/cm in the d.g.s. 
dx 

systeri. In experimental work the units employed ere kev^-cm^/mg . 

The conversion relation betv;een the values of straggling as exnressed 
in these two systems is: ^ 






. 1.602 



cga 



10 



15 



r 



A iR”) 

dx exp. 



(l-A-5) 



B, Revier/ of Stopping povrer Theory. 



A complete revie;v of stopping pov/er theory would be almost im- 



possible and certainly impractical for a report of this nature. It 
is necessary therefore to present only a brief outline of the theory. 
Much discussion can be eliminated by concentration only on those 
aspects bearing on our particular problem. For example, relativity 
considerations may be Ignored, since v;e are interested only in a 
range of velocities well below that of light (v/c 1/30). 
Furthermore, we need consider only the stopping of heavy particles 
such as iirotons, and primarily by the nediuia or heavy elements. 

In a broad sense there are two approaches to the problem, and 

1 

each will be summarized in turn. The first, exemplified by the 
v;orks of N. Bohr and his follov;era, is classical or seml-clessical 
in spirit. This is advantageous for ease of understanding, mathe- 
matical simplicity, and flexibility In application to a wide variety 
of conditions. The most rigorous formulation of the subject has been 
achieved, hw/ever, by modern v?ave-mechanical methods, with their 
mathematical complexities. Such methods may be difficult to ap- 
preciate intuitively, but they yield more acctirate results in those 
realms v;here the beslc assumptions are valid. 

In presenting the simpler theory we shall follow primarily the 
excellent treatment given by H, Bohr (1948), 

1, Classical Fundamentals. 

Let us consider the elastic collision between an incident 
particle (number one) and a second particle which is at rest in the 

For the sake of clarity and brevity we are not following strictly a 
historical, or chronological, scheme of presentation. All major 
contributors to the theory however will be duly acknowledged where 
their ideas best fit into this summary. 



laboratory frame of reference. Fig^r^e 1 shows the geometric rela- 
tionships existing between the velocities in the center-of-mass 
system. 




In this figure, 

"r - velocity of particle 1, laboratory system, before collision; 

V = vel, of center-of-raass, laboratory S 3 ^tem, before collision; 

V 

-V ■ vel, of particle 2, center-of-nass system, before collision; 
c 

v-v s vel, of particle 1, center-of-mass system, before collision; 
c 

9 m angle of deviation of each particle in the c,-m. system. 

(For direct collision, the particles ore always oppositely directed in 
the center-of-mass system, so that each particle undergoes the sane 
angular deviation in this system.) 

Each particle behaves in the center-of-mass system as if it 
rebounds elastically from the immovable center of mass without loss of 
energy. Then, 



V ’ ■ vel, of particle 2, center-of-mass system, after collision 



and 



hsi = Kl = 



Also, 



^ n _ 



(I-B-l) 



Tg s Tel, of particle 2, labOMtory system, after collision, 

s T2*+ ▼c • (I-B-2) 

We see that Tg foms the base of an isosceles triangle with 
rertei angle 0 . Therefore 

Tg = 2 Tg sin (6/2) . (I-B-3) 

From the definition of center-of-mass it is easily shown that 

Tg « ^ y (I-B-4) 

mi-MDg 



where and nig are masses of the tv/o particles. Thus 

Tp = 2 t sin (0/2) . (I-B-5) 

Hovr in the collision particle 2 is given energy at the expense 
of particle 1 (laboratory coordinates) , Let T be this energy trans- 
ferred. 



T - ^ 1H2 



(I-B-6) 



and from equation (I-B-5), 

T . ^ . T^ sin^(e/2) 



s T^ 8in2(6/2) 

tq2 










(I-B-7) 



where 1s the "reduced mass" and equals m].ng/(mj^ -h mg) • 

The maximum transfer of energy occurs on direct collision, when 
0 equals 180® • Let T^^^ be this amount. Then 



m 



and 

T 



^ y2 . 

®2 

Tj^ sli^(0/2) 



(I-B-8) 

(I-B-9) 



If the law of attraction or repulsion between the two particles 



is CoulOTibian, one nay derive the well laioi’m Rutherford scattering 
forraula: 

e 

CBC^ (e/2) • d<o , (I-B-10) 

\2 vy 

v/here 

do- = differential probability cross-section that the particle 
is scattered into solid angle, do> , measured in the 
center-of-raass coordinate system, 
dco a total solid angle corresponding to scattering directions 
between 6 and e + d9 , 

> 2-rsine de ■ 4'Tr sin(e/2) • cos(e/2) - de . 

(I-B-11) 

■ charge of the first particle, 

6g = charge of the second particle. 

Substitution of the expression (I-B-11) into equation (I-B-10) gives: 



d o- 



4TT 




2 

- csc®(d/2) • cos(G/ 2) • dO , 



(I-B-12) 



From equation (I-B-7) we obtain by differentiation: 
dT * T cos(e/2)' csc(e/2) • de . 

Then from the last two equations we obtain: 




(I-B-13) 



(I-B-14) 



vrtiere D is the quantity witliln the parentheses. 

We may use this result to obtain stopping power relations in a 

highly idealized case. Let us suppose that an incident ^rticle 

undergoes a succession of energy-diminishing collisions in traversing 

a thin layer of matter in a direction normal to the layer, 

_ a 










■K-' - I iiM 

piJillMMiHM. tiB iBt ii 1 1 i'i'i'ilwta' • 
Ej*!lPwiWlj>»ff |l:||^.'i 1 1 1 1 



[v", inJ: 






Let 

/HiX s thickness of layer of matter; 

N* s number of jwrticles per unit cube of the stopping 
medium that absorbs the energy. 

It is necessary to assime here that: 

(a) The particles of the stopping medium are free particles at rest, 
interacting with the incident part, through the coulomb law of force 
between charged particles. 

(b) llie energy transfer in an individual collision Is small compared 

with the total energy of the incident particle. This is true provided 
Tjjj « ^ , v/hlch according to equation (I-B-8) is true when 

(c) The medium layer is very thin so that v is essentially constant 
throughout the progress of the incident particle v/lthln the layer. 

Assumption (b) indicates that manentum change of the incident 
particle is small in a single collision, vftiich permits one to state 
that the i»rtlcle velocity is almost constant in direction as v/ell as 
in magnitude. 

Now consider that we have made a series of hypothetical experi- 
ments in which we were able to folio?/ each particle in an Incident 
beam with common initial characteristics, and further that we could 
obtain and tabulate the data for each collision each particle tmder- 
goes with particles of the stopping material. Divide the possible 
values of energy loss in a single collision into equal intervals, 
with Tj^ corresponding to the average energy for the 1*^ interval 
and being the probability that each collision will provide an 



Since the 



energy loss within the interval characterized by 
stopping process is a statistical process because of the randomness 
of location of particles of the medlxim in relation to the path of 
some incident particle, v/e cannot always get in different trials the 
same nimber of collisions characterized by . Define 

n^^ * nvimber of collisions characterized by energy loss Tj , 

for the case of any one incident particle; 

■ average number of collisions characterized by Tj^ , for 
the total population of Incident particles in the beam. 
From the given definitions it is evident that 

= N»AX0i . (I-B-15) 

In a practical case the number of collisions is large and 
Individual "single shot probability" is small. In such cases the 
probability distribution function for n^^ can be closely apjaroximated 
by a Poisson distribution fonaula (Hoel, 1947) : 

P(nj^) 5 e“^i / n^! , (I-B-16) 

having mean value as required and having standard deviation equal 
to , 

Now if AE is the change in energy of a particular particle 
when passing through the given layer of the medium. 



-AE z 



T. n, 



1 1 

The average value of energy loss is 

-AE = ^ T, w. ; 
i 

and the mean square fluctuation is by definition 



(I-B-17) 




(AE - *3e 



* / 



( I-B-IB) 




I 



a . (I-B-19) 

i ^ ^ 

eqiiation (I-B-15) we obtain: 

, (I-B-20) 

. (l-B-21) 

Under the rather idealized conditions assimed we can obtain a 
formula for stopping pcv/er and the statistical variation thereof by 
substituting equation (I-B-14) in equations (I-B-20) and (l-B-21), 

In doing so we may make sane assumptions based on the particular 
physical sitiiation we are studying: 

(a) The primary type of particle in the medium responsible for 
stopping effects is the atomic electron. Nuclear collisions do 
contribute slightly to stopping ixDwer but such contribution is 
usually small (as v/ill be shown later for our particular case) , 

(b) The type of incident particle considered will be heavy, ns our 
v/ork is concerned with protons. 

As a result of these assumptions we may make certain modifica- 
tions in our formulas. If we designate the mass of the electron by 
m , then we may replace lUg by m , Tlie reduced mass m^ is also 
approximately equal to m . If vw designate the charge of the electron 
by e , then bee ones ze and Og is simply e . 







¥/e see thet the criterion in sub-paragraph (b) on page 10 is 
quite valid. 

If N is defined as the number of atons per unit volume of the 
stopping medium, vie may use it Instead of N* , provided we also sum 
over all electrons in the atcm. Thus N* ia replaced by N and a 
summation sign over s , v;here a is the index for each electron 
of the atom effective in the stopping phenomenon. If there is no 
reason to separate the contributions of various atomic electrons, or 
if their properties can be characterized by overall average values, 
then K* can be replaced by N Z , where Z Is the. atomic number of 
the element constituting the stopping medium. 

We obtain then for stopping power: 

In 



- ® = D • N 



dx 



^max 



^nin 



2 7T *2 ^ 



(I-B-22) 



m v^ s 

For stopping power straggling we obtain: 



In 3222 
\in 



^ ( XT) s D N Z (T - T ) 
ux max min 

s' D N Z T , 

max ’ 



(I-B-23) 



provided (as is visual) T » T . , and T is the same for all 

max min max 

atomic electrons. 

2. Stopping Power - Classical Formulation. 

Inherent in the above analysis is the assumption that the electron 
behaves as a free particle, v/hereas it is actually bound to the atom. 
Nevertheless under many circvcnstances this assumption may be nuite 



reasonable. For the present discussion let us assume that the particle 



velocity, V , Is greater than the "orbital velocity" of the electrons, 
(The classical approach is essentially the same whether the electron 
is considered as bound by elastic or Coulonb forces to the atom.) If 

u is the orbital velocity of the s^^^ electron, and I_ its lonl- 

s s 

zation energy, then 

I s' m u ^ , (I-B-24) 

s s 

It Is possible to divide the stopping power into tV'fo parts - 
(a) that portion due to close collisions of the particle with atoms, 
and (b) that portion due to more distant collisions. These must be 
considered separately, 

(a) Close collisions are those for vrtilch T>I . If T I , 

then for most close collisions, the ionization energy is so much 
smaller than T that it may be disregarded. This makes reasonable the 
assumption that for close collisions the electron may be considered 

as free. Under such clrcvmistances it is clear that in a direct 

collision, giving T , the electron is ejected v;lth a velocity 

max 

practically eqxial to 2v , so that 

T = 2 m v^ (l-B-25) 

%x 

Vie may note that oiir original assumption that v > u leads directly 

s 

to the conclusion that 2mv^ » feu_^ , or T I - thus 

confirming this assumption. 

T ^ for these collisions is of course I , Then equation 
min s 

(I-B-22) leads to 

. (I-B-26) 

I, 

(b) For more distant collisions the energy transferred by classical 




collisions v;ould be smaller than the ionization energy. Classical 
models consider the electron as an oscillator (real or virtual) 
vrtiich can absorb any amoxint of energy, hor/ever small, when set into 
vibration by a sudden impulse from a cha'nging force field, such as 
would be provided by a passing charged particle. The fact that v 
is greater than indicates that the oscillator has moved from 

its equilibrium position only slightly by the time the particle has 
possed. Thus the "restoring force" has not been mobilized to any 
extent dvirlng the collision time. Once more therefore the electron 
may be considered as undergoing a free collision, imder a situation 
however which peraits energy absorption v;itho\it ejection from the 
at<xn. 

The vnlue of T used in this region of stonning is obviously 
max 

the value of used in the previous case; to wit, 1^, 

The value of T . cannot be zero, since this would give an in- 
nin * 

finite answer for stopnlng pjower. It has been deemed necessary to 

presume the existence of an "adiabatic limit" to the distance from 

the atom to the prarticle path. It Is pxresumed that if a particle 

passes at a distance greater than this limit, the electric field is 

applied and released so slov/ly that the electron is eased out of and 

back into position v/ithout being left in an oscillatory (excited) state. 

The duration of the impulse from the particle field is on the 

order of d^/v , v/here d is the adiabatic limit corresponding to 

the s'^^ electron. If this duration is less than l/o_ , where oj 

° & 

is the oscillatory frequency, then the impulse is faster than the 
ability of the electron to absorb and return it in adiabatic fashion. 



-'3 - 



Thus we make 



(I-B-27) 



dg = t/c,3 , 

Then will be the value of energy transfer for collisions at 

this distance. It can be shown that for such distant collisions be- 
tween free particles (H, Bohr, 1948, Eq. 1,1,10): 



T 2 2^ , 1 

m 






(l-B-28) 



We see then with the use of equation (I-B-27) that 



min 



2 z® 



m V* 

Equation (I-B-22) leads to: 
F- is] = D N 

L ^Jb » 



( I-P-29) 



In 



Ig_m 



2 z2 g4 <0^2 



( I-B-30) 



The total stopping povrer Is given by adding the results of 
equations (l-B-26) and (l-B-30): 



dE 

dx 



2 D N 



In 



k m 

z e^t*>» 



(I-B-31) 



This is the classical formula proposed by N, Bohr (1913) . It 
should be noted that the constant k has been introduced arbitrarily 
into the formula. Insofar as the above considerations apply the value 
of k is unity, Bohr’s more rigorous approach gave to k a value of 



1.123 . 

3, Elementary Ouantxmi Considerations in Stopping Power . 

Unfortunately Bohr's formula gave values vdiich v;ere not in very 
good agreement with the experimental data on alpha particle sto-ining 
which v/ere available at that time. Early attempts to improve the 



formula were unsuccessful. 



The older quantum theory of atomic structure vdiich Bohr developed 



at about this same time established that the boxmd energy states of 
the atonic electrons are discrete. This threw doubt upon the classical 
oscillator concept vxhich permits the electron to absorb any emoimt of 
energy (within the restrictions of the adiabatic limit) and leads to 
formula (I-B-30), Henderson (1922a) suggested that only close , ion- 
izing collisions should be considered, which led to a formula eaulva- 
lent to equation (I-B-26) for the complete stopping power. This scheme 
liketvise failed for alpha port ides, but it is interesting to note that 

it gives anspjers equal to just about one-half of the experimental 

* 

values (Fowler, 1923). 

V;ith our present comprehension of the limitations of classical 

mechanics at atomic distances, it is not difficult to deduce wherein 

formula (I-B-31) fails. Vie must accept the quantum requirement that 

energy transfer in a collision raist be at least as large as I for 

0 

the s*^ electron, (For simplicity, we disregard excitation possi- 

bilitiesj Thus classical theory can be expected to hold only if 

T . , as given by equation (I-B-29), is greater than I . This 

min o 

criterion can be transformed by use of the relation 



Now if we define v as the velocity of the electron In the ground 

o 

state of the hydrogen atom, we know that 




(l-B-32) 



and the equation (I-B-24)^into the form 




1 



(I-B-33) 



V 



o 



e^/ Jf = 0,219 X 10® cm/sec , 



/7 - 



(I-B-34) 



By substitution in equation (I-B-33) and taking the square root, the 
criterion becomes: 



z V / V 
o 



V / u 



(I-B-35) 



Our original assumption v/as that V P>-u_ , so that all criteria for 

s 

applicability of the classical formula can be expressed as: 



z 



> ~ ^ 1 

s 



y ^ ^ - . (I-B-36) 

Bohr (1948) showed that this criterion is a little stricter than 

necessary, and that correct results can be expected fron the classical 

formula for values of z v^/v down to about unity . 

Unfortunately, with respect to protons or alpha ^articles, the 

whole stopping pov/er theory herein derived fails v/hen z v^/v is 

appreciably greater than unity, because the particles tend to lose 

their charge by electron capture under such circumstances (Bohr, 1948; 

also see Section 7, belov/). Thus the failure of the classical 

equation for such particles is obvious. 

For analysis of stopping power vSien v >z v (or for protons 

o 

V >Vq), we must once again divide the effect into Uio parts. 

(a) Consider the s^^ electron, for which equations (I-B-24) and 

(I-B-32) are still valid. This electron carin-^t be nrecisely located, 

*• 

but has a position v/hlch is smeared out over a distance having 



• order of magnitude a , where 

8 

a = X / nrug . 



(l-B-37) 



a is fairly conr)arable to atonic dimensions. P/e first consider the 
s ^ 

contribution to stopoing pov;er of those collisions v.'here the minimum 
approach distance of the Incident particle is less than a_, which 



means that this portion is comparable to part (a) of the classical 
approach. 

For .those fast electrons of the stopping material, in case of 
which u ^ T , v;e note that 

-L ^ JL ^ JL 

mu m T m - Y * 

s o 

or a < , where A is the de Broglie wave-length of the reduced 

8 

particle. ( ^ eqmls P\. /ZTT , just as equals h/2 ‘77' ,) 

Since in this discussion the approach distance is less than a , it 

s ' 

is always less than /A ; and we no longer have a well defined collision 
in the classical sense. Such electrons therefore have little stopping 
effect. 

For electrons in the case of which v^u„ (giving ®-) 
a situation which still permits a classical approach within certain 
limits. 

V/hen the incident particle penetrates the s*^ electronic cloud 
of the atom (that is, roughly, within the distance a ) the collision 

o 

may be a direct one, giving * 2m v^. However, many collisions 
will not be quite direct. The least amount of energy able to be 
transferred within this region is indicated by the uncertainty prin- 
ciple. Since we are considering only those collisions where approach 
is closer than , our uncertainty in electron position cannot 
exceed this amount. Then equation (I-B-37) Indicates that uncertainty 
in momentum cannot be less than m u^ . Tlxia in turn means that 
uncertainty in energy is not less than ^ (=Ig)» 'Ihus we can use 

the classical formulas v;lth T^.,^ = that 3^ 



f 



Just as In part (a) of the purely classical approach.) 

V/e thus obtain for stopping power in this region a formula which 
is the same as (I-B-26): 



[■«].■ f 



In 



2 m 



(I-B-38) 



Collisions in this category are called "free collisions" since 
most of them provide energy transfers greatly exceeding the binding 
energy of the electron, permitting it to be considered as free. 

(b) For collision approach distancds equal to , one can find the 
energy transferred by a formula similar to equation (I-B-28): 

T : 



*s 



2 2 ^ 1 

— — ? 
m v2 



■M 



(I-B-39) 



obtained with the assistance of equations (I-B-24) and (I-B-37). 

This indicates" that, for heavy particles with velocities well above v^ , 

T < I_ . Thtis, if we take a classical point of view, few, if any, 

®s ° 

collisions where the particle passes outside the atran can cause 
ionization. It is this fact which led to the difficulties of 
Henderson (1922) previously mentioned, and such difficulties were 
unresolved until modern v/ave-mechanlcs established the possibility 
of energy transfer through the "resonance" process. 

Previous to this, however, Fermi (1924) had espoused an idea which 
was to be fruitful. He proposed that the perturbing force exerted 
by the particle on the atom could be analyzed as a function of time 
into Fourier components, and the effect of each component on the 
atoms could be compared to the effect of electromagnetic radiation 



of the same wave-length. This analog 7 between particle collision 
theory and eleetronagnetic absorption theory nay not be valid for 
close collisions, but for distant collisions where the particle 
field is almost constant over the whole atom the analogy seems 
promising. Even the idea of the "adiabatic limit" can be carried 
over, since at that remote distance the perturbing field no longer 
contains oscillatory components of anpreciable magnitude which can 
effect any sort of resonance interaction with the atomic electrons. 

It is well known today that the classical theory of electro- 
magnetic absorption gives results almost the same as those obtained 
by a strict qtiantum approach and confirmed by experiment (Compton, 
1919; for general discussion see Compton and Allison, 1935), One 
may presume then that for distant collisions we may obtain a correct 



result by continuing the classical approach. Uiis gives us then, 
from equation (I-B-22): 




obtained with the use of equations (I-B-24) , (I-B-27), (I-B-32), and 
(l-B-37). 

Collisions leading to this component of stopnlng povrer are 

Called "resonance collisions". Just as in the case for free collisions, 

those electrons for which u >v are to be ignored, as their con- 

s 



i 

1 






tributlon aecordlng to the formula is quite small. In fact it can 



be seen that the formula even provides a negative (and ccinpletely 
unphysical) resvilt when 1^ > 2 mv^ , 

In general, then, the formulation for total stopping power is 
simply: 

-•g = ED-N^ In-L®-!^ , (I-B-41) 

s -^s 

being the sura of contributions given in equations (I-B-38) and 
(I-B-40) . 



One refinement may be inserted here, for which we shall give 
only a qualitative Justification. Hitherto the Interaction of the 
atomic electrons has been ignored, but because of the screening 
effects of the electrons on one another they cannot be considered as 
independent oscillators. It is proper to weight the contribution of 
each electron by a factor f , called the "oscillator strength" of 

o 

the s*^ electron. The value of each f_ is still on the order of 

S 

unity, and it has been shown (Kuhn, 1925; Thcsnas, 1925; and many 
others) that 



^ ^s = Z 

8 



(I-B-42) 



Thxis we see that equation (I-B-41) is still approximately correct, 
but more precise answers would be given by: 



~ M s 2 B N f In 

B 

. 4irz2 ^ ^ f in^nivS 



dx 



( I-B-43) 



m v2 8 -B 

Except for a slight modification of the term 1^, this formula checks 



with that found by Bethe (1930) , using a rather strict quantum- 

nechanlcal approach, similar to that given in the next section, 

- 



We should note here that KLoch (1933) t (7 another rigorous quantun- 
mechauical approach obtained a general formula which contained 

Bohr*s classical eouation (I-B-31) and Bethe*s equation (I-B-43) as 
limiting cases, according to whether 2 zt^/y is much greater or 
much less than unity. (He retained the qualification that u 

fl 

for all electrons.) As we have seen, the conditions we are Interested 
in are those which put us definitely in the quantvun region of appli- 
cability, so that the refinements of Bloch*s formula need not be 
considered. 

4. Strict Quantum-Mechanical Derivation of Stopping Power Formula. 

The derivation of stopping power formula vising quantvm 
mechanics nay vise any of several different approaches. The one 
outlined here is equivalent to that given by Bethe (1930, 1933). 

The system to be considered consists of an atom and a moving 
charged particle within a large cube of dimension L . The Interaction 
between the atcni and the moving particle is to be considered as a time- 
dependent perturbation. Before the interaction occurs, the particle 
and atom are completely Independent, v/lth the former having the 
wave function fca? a free perticle. The vreve-functlon for the initial 
state of the system can be written therefore as: 



where 




eag^ ( ^ Po*R ) 




Po * 



(I-B-44) 



M 



the initial momentum of the ''reduced" particle, with re- 
spect to the center-of-mass of the system, = U t ; 
the reduced mass of the particle-atom system 



R s the position vector of the particle, referred to the 
atonic nucleus; 

Tj s the position vector of the electron, referred to 

the atomic nucleus; 

(Tj) m the vrave function of the unperttirbed atom in its 
ground state. 

Uo will concern outselves with the probability of transition 
of this systen to a final state in vdiich the atom is left In a 
certain state (excited) characterized by the subscript n (this 
n*^ state of the atc^ nay be one in which the electrons are still 
bound or nay be in that continuum of states vdiere one of the elec- 
trons is emitted frcm the atom) , and the particle is moving in 
another direction v;ith momentum p . After the collision is 
essentially completed, the wave-function of the system is again 
that of atCTa and free particle without further interaction: 



U 



n 




ex£ ( ^ p.R ) • ^^(rj) 



(I-B-45) 



Let us use the v/ell known formula for first order transition pro- 
babilities (Schlff, 1949, p. 193; Land&, 1951, p. 172) : 



W 

where 

V/ 




(I-B-46) 



the number of particles per unit time scattered into a 
differential element of solid angle (to be precisely 
described belov;) with final momentum p , presuming that 
incident particles arrive at such a rate that the re- 
quirement of one incident particle per cube of 

- 




dimension L is nelntalned; 

the ddnsity of final energy states of the system (since 
the atran is left in a definite energy state, the final 
density of states of the system xaust be obtained by 
consideration of density of states available to the 
scattered particle hevlng momentum p ; the density is 
further cxirtailed by the fact that v;e must consider only 
possible directions of the final momentum vector con- 
tained v/lthln the element of solid angle specified in 
the definition of W) ; 



H* z |U_H’ U dx « v/here the integration is carried 
on J o n ' 

out over all space coordinates of all narticles in the 
system. 

To find , we note that the number of momentum states avail- 
able to a free particle definitely contained in a cube of dimension 

L, with momentum within the increments dp , dp , and dp is 

X y z 

(L^/h^) dp dp dp . This is equally applicable to any element in 

X y z 

momentum space, so that if vie ore interested In the number of states 
of the particle v/ith an absolute value of momentum betiveen p and 

3 2 

p ^-dp , the number of momentum states available is (L/h) p dp 
per unit solid angle. Our solid angle shall be restricted by con- 
sidering only those collisions wherein the scattering angle in the 
center-of-Huass system is betvjecn 0 and 6 +- d6 , so that the solid 
angle is 2 TT sin 0 • d0 . 

The number of energy states In the energy interval dE^ is the 
same es the nvunber of momentum states in the momentiun interval dp , 



\ 

\ 



provided the intervals correspond. Thus 

d£^ = (LA)® dp ( 2 'Tr sin « • d©) . (I-B-47) 
“Hie relation betv;een dE^ and dp can be obtained frcan the expression 
for total final energy of the system; 



and 

dp 

Therefore, 

/^E 



■— -f E (atom) ; 
811 n ' * 



P 

IT 



2'rrii p 



( I-B-48) 



(I-B-49) 



sin © d0 



( I-B-50) 



Substitution of equation (I-B-50) into equation (I-B-46) gives; 



Vi 



27r 



p sin 0 d9 



I"®-'!" 



(I-B-51) 



Let us now define d $(0) as the differential acattering cross- 

n 

section v/ith angle betv/een 0 and 0 + d0 , for one particle inci- 
dent per second upon one square centimeter, and leaving the otom in 
the n^^ excited state. This differs frcm W in that vre new have 
l/v particles per cubic centimeter whereas before we had 1/L^ 
particles per cubic centimeter. Then 

d^(0) . (L^/v) W ■ It-Ji v; 



. S. Bln e d9 
27rjf(* Po 



Kr 



(I-B-52) 



This result is the sarae as that obtained through use of the Born 



approxination, to v/hich this approach is eauivalent; and therefore for 
for its validity it depends upon those criteria v;hich make the Porn 

-^6 ' 



B 







approiiniation valid, (Sohlff, 1949) 

Let us reduce . The perturbing portion of the Hi^miltonian 

is the interaction potential energy betv/een the particle and the atom: 




(I-B-53) 



The first tern in the bracket can be neglected, as it gives a result 
for elastic scattering (corresponding to what Bohr colls nuclear 
scattering) wiiich may be ignored for protons or alpha particles 
insofar as it affects energy loss. 

It can be shovm (Bethe, 1930) that by integration over the 
space coordinates of the particle, can be reduced to the 

folloiving form: 



H* 



on 




4ir 




^n 



(l-B-54) 



where 



£ 



n 




> fo'in • 

( I-B-55) 



Substitution of the above expression into equation (I-B-52) gives; 

2 



8 TTI I^ 



= , ,, 
l2c -P| 

Define a flev; parameter: 



(z e^r ■ sin © de |£„|2 . 

Po ' 



(I-B-56) 
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Figure 2 
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Fran the figure vie see that 

q2 s - 2 p cos © . (I-B-57) 

For a given p^ and p , q is a function of 6 only. By dif- 
ferentiation we can derive that: 



sin © d© s SJlSL . 

P Po 

Then (I-B-56) can be ivritten: 



(I-B-58) 



- ~(-r - ^ l^«l" • 

q V 

This expression can be xvritten more simply and in a form per- 
mitting easy coninarison with the classical approach given above, if 
v;e let: 



Q 



. l£o-2l 



2 m 



2 n 



(I-B-60) 



Introduction of Q in equation (I-B-59) to replace q as parameter gives 



= 



2 TT z2 -4 



e" dQ 



K\ 



(I-B-61) 



i^'e recognize this as being in correspondence with enuation 
(I-B-14) obtained by strictly classical means, provided ve relate Q 
to T and provided is approximately unity for a one-electron 

atom. It is easily whoi'/n that Q , fran its definition, is the energy 
the secondary electron would receive from the collision provided it 
was not bound to the atom. 

To obtain the total probability of a transition which leaves th6 
atom in the n^ excited state, v/e must integrate over all possible 



values of Q : 



( I-B-6S) 



f 



n 




Let us define as the energy of the atom In its ground state, 

and remember that N equals the number of atoms per cubic centimeter 

of the stopping material. The total energy lows of the particle ner 

atom (atomic stopping pov/er) is the average energy gain per atom, 

obtained by multiplying i by (E_ - E ) and then sunning (in- 

^ o 

eluding integrating in the continuum of energy states) over all 
values of n available for excitation. If v/e multiply by H, we 
then get the total stopping power: 



dx 



O-N (En - V \ 

“ ■'r. 

%ita 




(I-B-63) 



For high values of energy transfer, the ionization energy can 
be ignored; and as vie have previously seen, Q T . Ihen for 

we may use s’ 2 m , by equation (I-B-25), (This is 

not accurate for Iw values of n since the binding energy of the 
electron has too strong an Influence to be Ignored in such cases. 
Under normal circumstances these low levels give only a very minor 
contribution to stopping power, so that the ap’^roxlmatlon for 0 

max 

gives a negligible error overall.) 

For Q , we can easily derive on axnression from (I-D-60) : 
min 



min 



2 m ^ 



(l-B-64) 



ra 4 E 

F\irther development wotild be facilitated by reversal of the 






-tS,: 











order of- sxiinmation and integration. As the eqiiation stands this is 
not possible since is denendent on the particular value of 

the quantum number n . However let us split the integration into 
two perts - vjitfi the middle limit, being independent of n . 

The only other restriction on 0^ is that it be less than the 
average Ionization potential of the atomic electrons concerned in the 
stopping. 

Let us determine first of all the stopping pov;er formula, 
considering hydrogen as the stopping medium. If Ry is defined as 
the ionization energy of the hydrogen atcm (13,6 ev) , we can use 



Q = (Ry)^/ 2rav^ 



(I-B-65) 



Then we obtain: 
dE 

"■S 



D-N 



j 

f ^ ("n-V K\ 

Qn 



+ 



n 







(I-B-66) 



In the evaluation of the first integral, we can use a theorem 

2 



proved by Bethe (1930), v/hich indicates that (E^^ - E^) j£ j - Q 

Evaluation of the second integral con be simplified since 0 and 

o 

are both much less than Ry , This indicates that Q is very 

small; condequently J^o “ £ small and approximately perpendicular 

to the path of the incident particle. The situation becomes closely 

equivalent to the case of energy absorption from electromagnetic v;aves. 

In such cases the sauare of the matrix element can be shopm to 

n 



equal the optical oscillator strength for the transition from ground 

-30 - 



to excited state of the at cm. Bethe (1930; 1933) uses this 

simplifying fact and carries out a conplete numerical solution for 
the stopping poiver of hydrogen. In doing so he uses the well known 
wave functions for hydrogen - expressed in polar spherical cocardinates 
for the discrete levels, and in parabolic coordinates for those levels 
in the continuum above. He obtains: 



n 




l^nl ” “ 0*098 



(I-B-66) 



Then, 




Hydrogen 




- 0.098 



] 



e 4 rr N ^ 2 . (I-B-67) 

m 1.103 Ry 

This formula is fairly acciurate and is confirmed approximately 
by experimental v/ork (Hano, 1934) . It should be noted that in its 
derivation it has been necessary to do the same thing as vms done in 
sub-sections 2 and 3 above - that is, to divide the stopr)lng effect 
into two parts, that involving close collisions vjhere the binding 
energy of the electron is negligible, and that involving distant 
collisions v/here analogy to electromagnetic dispersion and absorption 
theory is utilized. 

For heavier elements the sitviation is a little more complex. 

Bethe (1930) shoi/s that for each electron of the stopping material, 

the above formiila is valid, provided that the factor 1.103 Ry is 

replaced by a term I* , which is somewhat larger than the ionization 

8 

energy of the electron. Further, in summing the effect of each 



- 3 ' - 



electron, it must be weighted by Its oscillator strength 



The determination of I» is rather complex, as its precise value 
depends upon a knowledge of the wave-functions for the atom concerned. 
The formula for complex atoms nay be stated as follows: 

-i -- ^ f, lol^ . (i-B-ee, 

^ 4?: . ^ e« , N 2 , (i.b-69) 

n v^ ^av 

where I , called the "average excitation energy," is theoretically 

av 

defined by: 



Z In I ^ a In I* . (I-B-70) 

av — ' a 0 ' ' 

8 

Formula (I-B-68) is of course identical with (I-B-43) above, except 
that I* is probably a little larger than Ig , the true ionization 
energy. In practice it is necessary to select empirically the value 
of or the values of one or more to make the formula agree with 
experimental data. 

5» Modifications to Theory for Slo v ? Particles Through Pleavy ITedia. 

IJost of the above theory v;as developed on the basis that the 
charged particle_has a velocity faster than the "orbital velocity" 
of each of the atcmiic electrons in the stopping medium. That the 
criterion falls for protons having energies in the region we are 
studying is easy to show. If we consider the velocity of an atonic 
electron given by equation (I-B-24), then the energy of a proton 
having the same velocity as an atomic electron of ionization energy 

I is given by {ll/m)I . The criterion for validity of the vtnmodi- 

8 6 

fied equation for stopping power becomes: 






E >> (M/ni)I for evfery electron of the atom, where E Is the 
s 

actual proton energy, 

Hjom tables of ic^nization energies (Slegbahn, 1931), it is easily 
computed that the criterion fails: 

For yN , K' electrons, vAien E < 0.72 Her; 

For 0 , K electrons, vrhen E ^ 1.0 Mev; 

O 

For 2 o^®» ^ electrons, v/hen E <1,6 Mev; 

Fot A ,K electrons, when E < 5.8 Mev; 

18 

For , K and L electrons, when E < 0.35 Mev; 

For g 0 Ni, K and L electrons, when E <1,5 Mev; 

For „-Kr, K, L, and M electrons, v/hen E < 0,48 Mev; 

For K, L, and M electrons, when E Cc. 1,5 Mev, 

What does one do if the criterion fails? One approximation 
(Livingston and Bethe, 1937) involves the assumption that all elec- 
trons for which E >{M/m)I are fully effective in stopping, and 

0 

the other electrons do not contribute at all. This means that 
formula (I-B-69) is usable provided Z* is used in place of Z, 
where 



the summation being taken only over those electrons slower than the 
incident particle. is also different in such cases, since only 

those outer electrons conforming to the criterion are used in the 
averaging. 

This method gives a curve of proper shape but requires experinen- 




f, 



n 



(I-B-71) 



n 



tal values at various energies to "provide precise acciiracy. Further- 
more, unless I is readjusted accordingly, discontinuities are likely 
av 
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to appear In the curve at the points where E = (M/n)Ig , since 
Z* changes at these points. 

Bohr (1948) , who was less interested in precise results than 
in general trends, evaded this difficulty. He recognized in his 
analysis which led to the equation called herein equation (I>^-41) 
that those electrons for v/hlch u >v should be ignored. To 

o 

determine those electrons xvhlch are effective in the stopping process, 
he used a foraiula based on the Thomas-Perral statistical model of the 
atom { 'Thomas, 1927; Fermi, 1928); 

n(u^) s the mmiber of electrons in the atom with velocity 
smaller than ; 



S Ug/Vjj . (l-B-72) 

Equation (I-B-41) is modified by integrating over dn(u_) instead 

8 

of summing over s ; 



- ^ = 2 n N 



•2v 



In 



2 m v^ 



frm Ug^' (= Ig) 



<3n(u_) 



(I-B-73) 



The upper limit is that value for which the logarithm goes to zero. 
Simplification and integration then gives; 



8 D N (▼/▼«) 



dE 
dx 

_ 16 Tf z2 N 
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(I-B-74) 



Bohr states that the formula he xises for n(u_) , is not very accurate 

O 

if u_ represents the velocity of the very outer or very inner elec- 

O 

trons. Thus formvila (I-B-74) may be expected to be most valid for 
those particle velocities which are comnorable to electron velocities 
in the intermediate shells. 



Neufeld (1950) attempted to improve Bohr’s approach so as to 



provide ansivers more accurate than that given above. He went back 
to the basic formulation wherein stopping effects are separated into 
two parts - the contribution from "free collisions" and that from 
"resonance collisions". For the latter he used a more elaborate 
formula which had first been proposed by Fermi (1940) . Itis formula 
is noteivorthy since it provides positive contributions from eveiy 
electron, regardless of velocity. Furthermore it is convergent if 
all distant collisions are included, even out to a collision distance 
of infinity. This obviates the need for the "adiabatic limit" 
concept. 

Neufeld’s approach is tedious and requires a separate set of 
calculations for each type of stopping element. His calculated 
values for a feiv metals seem to be about lOf) less than the best 
experimental restilts available at the moment. The advantage of his 
approach is that it requires no previous experimental work to 
dettonlne the parameters. 

V/e have reserved till last the discussion of a more exact 

% 

treatment of this problem, which v/as initiated by Livingston and 

Bethe (1937), As a prelude to this contribution, it should be noted 

that shortly after Bethe (1930) gave his stopping potver formula, Kott 

(1931) and Henneberg (1933) shov/ed that the Born method of treating 

collision problems could be applied with validity to collisions of 

heavy particles with heavy atoms having electronic velocities faster 

than the particle velocity. This is true for atoms of high atonic 

nximber, since the perturbing potential of the incident particle is 

small compared to the binding potential for the fast, inner electrons. 

- 3 ^- 



(see also Mott and Massey, 1949.) This fact gives one faith in the 
use of previously derived formulas even when the original criterion 



that V > u„ is not completely adhered to. 

Suppose one v;iahee to obtain results which are valid for par- 
ticle velocities dovTn to the order of magnitude of the velocity of 
the K-shell electrons, or somewhat lower. An exact solution of the 
contribution of the K electrons is obviously much better than an 
assumption requiring either fvill contribution or zero contribution 
of the K electrons to the stopping effect. We may note that fairly 
simple wave-functions are available for use here, since even in a 
complex atom the innermost electrons obey wave-functions vrhlch are 
almost hydrogen-like. The formula used by Livingston and Bethe is 
thus simply a modification of the equation used by Bethe to derive 
the stopping power of hydrogen (equation I-B-63) : 



In this formula one integrates instead of summing over the 
various energy changes, since most of the discrete levels well below 
the continuum can be considered as filled, llie scale of the values 
of energy and Q are changed to put all calculations on a uniform 
basis for different values of Z : 




(I-B-75) 



where 







(I-B-76) 
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whera 
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effective nuclear charge as felt by the 


K electrons, 
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Z - (inner screening constant) 


(Slater, 1930) 
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which is actually the sane as eouation (I-B-64). 



is the Icfwer limit of E* , inasmuch as we consider transitions 
to the discrete levels as impossible by reason of their bein'? filled. 
Infinity is permitted as upper limit for Q, since the computations 
of the matrix element make the probabilities exceedingly small for 
any contributions by values of Q above 2mv^, 

Livingston and Bethe have evaluated this formula for by 
numerical Integration. For -B- equal to 0.7 , v/hich is ap’^roximately 
correct for values of Z from 6 to 13, they have plotted Bjr against 
^ as argument. can also be expressed as: 

«-r 0.7 , ^ ) r 1.81 In 3.63/^ - Cg-(^ ) . (I-B-82)^ 

where C-^ is a corrective tern to the basic formula (I-B-68 or 

JX 

I-B-69) which approaches zero as particle energy (or ^ ) increases. 
Ctirves for C„ are also given by Livingston and Bethe. 



It is easy to shay that the factor 1,81 is approximately the sum 
of the oscillator strengths of the tv/o K electrons, while 
2 m v^ 

1.103 ‘Ry 



3.63 
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They have also given a curve for atomic stopping ijower of air for 
protons, data from v/hich is tabulated in Table I, over the energy 
range in v/hich ive are interested. 



This formulation of K electron stopping povjer has been ex- 
tended by Brcivn {1950) and i'Jalske (1952) to values of •©• higher 
than 0.7 , v/hich In effect extends the theory to heavier atoms. 

They have made the computations and supplied graphical results for 
Bg( -e-,^) and ) . 

It is pertinent to ask about the corrections required for 
L, H, and higher electrons as the velocity of the incident particle 
approaches their orbital velocities. This is particularly Impor- 
tant, as vie have seen, for the heavier elements. The exact analj'-sls 
for stopping pov/er of electrons in these shells is v/ell nigh im- 
possible to carry out because of the ccmiplexity and lack of knowledge 
of the v/ave-fvmctions at these shells. 

Hirschfelder and Kagee (1948) generalized the foregoing method 
to apply to shells higher than the K shell. They utilized the 
original curves provided by Livingston and Eethe since those of 
Brown or Walske were not available at that tine. In effect they 
noted that the stopping number B , v/here 
2 m 



B = Z In 



- C - (higher shell corrections), 
av ^ (I-B-83) 



can be expressed as: 



B = , 



(I-B-84) 



where 



th 



X number of electrons in the 1 shell; 



~^S 



r effective oscillator strength of nn electron of the 



i^^ shell; 



th 



= the essential variable for the shell, depending 

upon ^ , 

For the K shell it is seen from eouation (I-B-82) that 



V = ^ 

K 1.81 



In (3.63 ~ 
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and, as before, 
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(l-B-87) 



• The generalization was made in the assumption that, in form , 

b is the same as b for ’all values of i , It is only necessary 
1 K 

to use a different value of I** (the ionization energy v;ithout 
external screening for the shell considered) in computing the 
argument ^ v/ith which to use the curves for . For all but the 



outer shell, the value of !*• was taken to be: 
Z •Ry 

. 



(I-B-88) 



with n^ being the effective principal quantum number of the 1 
shell. 



th 



It is hardly necessary to point out that such a procedure is 
approximate - primarily because of the marked deviation of the wave- 
functions of the middle shells from a hydrogen-like character. The 
originators of this approach hoped to predict the trends rather 
closely however; and the semi-eraplrical character of Bethe’s original 
formulation was retained by adjusting the value of I** for the outer 



Table I 



Semi-emplrical Calculations of Atomic Stopping Fewer 
Air - Livingston and Bathe (1937) 

Argon- Hirschf elder and I'agee (1948) 

Xenon- Hirschf elder and Hagee (1948) 



E 

(Mev) 




cr - ( io “^ 


er-cm^ ) 


Air 


Argon 


Xenon 


0,3 


12,93 


20,834 


35.5828 


0.4 


10.63 


17.601 


30.6465 


0,5 


8.86 


15.432 


27.2376 


0.6 


7.70 


13,833 


24.7589 


0.7 


6.90 


12.604 


22.7846 


0.8 


6.30 


11.610 


21.1880 


0.9 


5,81 


10.792 


19,8589 


1.0 


5.45 


10.108 


18.7458 


1.1 


5.11 


9.522 


17.7874 


1.2 


4.80 


9.011 


16.9314 


1.3 


4.53 


8.563 


16.1795 


1.4 


4,29 


8.160 


15.5181 


1.5 


4.07 


7.807 


14.9227 






shell of electrons as necessary to confom to experimental values of 
range in the media (gaseoiis) selected for study, ^ 

The stildies of Hirschfelder and Magee were applied to two gases 
of partictxlar interest here - argon and xenon. The computed results 
for atomic stopping po»/er of these tv/o gases for protons are tabu- 
lated in Table I , over the range of energies in which we are interested. 
6. Effect of the Physical and Chemical State of Media on S topninQ: '^ower. 

In all of the theory presented previously the atom which absorbs 
the energy has been considered to be in a free state, such as a 
monatomic gas. Any sort of chemical combination or condensation of 
the medivm into a liquid or solid state might be expected to shift the 
energy levels, at least of the electrons in the outer shells, and thus 
affect the stopping poi'/er by changing the value of , the average 
excitation energy. 

✓ 

As far as chemical combination in gases is concerned, the effect. 

If present, must be ouite small. If the atoms act independently, the 
stopping povver per molecule should be a simple sum of the atomic 
stopping powers of the individual atoms. This rule seems valid in 
gases though not always in liquids and solids, (See the excellent 
revieivs of this aspect of stopping pov/er given by Oeiger (1927) and 
Taylor (1952),) v7e aro concerned with the stopping power of com- 
pounds in two respects. First, seme of the impurities found in the 
gases are compounded of more than one element; second, one of the 
primary gases studied in this research (nitrogen) occurs in a diatomic 



The relation betv/een range and stopuing pov/er is easily seen to be 



R(E„) : 






Gtate 



Theory hints that tiore appreciable effects night be forth- 
coning if the stopping nediim is in a condensed state, because of 
the influence of the nany adjacent atons on the outer electron energy 
levels. Sone early experiments seemed to indicate a consistent 
difference between the stopping of alpha particles by elenental 
solids and gases (Brings, 1927; Geiger, 1927). These seeming 
differences v/ere later shovm to be due to experimental factors 
(Livingston and Bethe, 1937). Experimental data - not entirely in 
agreement among themselves - still accumulate however to indicate 
that special effects raay be present in certain coses (Taylor, 1952). 

Both theory and experim.ent have shoim i,t rked effects in 
condensed substances at relativistic velocities of incident pf-rtloles, 
but v/e are not interested in that energy region. 

Theoreticians have given special attention to stopning by metals. 
The outer electrons in metals are practically free and their behavior 
is radically different from those electrons bound to individual atoms. 

As early as 1933, von V/eisacker (1933) derived a formula which strongly 
related stopping pov;er to conductivity. His theory was sho\'m by 
subsequent experiment (Gerritsen, 1946) to make incorrect predictions, 
and Kramers (1947) proved theoretically that special conduction 
effects v/ere practically cancelled by polarization effects. A. Bohr 
(1948) and Pines (1952) published further estimates of small reduc- 
tions in stopping pov/er by free valence electrons for incident particles 
of moderate velocity. These corrective factors are not large, however, 
and for heavier elements v/herein the proportion of conduction elec- 






irons is small the corrections predicted are auite negligible. 

This is confirmed experimentally by Heller and Tendon (1951), v/ho 
found that for Z ^ 13 no consistent differences between stopping 
povrer of metals and semi-condhctors are found except the expected 
variations v/ith atomic number. (Note should be made that such 
results v/ere obtained only for particle velocities in excess of the 
velocity of 1 Mev protons.) The only experimental confirmation of 
special effects due to conduction properties of metals is the work 
of Hadsen and Venkateswarlu (1948a) on stooping by beryllium for 
protons in the energy region from 500 to 1500 kev. The effect is 
expressed as an abnormally large value of I for beryllium to be 

8, V 

used in the stopping pov/er equation. 

A. Bohr further indicated possible polarization effects may 
affect stopping povrer for condensed media even v;hen non-conducting. 
Halpern and Hall (1948) also suggested the possibility that the more 
loosely boimd electrons in condensed media nay have different oscil- 
lator frequencies associated with them than they would have in the 
gaseous state. The special effects predicted are still small except 
in those cases where the outer electrons constitute a large pro- 
portion of all electrons effective in stopping. 

7. Capture and Loss of Electrons. 

In 1922, Henderson (1922b) observed that a beam of alpha p-irticles 
contains a portion of singly charged particles, which portion increases 
with decreasing bean velocity. Since that time experiment and theory 
have confirmed that for all charged particles there is a pronounced 
tendency at low velocities for the particle to change its charge by 



capt\ire and loss of electrons. It is desirable to consider the 
X>henaEienon briefly in respect to stopping power of protons, since 
the stopping power theory presented above depends on the retention 
of positive charge by the proton. 

Theory (Bohr, 1948) has given an approximate answer to the 
cross-section for capture of on electron by a proton, giving neutral 
hydrogen; and also a formula is provided for the cross-section for 
loss of this electron by the neutral hydrogen in flight. A proton 
undergoes in its passage through matter a succession of neutralizing 
and re-ionizing events, so that the proportion of its flight which 
occurs in the charged state may be obtained by taking the ratio of 
oross-sections. 

Let y z the proportion of flight as neutral hydrogen; 






(I-B-90) 



where 

0 ”. - the electron capture cross-section ; 

c ^ 

I 

(7“^ s the electron loss cross-section . 



If , then 






y ^ (T/<ri 




(I-B-91) 


According to Bohr (1948) , 
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Since z eouals unity, 


y = 4 Z"^/^ • 




(I-B-94) 


irow in our experimental work the lovvost value of v 


is 0.81 X 10 



cm/sec (corresponding to an energy of 340 kev) , the lor/est value of Z 
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f 



is 7, and the value of t is 0,219 x 10^ cm/sec , so that y equals 

o max 



0,0031 , Thus y is so low and decreases so rapidly as particle 
Velocity increases that we may ignore entirely the possibility of 
electron capture, 

A confirmatory note should be added, T, Hall (see Warshaw, 1949) 
has found that experimental results are in substantial agreement with 
Bohr's theory used above, 

B, Straggling In Stopping power, 

Bohr's equation for straggling in stopping pot/er is given by 
substituting equation (I-B-25) Into eqviation (l-B-23): 



•niiB equation is valid for fast incident particles. For slower 
particles, such that some atomic electrons are no longer slow compared 
to the particle velocity, Bohr returns to the Thomas-Fenni model and 
uses only those electrons slovier than 2v . That is, Z is replaced 
by Z * , where 



os determined from equation (I-B-72), 

Livingston and Bethe (1937) outline an apnroach to this prob- 
lem which parallels the analysis of Bohr, previously given, in much the 
same way that Bethe's derivation of stopping poi?er formula is analo- 
gous to Bohr's classical apnroach. They obtain: 



4- (il^) = 4lrz2e^NZ . 



(I-B-95) 



= (2v/Tq) , 



(I-B-96) 




(I-B-97) 



frtiere Z* 



is defined as in equation (I-B-71); k is a nvimerlcal 



n 



factor usxially equal to about 4/3 ; is the number of electrons 
in the n^ shell; and P is the average excitation energy of the 
n*^ shell, as usual. 

This formula is somewhat difficult to use, as the values of 
can only be estimated. Purthemore the value given is only 

approximate and its true value is difficult to determine. 

The extra factor in the bracket is the essential correction to 
Bohr*s foTBiula, and actually seems to be not so much a true quantum 
correction, as a term to io take into account the motion of the 
orbital electrons. (The factor in the numerator is actually 

only an approximation for the kinetic energy of the electron.) 

Bohr*s theory, it must be remembered, assimes that the electrons 
are motionless before collision. Averaged over all particle- 
electron collisions, this is approximately true, so that the stopping 
power relation is not affected by such an assumption; but energy 
transfer variations in individual collision cases are enlarged be- 
cause of this electron motion, so that the stragrrllng is Increased. 

C. Review of Existing Experimental Data. 

1. Stopping Power. 

To review all experimental data which has ever been published 
on stopping power would be very tedious and quite unnecessary. At 
the opposite extreme, to list available experimental results only if 
they cover precisely tho same set of conditions and materials v.'hich 
we ere studying would be rather simple - there is very little to date. 

It is profitable to take a middle course - to list not only that 
data which may duplicate ovu: oivn, but also to list data in either one 



of the following categoriee: 

(a) stopping power for protons of the raoteriels v;e are study- 
ing in energy regions to either side of the regions v;e are 
ooncerned V7ith (400 - 1100 kev) ; 

(b) stopping paver for protons of other elements, in the energy 
region which coincides v/ith ours. 

In this vxay we may synthesize our results with others into the 
overall stopping power scheme. 



Table II. 



Experimentally Determined Values of 
Average Excitation Energy 



Element 


2 


2* (if 
tised) 


I 

av 

(ev) 


Particle and 
velocity or 
energy 


Reference 




V 


81 




♦ 


rano, 1934 


lie 


10 


132 




* 


ti If 


A 


18 


195 




* 


ft ft 


Hi 


28 


325 




* 


ft If 


Cu 


29 


320 




* 


If fi 






279 


300 


Mev prot* 


Bakker & Segrh, 1951 






307 


240 


ft ft 


Hather & Segri, 1951 






375 


70 


ft ft 


Bloerabergen & van 






365 


100 


ft ft 


Heerden, 1951 






27.4 


293 


2 


tt tf 


Liv. ?c Bethe, 1937 


Er 


36 


390 






Mano, 1934 


Xe 


54 


530 




* 


It IT 



Alpha Particles from natural radioactive elements (6-11 I.'ev) , 
which corresponds to protons having velocities equivalent to 
1,5 - 2.8 Mev of energy. 
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Let first list all data available for energies higher than 

the range studied herein, 'Hiis is best suimnarized by indicating the 

value of I found necessary to rnake Bethels formula (eauation I-B-69) 
av 

cheek with experiment (see Table II, above). 

On the low energy side of the energy spectrum, some data is 
available relating to the materials in which v/e are interested. A 
eeirtain amount of work has been done by the "Kevatron" group at 
Chicago. V/llcor (1948) published results at lovr proton energies for 
aluminum and gold, but his results have been superseded by better data, 
itfarshaw (1949) has published data for protons of energies from 40 to 
360 kev for aluminum, copper, silver, and gold. He extrapolated his 
curves up to several Hev, where the theoretical formula of Betho would 
be expected to be valid. 

Work for gases has been reported recently both from the Univer- 
sity of Chicago and the California Institute of Technology. From the 
latter institution Dunbar, et al. (1952) have reported on stopping 
power of some gases for protons of energies from 25 to 600 kev. 
Additional data has been recently provided by this group of workers^ 
which means that most of the simple gases, including all those we 
are reporting on herein, hove been covered for protons up to 600 kev 
of energy. Their curves ere rejjroduced in Figure 13, to show the 
match with our data. Weyl (1953) at Chicago has jrovided data for 
stopping pov/er of various light gases, including argon, for protons of 
energies from 50 to 500 kev. His data is consistent with that from 

^ Up-to-date information from the Cal. Tech, laboratory has been kindly 
transmitted by Prof. V;. V/haling (personal communication). This is 
expected to be published later under authorship of Dunbar et al. 
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the California Institute of Technology and is therefore not plotted. ^ 

Stopping data prcrviding a complete coverage of the particular 

region of interest to us has been provided by only a fev workers. 

Madsen and Venkateswarlu (1948a), being the first to employ the 

experimental technique we employ with the Van de Graaff generator, 

published results for berylllxm stopping of jarotons with energies 

from 500 to 1500 kev. Using Bethe’s formula, and Including the 

correction term C , they found the data to be consistent v»ith a 
K 

value of I of 64 ± 5 ev. 
av 

Using the same technique Huus and Madsen (1949) obtained a 
little data In this region for gold. Madsen has also recently sub- 
mitted seme data for copper.^ His results for the latter element are 
not consistent vrlth other eocperlmental work, which is possibly due t© 
difficulties In obtaining foils of consistent thickness. 

The greatest amo\mt of v/ork in this region Is provided by 
Kahn (1953), v/ho has obtained cxirves for protons of energies from 
400 kev to 2,0 Mev for beryllium, aluminum, copper, and gold. 

Figure 12 contains his results, which agrees T'/lth the earlier data on 
beryllium and gold mentioned above. 

2. Sferaggling. 

Very little previous experimental work has been published on 
straggling In stopping power In the energy region of interest to us. 
The reason is easily seen. Even the techniques for stop-oing power 

^ Detailed Information on the v/ork at Chicago, os v/ell ns other 
information available to him, has been kindly provided by Pfof. 

S. K. Allison (personal comtiunicatlon) . 

2 

I7e acknowledge appreciatively a personal communication from Dr. 
Madsen, giving detailed information on all hie v/ork to date. 

- - 



determination in this energy region are sanevfhat crude, and eruerl- 
nental variations scmetiiaes wide. It can be appreciated that the 
theoretically valid deviations in stopping ijower con be easily masked 
by the same factors which cause inaccuracies in stopping newer ex- 
perimeots. 

The only experimental v;ork published on straggling is by Madsen 
and Venkateswarlu (1948b), vdio were able to analyze the results of 
their stopping povrer data (1948a) to give straggling in addition to 
mean stopping povrer. Results are for beryllium, V7ith proton energies 
from 400 to 1200 kev. The results vary v7idoly, but seem on the average 
to be in reasonable agreement with the theory of Bohr (1948) , given 
herein as equation (I-B-95), The value of il/(Ax)^ v;as found to 
average 8,99 kev-cm/mg^ , with individual readings ranging from 6.9 
to 10.4 . In these units, the theoretical prediction is 8.3 . 



S' 



II. Bostc Eroerlmental TechnlQU®. 



A. Basic Principle. 

If a beam of monoenergetic pajotons impinges upon a target of 
a certain material, a nuclear reaction may be induced which pro- 
duces quanta or particles as a result. The investigation of proton- 
garnna reactions requires the determination of the number of gamma 
photons produced (and detected) per designated number of protons 
incident on the target material. For a given target, such a figure, 
which we may call the yield , is a function of the energy of the in- 
cident proton. By varying the proton energy in small Increments and 
measuring the yield at each step, one can accumulate data on yield 
versus proton energy. For mony light isotopes as targets, a plot of 
these results in the Van de Graaff energy region (a tm hundred kev 
up to several I'ev) gives a type of curve which is low and slovjly 
rising, on v;hich is superposed a n\mber of "resonance ^^eaks" whose 
height, width, and position in the energy spectrum are characteristic 
of the isotope (s) involved. 

IIo;v if a given thickness of matter is interposed in the noth of 
the beam after it leaves the accelerator but before it hits the target, 
the protons are slowed by amounts, the average of which is determined 
by the overage energy loss of the protons in the interposed material. 
Then to duplicate a "peak" in a yield-versus-proton energy experiment, 
the protons obviously must be given an initial energy greater than that 
characteristic of the target material by just the amount of energy loss 






in the interposed noterial. 







m 




I 



If, in such ar experlnent, one ineaGUires proton energy before 
the beam strikes the interposed stopping material, a yield-energy 
curve nay be constructed, called hereafter a shifted or displaced 
spectrum. The peak-to-peak correspondence vfith the imdeviatefl spectrum 
is easily established by inspection. 

The undeviated peak energy values may bo determined by experi- 
ments without stopping media interposed or, if already xvell estab- 
lished, they may be determined from the literature. One then can 
subtract the energy values for the corresponding peaks in the 
deviated and undevlated spectrum to obtain energy loss 4 E iu the 
material. Dividing by the path length of the protons in the medium 
gives ^E/Ax . The theorem of mean value indicates that this v^lue, 
which depends on the energy of the beam as it enters and as it leaves 
the medium, is the same as the value of dE/dx for some value 'of E 
between the Incident and exit beam energy. This mean value of E is 
not difficult to determine when AE is small compared to E . 

It is found that the peaks in the displaced spectrum are broader 
than those in the basic, unshlfted spectmira. This is caused by 
individual proton variation, or straggling, in stopping pcr./er. Ihat 
is, when deviations in energy loss are possible for individual protons 
passing through matter, then beams of average energy aoprecinbly 
different from the proper displaced peak value will still contain 
protons of the proper energy to excite the target atoms In the 
resonance region. 

The technique briefly outlinei/obove v;as first used by JTadsen and 
Venkateswarlu (1948a, b). A more detailed analysis is noiv presented. 
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showing the corrections neceos'^ry to get r. tr';e onsvor fron oxocri- 



nental results, 

B, Detailed Analysis. 
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Figure 3 is a Bhhematic diagram of the atmaratus for measure- 
ment of stopping pa.7er, with the voltages specified for a proton 
from the moment it leaves the accelerator section until it enters 
^the target nucleus to release a phtton to be counted. In this 
figure. 



s energy of a phot on- producing proton as it leaves the 
accelerator; 

Ej^ = energy of the photon- producing proton after nassina 
through the analyzing field; 

E s energy of the photon-producing proton after pas'^ing the 
defining slit and as it enters the foil; 

E„ e energy of the photon-producing proton as it leaves the 
foil and enters the gas space; 
r energy of the photon- producing nroton os it leaver tlic 
gas space and enters the target material; 



E_ = energy of the photon-producing proton as It enters the 
o 

nucleus to produce the reaction. 

V.'e note that if the gas chamber is eracuated, . 

In order to find the energy Eq of a aarrfia-producing proton os 
it comes frcun the accelesrator, it is possible to start v;ith Eg and 
successively add to it the energy losses suffered in all parts of the 

system between E_ and E . Thus 

o o 

= %*l'iEjMAE^|+lAEyl+|AE^l+(AE^ll . (II-B-1) 

vjhero 



/lE^ = the energy change of the proton In passing through a 
certain amount of the target material; 

AE r the energy change of the proton In passing through 
G 

the gas; 

AEp s the energy claange of the proton In passing through 
the foil; 

AE s the energy change of the pooton in passing through the 

s 

defining slit, or as a result of the regulatory action 
of the defining slit; 

the energy change of the proton in passing through the 
magnet, or as a result of any niscellaneoios effect In 
the system. 

V7hat v;e are particularly concerned \vlth, however, are the 
relationships among the distributions in energy of the probability of 
detecting a gama photon for each proton accepted from the accelerator. 
These distribution functions depend not only upon proton energy, but 
also upon the position in the apparatus where the energy is specified, 






I 






so that we conclude that there is a probability distribution for each 
of the variables in equation (II-B-1). 

The statistical variables of Interest are the nean and the 
variance (sqxmre of the standard deviation) of the distributions. For 
the mean value: 

= *5*|4%i+|^l+fe( + |A%|+(^,| : (II-B-2) 

and for the variances, the following formula is valid (Hoel, 1947) : 



_2 ^2 ^2 , ^2 ^2 ^2 «2 

_Q = +_n +/I +/1 +jl /it n 3\ 

0 5 T G S F M » tII-B-3) 

where the subscripts refer to the same portion of the apparatus as 

in equation (II-B-1). The precise values of sane of these statis- 
tical variables v/ill be determined later, but others can be simplified 
in such an obvious v;ay that we may convlently do so at this point. 

We assume that the resonance peaks are so ivldely separated that 
they do not interfere v/ith one another, and thus the statistical 
analysis may be carried out on each peak alone. Furthermore, the 
frequency distribution is considered sufficiently symmetrical that the 
maximum value can be taken as the mean value without appreciable error. 
To simplify the notation let us assume, unless otherwise stated, that 
when we refer to energy henceforth, we mean the mean value of the 
energy of the whole proton bean, without the necessity of placement of 
a bar over the symbol. 

Eq becomes the average beam energy from the accelerator giving 
the resonance maximum. Eg is the value of proton energy as it enters 
the target nucleus to give the resonance maximum, and we will therefore 

AE^ is the average value of energy loss of 



re-designate it as E, 



protons in the target material, and is therefore ^ , where T is 
the target thickness expressed es the average energy lose for 
protons passing through the target. If we assume that the effect 
on individual protons by the slit and the magnet are os likely to be 
in a positive as a negative direction (see discussion below) then 
Aig and are zero, 

JTq is the standard deviation of the yield curve about the 
resonance maximum point, and is determined from the experimental data; 
whereas ilg is the standard deviation of the natural, or theoretical 
curve and will be called jjg-t henceforth, and mensure 

the straggling effect in energy loss in the gas and foil respectively. 
The deviations Caused by the slit and magnet may be explained somewhat 

(a) The slit is of finite vd.dth, and the voltage regulation system 
of the Van de Graaff generator (Smith, 1952) is such as to keep the 
I^oton beam betv/een the sides of the slit. Then the maximum swing 
in beam voltage is on the order of the amount required to siving the 
beam from one side of the slit to the other. If the value of average 
beam voltage is measured v^ith the beam centered in the slit, the 
deviations are equally likely to be plus or minus. 

(b) VJe may expect other miscellaneous variations of a presumed random 

nature: inhomogeneities in the magnetic field, variations in the 

vertical and horizontal position of the principal axis of the focused 

beam, slight variations in the electranagnetic field current, and 

» 

miscellaneous variations of a human or instrumental nature. These are 
conveniently grouped together with the subscript . 

Using the modifications given aJ?ove v;e can now re-write equations 



(II-B-2) and (II-B-3). Let us first do so for an experimental 
curve XTlthout any stopping media at all: 






o 

^2 -2 ^2 ^2 
n =il +il 

o nat T S M 



(II-B-4) 



. (II-B-5) 

Then we can v.Tlte the eouations valid for experimental curves where 
stopping media are interposed in the beam. In doing so those values 
which are the seme for the corresponding peak without media inter- 
posed are left unprimed, values which are shifted are primed: 



and 



,2 2 _2 2 ^,2 ,2 

* n + +A, 

*o nat T S H 



(II-B-7) 

In the latter formula, signifies the sum of the variations in 

energy loss through both foil and gas. In practice it will denote 
only the variation for the foil, inasmuch as we have obtained no 
resists in cases where both gas and foil are used. 

Subtracting equations (II-B-4) from (II-B-6), and ( II-B-5) from 
(II-B-7) gives the basic eouations which permit determination of 
stopping power and straggling from experimental results for the 
resonance peaks: 



V/ithout gas. 



A'E 



E* 

o 



vlth gas and foil. 
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(lI-B-8) 



(II-B-9) 



and without gas, 
2 



#2 



ji = (/i; - /!,) - (ii; . Ag) - ui;, - i2„) _ 

(II-B-IO) 

That the primed and unprimed values of Jig and do not alrrays 

cancel may be understood from the discussion ftn sections in.F.3 and 4, 

-S'y- 



III» Erperinental Details and Preliminary Calculations. 



The principle of operation of the Van de Graaff type of 
generator-accelerator is too well knovm to require any explanation 
here, A description of the particular one at the Ohio state 
UnlTersity and its associated oxixillary equipment may be found 
elsewhere (Grove, 1947; Cooper, 1949; Grove, 1950; Taylor, 1952; 

Smith, 1952). Certain improvements and modifications ware 
required in order to fit the machine to the purpose of the present 
research. They will be described in the sections immediately following. 
A, Target-Faraday Cage Assonbly. 

In order that the protons may be collected and counted, the 
target on vihich the protons impinge is mounted on an assembly called 
a "Faraday cage," This is Insulated electrically from the rest of 
the accelerator, so that the collected positive charge can drain to 
ground through an "integrator" v/hich measures this charge. 

It Was necessary to provide a mechanism for introducing a foil 
into the proton beam; also a gas-holding chamber was required with a 
thin foil-covered window so that protons could be introduced to 
traverse the gas. For this purpose a type of valve box^ was 
modified to perform all the functions of target holder, Faraday cage, 
gas holder, and foil retainer. This assembly is depicted in Figure 
4, and is almost self-explanatory. 

The system may operate with the valve down for basic calibration 
runs. If stopping power of a foil is to be investigated, the foil 

^ This box was originally mode after a type used on the University of 
Wisconsin Van de Graaff generator ^tfie pattern for which was 
generously provided by Prof. R, G, Herb . 
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Is attached (with the use of clear glyptal) across the slit in the 
Talve, which is then raised into a vertical position. In this position 
the system is aligned so that protons passing first through the beam 
defining slit strike the foil at approximately the middle of the 
valve slit. 

If stopping iKraer of a gas is reoulred, the foil across the 
valve slit is made gas-tight by painting with glyptal around all the 
foil edges. The system is placed on the vacuum system and the valve 
is raised end clamped against the valve seat to gas-tightness, ("-the 
valve may be manipulated from the outside by a shaft through a Wilson 
seal.) Gas may then be introduced to the desired pressure, and the 
experiment performed. This technique Interposes both a foil and the 
gas in the path of the beam at the same time, so that in determining 
the stopping power of the gas the energy loss in the foil must be 

i 

known and allowed for. 

B. Preliminary Measurements on Foi3s . 

The foils used i7ere of copper and nickel, of norainal thickness 
betv;een 0.03 mil and 0.1 mil. This amounts to a surface density of 
from 0,6 to 2 mg/cm^. This is sonep/hat thicker than is usually used 
for stopping power experiments; but in selecting such foils the 
greatest Importance was placed on consistency in thickness - the 
lack of which has often plagued experimenters viio used thinner foils. 
Furthermore certain fixed errors, such as that caused by carbon 
accumulation on the foil, are less percentagevjlse When thicker 
foils are used. (See Section IV. A, below.) 

The foils vjere obtained from the Chromium Corporation of 

- 60 - 
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>\merlca, ?'aterbury, Conn., who prepared them by electrolysis. 

They were guaranteed by the coupany to have an average thickness with- 
in lOjT. of the nominal thickness, and to have a spot-to-spot variation 
in thickness of less than 5^. The foils used vjere judiciously selec- 
ted to obtain the best consistency possible; and a technlaue was 
developed to check thickness variations to v/lthin about 1^?. 

This technique depends on the range-straggling of alpha particles. 
The apparatus is quite simple and is Illustrated in Figure 5. The 
principle of operation is explained by following the procedure out- 
lined below. 

Vfith the foil removed, data is taken of counting rate as a 
fvmction of height of the movable platform from the fixed surface. 

This gives a cvirve similar to curve (a). Figure 6. The steepness of 
rise of the curve is determined by the range-stragPiling of the ali^a 
particles and is important in determination of the sensitivity of 
measurement of the relative foil thickness. Now if the experiment 
is repeated v/ith the foil interposed over the 1/32" hole, a curve of 
similar shape, but displaced by amount A, is able to be detemiined. 
This distance A is the air-equivalent of the foil for alpha 
particle stopping. 

If we set the apparatus at height H, to give an operating point 
on a reasonably linear part of the steep portion of the curve, we may 
then take readings with various portions of the foil positioned over 
the hole. If the counting rate changes by amount Ac, then this 

^ Information on khese foils was kindly provided by Itr. IJ, M. Stem- 
fels, Chemical Engineer for the Chromium Corp, 
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APPARATUS FOR CHECKING 
THICKNESS VARIATION IN 
THIN FOILS 



THIN WINDOW (^1.5MG/CM^) 

GEIGER TUBE , SUCH AS 
TRACERLAB - TG C - 2 
-HEIGHT FIXED- 



THIS APPARATUS 
MAY BE RAISED OR 
LOWERED BY TURN- 
I NO SCREWS 



3 ^^ 



FOIL TO BE STUDIED 
OVER 1/32" HOLE 

—THIN DEPOSIT OF 
POLONIUM (~IOpC) 



FIXED SURFACE 
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FOR CHECKING VARIATIONS IN 
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indicated a change in air-equlOTlent of Ah . Ihe value of Ah/A 
indicates the fractional variation in thickness frm the original 
spot to the nev/ spot studied. If a large number of spots on the 
foil are studied thus, T/e can determine the mean value of H by 
knowing the overall average counting rate. The statistical variation 
in thickness is easily computed as the average root-mean-sarare of 
^h/A • Such a variation corresponds to a counting variation 
(fractional) of 2ie/C . To translate a counting percentage devia- 
tion to a thickness percentage deviation, we obviously have to mul- 
tiply the former by 

PBCtor r - , « (m-A-D 

A Ac A ‘(slope cvirve at opert’g pt.) 

If a finite number of counts are taken at each point, the count- 
ing rate is subject to a statistical variation due to the randomness 
of the radioactive emission process. Such statistical variation must 
be subtracted from the total variation to give the most probable 
variation due to thickness variations in the foil. Furthermore the 
statistical variation from the radioactive emission -orocesa nlaces 
a practical limit on the accuracy of this technique. 

Table IV (a) gives a set of sample data taken on one of the foils; 
and Table IV (b) summarizes the results of measio'ement on all the 
foils. 

Y.q conclude from the data that the standard deviation of the 
thickness measured over spots 1/32" in diameter is about the order of 
magnitude of the precision of the measiirements. V/e estimate that we 
can cite for all the foils a standard deviation of about 1.45; as 
measiire of the spot-to-spot variation from the mean. Thus, if we 






Table IV 



Spot-to-Spot Variations in Foil Thicknesses 
Spot Size : 1/32 " 

Part (a) ; Sample Data - 0,1 nil Copper Foil. 



Spot No, 


Counts 


Sec. 


Cts./nln, 


Ac 


(Ac)2 


1 


4096 


828,5 


297 


- 1.4 


1.96 


2 


rt 


770.3 


319 


+ 20.6 


424.36 


3 


ft 


800.0 


307 


+ 8.6 


73.96 


4 


ft 


831.0 


296 


- 2.4 


5.76 


5 


ft 


844.0 


291 


- 7.4 


54.76 


6 


ft 


827.0 


297 


- 1.4 


1,96 


7 


ft 


842,0 


292 


- 6.4 


40.96 


8 


tf 


853.0 


288 


-10.4 


108.16 






II 

• 

< 


290,4 




711.88 



Standard Deviation (small sample theory) 



711,88 ' 



Standard 



« 10.1 , or 3.4f, 

Deviation of Radioact, Decay Rate = 



= 1 . 56 /. 




C 

A 

Slope 



298,4 Cts,/mln, 
3,0 turns 
100 Cts./turn 



Conversion Factor » 

8 X 100 

• 0.37 



Part (b) ; Sirnmarj’ of Results. 



Foil 


No, of 

Spots 

I'easured 


Total 

Std. 

Dev. 


Std. 

Dev, of 
Radioac. 


Net Stt, 
Dev. in 
Count .Rate 


Conver- 

sion 

Factor 


Sti, 

Dev. of 
Thickn. 


,05nil Ni-1 


9 


5,e^o 


2,2^0 




.54 


2.9/- 


tt 


11 


3.7 


2.2 


3.0 


,50 


1.5 


rf 


13 


3,5 


2.2 


2.7 


.53 


1.4 


,05nil Ni-2 


10 


3.8 


2.2 


3.0 


.48 


1.4 


.04mil Ni 


5 


3.5 


2.2 


2.7 


.67 


1.8 


ft 


9 


2.9 


2.2 


1.9 


.71 


1.4 


,03nll Ki 


8 


2.6 


2.2 


1.4 


.83 


1.2 


,05nil Cu 


5 


1.3 


1.7 


0 


- 


small 


,1 nil Cu 


8 


3.4 


1.6 


3,0 


.37 


1.1 
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use average thickness for our canputations, while the proton beam 
passes through the foil at a spot selected at random, a probable error 
In stopping power results of about 1^ •xists because of thickness 
variation alone. 

We can drav/ the above conclusion because the 1/32” spot size 
used is about the same os the diameter of the intense "core" of the 
proton beam fraa the Van de Granff generator when it is well focused. 

We cannot as easily make the same conclusions with regard to strag- 
gling measurements, since variations of thickness v/ithin a slnrrle 
1/32" spot will Increase the straggling without affecting average 
energy loss. 

It is possible to make a hypothesis that the variations using a 
l/32" spot are due to more radical variations vrithin the snot of this 
size which ore not completely averaged out. This hypothesis is opnosed 
to the hypothesis that the variations are broad and slowly changing 
in comparison to the 1/32" distance taken as a unit. These vjere 
tested by making a series of runs in one of the nickel foils v/lth 
successive spots btudied being contiguous and along a straight line. 

The results v;ere quite consistent v;ith the hypothesis that the varia- 
tions in foil thickness ivere quite gradual, and highly inconsistent 
with the hypothesis of more rapid, stronger variations. 

This still does not obviate the possibility that very finely 
spaced variations nay exist which average out almost comnlrtely over 

a 1/32" spot. The possible existence of such variations therefore for 

-6 -2 

spot checks of size 10 to 10 centimeters must still be admitted. 
Such variations may cause errors in the exnerimental results on 
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stopping pov/er straggling and are discussed in the section concerned 
with such errors (Section IV,D.2) , 

The average thicknesses of the foils, measiared in units of 
were determined by weighing them and measuring their linear dimensions. 
Tbe iTeighing was accomplished by microbalance in the liicroanalytical 
Laboratory, Chemistry Department, Ohio State University, through 
courtesy of Professor W. M, MacNevin. 1116 measurement of the linear 
dimensions was accomplished under a lov/ power microscope, and checked 
by steel rule end reading glass. The foils v;ere rectangular in shape, 
as checked by assuring oneself of the eouality of the diagonal dimen- 
sions, and therefore the area was taken to be simply the nroduct of 
the two linear dimensions. Results of the foil measuirements are -Iven 
in Ifable V. 



Table V. 



Determination of Foil Thicknesses 



Foil 


Area (cm^) 


V'eight (mg) 


Thickness 

(rng/en?) 


,05 mil Nl-1 


7,05 ±:.03 (. 435 J) 


8.10 i.Ol (.184?;) 


1.150 Jt(.45f,) 


,05 mil Ml -2 


3.28 dt .02 


3.86 .01 (.26?;) 


1.175 -t(.65fj) 


»04 mil N1 


7.40 ^ .03 ( .4(P/) 


7.23 i-- .01 (.13^;) 


0.977 i(.45f4 


,03 mil Ni 


7.40 ±. .03 (,40r.) 


4.75 J: .01 (.21?;) 


0.640 ±(.45?;) 


,05 mil Cu 


7.495 02 (.27?;) 


8.65 .01 (.llGf^) 


1.155 +(.35/) 


.1 mil Cu 


7o42 i .03 (.40?J) 


16.30 .01 (.061^.) 


2.20 +(.45?;) 



Ue may conclude that the probable errors inherent in the use of 

the thicknesses given are the statistical sum of the 1^ due to non- 

uniformity and 0,45^ due to error in measurement of average thickness. 

This gives a total probable error of about l.lf> in estimating average 

- 



proton path length through the foils by the above values of foil 
thickness. 

C« Design of Magnet Current Regulator, 

It v®G found upon first operating the Van de Graaff generator 
that the current provided to the analysing electromagnet was not 
sufficiently stable for the close energy resolution desired. It 
xras found that for a given setting of the controls, the magnet 
current tended to sv/ing to either side by a variable amount, giving 
a deviation from the mean of as much as 0 , 15 ? , 

This amount can easily be translated to energy variation. If 
we accept a rough formula (see Section III.F, belcn?) : 

E = (4/3) I^ , (III-C-1) 

It is easily sham that a deviation of O.Lfo in the current, I, is 
equivalent to a deviation of 0,2f? in the beam energy, E, Thus at 
beam voltage of 1 Mev, it was iinoossible to prevent the voltage from 
wandering as much ns 2 kev to either side of the mean. 

It Is true that the operator has manual controls for setting 
the magnet cxurent; but the variation v/as so erratic and sudden at 
tines that the operator could not anticipate and counteract the more 
rapid svjings. The disadvantages inherent in this variation are manifest: 

(a) 'Rxey make it difficult to maintain and read accurately the 
mean setting of the electromagnet current, thiis preventing accurate 
calibration of the beam energy in terms of magnet current. 

(b) Since beam position at the slit system (see Figure S) , located 

after the magnetic analyser , is used to stabilize the accelerator 

Voltage (Smith, 1962) , magnetic current variations are reflected 
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in beam energy variationa. The beam stabilizer, tha‘k,is> can never 
present beam energy variations corresponding to magnet current 
variations* 

In terms of the variables v;e are trying to determine, this var- 
iation would mean an uncertainty in peak position of about 1 kev 
added to all other imcertaintiee and would add greatly to the fpctoro 
constituting theiljj discussed in Section II. B above. It v;as found 
expedient therefore to provide some measure of regulation of the 
magnet current. 

A brief discussion of how the magnet cvirrent is provided would be 
appropriate. A raotolJ-generator set puts out a d.o. sunnly of maximum 
rated values 360 volts and 4.75 amperes. The output voltage is 
controlled by varying the d.c. voltage to a separately excited field 
winding on the generator side. The d.c, voltage for the field is 
obtained by full-rrave rectification of a variac-controlled alternating 
voltage, obtainid in turn frcan a Sola-regulated 115-volt a.c. power 
supply. Thus the magnet cmrrent can be varied by changing the 
setting of the variac supplying the field current. This actxially 
varies the voltage from the d.c, generator; a fine control of current 
is obtained by means of a variable rheostat in series with the magnet 
coils. 

A capacitance of 250^fd. is connected in parallel »;ith the 
magnet. Since the magnet coils have a high inductive reactance 
(inductance estimated to be about 30 henries), fluctuations of 
frequency much faster than one cycle per second pass through the 
condensers rather than the coll. Aside from this, and the Sola 
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regvilation mentioned above, no other regulation was in the circuit 
prior to the time the regulator herein described was installed. 

Several possible schemes for regulation of current itself were 
considered and rejected. It was decided that uncontrollable resis- 
tance changes in the magnet circuit (caused by temperature changes in 
the coils) were slow enough to be compensated for manually by the 
operator. p\»rthermore fine control of current by a variable resistor 
was desired. Thus, a voltage regulator was actually used for stabil- 
izing the current. The d.c. generator voltage is affected not only 
by the field current, but also by the motor speed, v.hich may nossibly 
experience small variations. It was decided therefore to nut the regu- 
lator directly in the magnet circuit betvjeen the generator and the 
load, which consists of magnet coils and the c<aitrol rheostats. 

It was necessary under the circumstances to keep the system 
quite simple, since any complex electronic apparatus has a very low 
power efficiency and the output of the d.c. generator is rather limited. 
The very simplest type of regulator of the "transconductance" type is 
used (Hunt and Hickman, 1939) . 




Figure 7 shows in ehheraatic fashion the regulator In principle. 
The triode is in parallel with the load, and the current through it 
adds to the magnet current to give the total current required from 



the d,c, source. The amount of tube current depends on the output 
voltage and the grid voltage. If a small positive variation occurs 
In the input voltage, more current flows fran the source. This rise 
also causes a rise in grid voltage, v/hlch produces a greater current 
drain through the tube. With proper choice of the circuit parameters, 
the current increase through the tube matches the total current rise 
from the source - causing the same amount of current through the loed. 
Or, from another point of view, a current rise through the tube 
creates additional voltage drop across , which is designed to 
just counteract the original voltage rise at the input side - thus 
giving the same output voltage, 

A rigorous proof of the above statements, indicating how the 
circuit parameters are selected is given in the apnendix, 7\lso in 
the appendix is a circuit diagram of the complete magnet current 
control and regulating system, ns it exists at the present time, with 
further explanatory comments on its operation. 

Results v;ith the regulator have been excellent. With proner 
attention on the part of the operator to compensate for slow 
resistance changes and keep the regulator controls properly adjusted, 
the current variations are cut do\m to about one- tenth their value 
without regulation. This means a voltage variation of the bean of 
about 100 - 200 ov , which is sroall compared to other variations and 
errors . 

D, Technique for Taking and Correcting Yield-Energy Data. 

I’flth the chosen target in place and a vacuum on the whole system, 

the Van de Graaff generator is started and protons produced in 

- 7 / 



accordance with standard procedure. The iiiagnet Is then "cycled” 
twice by running the magnet current up to a relatively high value 
(about 3 amperes) and then back to zero. The subsequent data is 
taken v;ith ever increasing magnet current, so that one stays consis- 
tently on a reproducible magnetization curve (magnetic flux versus 
magnet current), v/ithout any hysteresis effects. 

Ab and if required, the stopping foil is inserted and clamped 
at this point in the japocedure, and the desired amotint of gas (vdiich 
may be zero) is introduced into the holder (Figure 4.) . 

V/ith the magnet current set at the Initial (lov/est) value, the 
accelerator voltage is set so as to send the protons through the 
defining slits. The proton accumulated charge is counted by on 
•integrator" with a neehanlcal register and the gamma-rays detected 
by the usual Geiger counter-scaler apnaratus. At each current set- 
ting the data taken includes: magnet c\irrent, total gamma-ray counts 

(usually set to a predetermined figure) , time to collect these 
counts, total integrator count of accumulated proton charge. Arrange- 
ments are such that the proton Integrator operates only while the 
Geiger counter-scalor is accepting and registering gamma-ray counts. 

The gamma-ray counts must be corrected for "background". This 
background is a function of the generator-accelerator voltage, since 
x-rays from the generator provide a sizable proportion of said back- 
ground. Experiments made v/ithout a target Indicate that accumulated 
background can be estimated by dividing the time (seconds) for the 
prescribed number of total counts in a single run at a current set- 
ting by a factor, given in Table VI. These factors are valid only 
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Table VI 



Factors for Detemlnation of Background 



V/hen Operating Van de Graaff Generator 
Type Tube; Technical Associates Beta Counting Tube TA-Bl. 



Location; As close to the target as possible - separated from 

same by thickness of target holder, 1/8" of lead, and 
about three millimeters of air. 

Position; (One tube) - Side to target, with center line at mid- 
target height, (Ti' 70 tubes in parallel) - One above the 
other, sides to target, v;ith mean of positions of 
centerlines at mid-target height, , 



Shielding; Tube(s) mounted vjithin a two-inch thick lead box, with 
hole in same Just large enough to permit ready inser- 
tion of target holder. 



Anal|tzing Magnet 
Current (amp.) 


Factor by bhlch to Divide Length of 
Run in Seconds to Get Background Counts 




One Tube Two Tubes in Parallel 


.55 


3.0 2.0 


.60 


2.8 2.0 


,65 


2.6 2.0 


.70 


2.4 2.0 


.75 


2.2 1.9 


.80 


2.0 1.7 


.85 


1.8 1.5 


.90 


1.5 1.4 


,95 


1.2 1.2 
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If the geometry, tube type, and other factors are reproducible; end 
these variables are indicated in the table. 

The proton integrator count for a single run is obtained by sub- 
tracting the register values before and after the run. The values thus 
obtained also require correcting, as a certain amount of charge leakage 
occurs through paths other than that which actuates the register. By 
experimentally checking the integrator register counts per unit time 
against galvanometer readings of current frcra the Faraday cage to the 
integrator, it has been found that integrator readings must be increased 
by an amount corresponding to about 35 counts per minute of time for 
the run. Such counts, when corrected, thus provide a more constant 
ratio between the figure for total integrator counts and the actual 
number of protons striking the target. 

The ’’yield", on an arbitrary scale, is the result of dividing 
the corrected gamma-ray counts by the corrected proton integrator 
count. 

These data are token at successively higher magnet cuinrent 
settings. The desirable current Increments between successive 
readings for our experimental vjork vjas as follows: for sharp 

peaks, current increments of 0,5 na. are necessary; for brooder 
peeks, data every 1 ma. are sufficient. 

If gas stopping is being studied during any seouence of indi- 
vidual runs, the manometer, thermometer (see Figure 4), and time of 
day are recorded just before gas introduction, just after gas intro- 
duction, just before the gas is finally removed, and just after the 
gas Is removed. Readings at Intervals betiveen start and finish may 



also be made if desired 



The yield can be plotted against magnet current, or If the 
analyzing current-proton energy relationship has been established 
(called "calibration" herein) the yield can be plotted versus nroton 
energy to give a more pronsr si>ectruiii. For our purposes only the 
first type of plot has been required, as it is ixjssible to read the 
peak maxlraum values and the statistical measurements of v;idth 
directly in terms of magnet current and then convert t"^ the proton 
energy scale. 

Table VII gives a sample set of data for a single peak, v.’lth 
yield computations. Figure 8(a) shavs the plot of the resonance peak 
(unnhifted) for this particular set of data; Figure 8(b) is a plot of 
the same peak, shifted by insertion of stopping material. The 
probable error of the position of each noint is indicated by the 
size of the crosses. The error in horizontal direction is based on 
personal observation and judgment; the error in the vertical direction 
is detemlned frora st-ndard statistical deviations of the total num- 
ber of gamm-ray counts (P,E, = 0,67 yuo, of Counts ). In the latter 
case, the actual probable error is somewhat greater, but the largest 
errors ca’ie frcxn counting statistical errors and give p rough basis 
for estimating the extent to v;hich our smooth curve has to folloi*? 
the experimental points . 

E, Target Rreporation. 

The target ronterials \ised were either lithivun fluoride or 
aluminum. Ttie targets v/ere formed by coating a blank tantalum 
disc with a thin coating of the target substance. ‘;ii 0 coating vms 



Table VII 



Sample Set of Data for Computation of Yield and Plotting of Peak 
Peak Covered: 540.4 kev - . Target Used: LiF-2 

No Stopping Nedium Interposed. 

One Geiger Tube Used, 



Magnet 

Current 

(amp) 


Proton 

Int. 

Counts 


Time 

(sec) 


C orrec . 
to Prot. 
Counts 


Net 

Proton 

Counts 


Gamma 
C ounts 


B'grnd. 

Corr.to 

Gam.Cts 


Net 

Gamma 

Cts. 


Yield 


.4900 


49 


28.6 


17 


66 


32 


9 


23 


.3 


10 


82 


45.0 


26 


108 


64 


14 


50 


.5 


20 


83 


56.6 


33 


116 


64 


18 


46 


.4 


50 


59 


34.0 


20 


79 


64 


11 


53 


.7 


40 


94 


57.2 


33 


127 


128 


18 


no 


.9 


50 


84 


46.8 


27 


in 


128 


15 


113 


1.0 


60 


80 


43.4 


25 


105 


256 


14 


242 


2.3 


70 


69 


41.8 


24 


93 


256 


13 


243 


2.6 


80 


61 


37.7 


22 


83 


512 


12 


500 


6.0 


90 


62 


44.0 


26 


88 


1024 


14 


1010 


11.5 


.5000 


43 


32.0 


19 


62 


1024 


10 


1014 


16.4 


10 


71 


38.4 


22 


93 


2048 


12 


2036 


21.9 


20 


64 


37.4 


22 


86 


2048 


12 


2036 


23.7 


30 


105 


60.0 


35 


140 


4096 


19 


4077 


29.1 


40 


136 


72.0 


42 


178 


4096 


23 


4073 


22.9 


50 


57 


33.6 


20 


77 


2048 


11 


2037 


26.5 


60 


53 


36.3 


21 


74 


2048 


11 


2037 


27.5 


70 


83 


53.6 


31 


114 


2048 


17 


2031 


17.8 


60 


61 


38.0 


22 


83 


1024 


12 


1012 


12.2 


90 


70 


46.2 


27 


97 


1024 


14 


1010 


10,4 


.5100 


83 


62.6 


37 


120 


1024 


20 


1004 


8.4 


10 


224 


84.0 


49 


273 


512 


26 


486 . 


1.8 


20 


85 


55.0 


32 


117 


256 


17 


239 


2.0 


30 


67 


39,0 


23 


90 


128 


12 


116 


1.3 


40 


77 


43.6 


25 


102 


128 


14 


114 


1.1 


60 


70 


41.5 


24 


94 


128 


13 


115 


1.2 


80 


110 


62.7 


37 


147 


128 


20 


108 


.7 



These results are plotted in Figure 8(a). 
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Table VIII 



List of Resonance Energy Values for Proton-Ganra Process in 

I'aterials Used as Targets 




(2) Brostrom, Huus, and Tangen, 1947, 



performed by n rather standard technique of evanoration in a 
vacuxim {Strong, 1945) . The materials selected vjere those for 
which the energy resonance levels and peak widths are known with 
reasonable aceuracy. The levels for lithium, fluorine, and 
aluminxm os obtained from the latest and best references are 
tabulated in Table VIII. 

A succession of targets were made end used, before some were 
obtained which were entirely satisfactory. satisfactory target 
must have the folloi^ing qualities: 

(a) Tlie target should be thick enough to give high yields of 
gamma-rays per proton count, thus providing good counting "statistics” 

(b) On the other hand, the target thickness should be small enough 

that target thickness, T, and (section II. B) can be minimized, 

thus minimizing the effect of errors in their estimation, jftnd under 
any circumstances T must be much less than the energy difference 
between adjacent peaks, to ovoid interference froKi other peaks vdien 
a certain one is analyzed. 

(c) A fair uniformity of taiget thickness should be provided, so 

that when a proton beam is somev/hot diffused over the target, the 
same values of T and will apply to all portions of the beam. 

F. /g ia lysis of A uxiliar y Statistical Variables. 

Preliminary to the primary calculations of stopping power are 
those calcoilations of the other Vfiriables v/hlch are a part of eoua- 
tions (II-B-4) to (II-B-V), inclusive. These calculations require 
a knowledge of the calibration relation betv/een analyzing magnet 



CTirrent and proton energy. Unfortunately an accurate determination 



of this relation requires in turn a knowledge of target thickness, T, 
which is one of the variables referred to above. Originally a rough 
assumption as to the nroper calibration relation had to be assumed 
and calcxilations made by the process of successive aoproxlmations. 

In this report v;e may start I'd.th reasonably accurate information on 
the mean calibration ciirve, which has actually been obtained at a 
later stage in the experiment (See Section III.H). 

Within a limited range it is found by plotting calibration data 
on log-log grajii paper that the calibration relation can be expressed 
by a formula of the type: 



where E represents proton beam energy (Hev) and I Is anal:,TSing 
magnet current (aiaperes), (See figure 9.) 

Besides K and a , we wish to find the value of F, ^rtiere 



which is used to translate magnet current variations into energy 
variations at some particular peak. It has been found, using data 
to be presented later (Section III,!!), that the values for K, a , and 
F given in Table IX are fairly accurate. 

Let \is nav analyze sane of the statistical distributions In 
energy a little more closely. It is necessary for our purposes that 
the statistical siiread in these distributions be expressed In terms 
of the standard deviation about the mean or its square, the variance. 
Many of the distributions are of such a type that the spread is not 
usually expressed thus, and calculations must be made to permit 



E = K I® , 



(III-F-1) 



F 




(lII-F-2) 



Table EC 



Values of constants K, a, and F { s dE/dl) in the 
Magnet Current-Proton Energy Relationship: E = K I® . 



Current 

(amp) 


Approx. 

Energy 

(Kev) 


K 


a 


F 


,50 


.344 


1,355 


1.975 


1.362 


,55 


.417 


1,353 


1.972 


1.490 




.494 


1.350 


1,970 


1.620 


,65 


,579 


1.348 


1.967 


1.748 


.70 


.669 


1.345 


1.960 


1.870 


.75 


.765 


1.339 


1.950 


1.986 


.80 


.868 


1,332 


1.935 


2.09 


,85 


.972 


1.323 


1.920 


2,19 


,90 


1.076 


1.315 


1.905 


2,28 


.95 


1,181 


1.305 


1.887 


2.35 



translation to those statistical variables we employ in common for 
all the distributions. Resonance curves for example have their 
spread expressed in terms of a width of peak at half the maxlmvim 
height, 

1, Method of Obtaining -A , 

As Just noted, it is necessary to translate the natural peak's 
v/ldth at half height, designated as Hiat * the standard deviation 
n V/e are immediately faced v/lth a somenftiat difficult situation, 

since the theoretical computation of jjg-t takes one into compli- 
cations of both a mathematical and a practical nature. 

The shape of the theoretical resonance curve is given by the 
well known Breit-VJigner formula (Breit and V;igner, 1936) , modified 
for the proton-gamma reaction (Bethe, 1937): 




cr - 



» 



(III-F-3) 



• . Ji_ - ^ 

^ /T (E - E + Kg 

where O" is the probability cross-section for an atom to undergo 
the proton-garnina reaction. is the barrier penetration factor 

for the proton in trying to pierce the Co\ila.ib potential barrier of 
the atom, and at high proton energies (compared with Coiilomb barrier 
height) it is dlose to unity, v/hereas at low proton energies it 
equals approximately exp ) • 

This expression is difficult to deal *ith mathematically, but 
we may make further simnlif Ications for the sake of discussion. Since 
most of the resonances used are extremely narrow conpared to the proton 
energy (E) exciting them, the factor Gp/i^ is very sloi^ly varying in 
the immediate vicinity of the peak, and can be taken as constant. V^e 
then have a resonance peak with a shape given by the formula: 



cr r 



K 

(E - + Kj, 



(III-F-4) 



By 



definition, If , is then given by 



fi' 



nat 



K dE 



1 Kp/{E - E )' 
« res 



(III-F-5) 




K dE 



The difficulty arises v/hen we try to decide v/hat limits to take for 
E^ and E_ , on either side of the peak. E_ cannot go below sero, 

but there is no natural upper limit to select. The value of the 

- - 



denominator is finite v;lth limits zero to infinity; but the value 
of the numerator increases Mthout limit as an’^roaches infinity, 
Uius the value of the standard deviation will depend unon ;vhere we 
choose to take this upper limit and becomes somewhat arbitrary. The 
inclusion of the factor G^/ /W does not simplify our task, since 
the numerator still diverges \vlth Increasing value of . 

- -v ’ 

To avoid the decision as to where precisely v;e cut off the 
"tails" of the distribution, we shall employ the same apnroximntion 
used by I.!adsen and Venkatesv/arlu (1948b), Noting that the resonance 
curves are rather similar in shape to a Gaussian, or "normal", dis- 
tribution curve, one may assume that the ratio of the standard 
deviation to the v/idth at half -maximum height is the same for a 
theoretical resonance curve as for the true Gaussian shaped curve. 
Fron tobies of the latter function it is easily found that as a 
consequence 

•C nat * Ct • (III-F-6) 

Table VIII, \7hich gives the resonance v;idths for the target 
materials used, also gives the values of ■^nat computed by the above 
formula, 

2, Method of Determining . 

The distribution curve for energy loss in the target material 
is quite simple, I'/e may assume that the ’>roton stopping ponder is 
practically constant over the small energy variation the beam under- 
goes in the thin target, so that this distribution curve has the 
shape of a rectangle. That is, if y(AE^) is the irequency veiriable, 
then it is defined as: 



- - 



y = - 


0, for AE^ < 0 and AE,p > T ; ^ 


(III-F-7) 




k, for 0 < < T . J 




Then, 


/T 




II 


: ^ \ k (^E - d(AE) , 






T^A2 , 


(III-F-8) 


so thnt 




= ,288 T 


(III-F-9) 



3* Determination of ilg , 

The distribution of energies of protons passing through the 
slit cannot be knovm with the degree of certainty one would like. 
Experience indicates that for a narrow slit the beam energy stabili- 
zer can tisually hold the center of the beam within the edges of the 
slit quite well. The maximum variation in proton energy which will 
still permit the beam center to pass throvigh the slit has been ccin- 
puted by Smith (1958) to be: 

• (m-F-io) 

;^ere w Is the width of the slit in millimeters, and and E 

are in the same units. This formula is probably not highly accurate 
in view of the assumptions and approximations involved in the derivation. 

To provide a check, the shift in location of a peak at 669 kev was 
determined experimentally when the slit position was shifted by 
millimeter. The shift in peak position was 1.2 ma, of magnet current, 
which corresponds according to Table IX to a bean energy shift of 
lo87 X 1.2 m 2.25 kev. These values are consistent \d.th an equation: 






(III-F-U) 



^ • 

This e:q)erinent probably does not give a highly accurate 
answer, but it is considered an improvement numerically over the 
theoretical formula’s prediction; and ve shall use equation 
{ III-F-U) therefore. 

If \ie assume that energy variations in the accelerator are 
such that a proton is equally likely to pass through the slit inth 
any value in the corresponding spread of , we have a rectangu- 

O 

lar distribution; and as in equation (lII-F-9) we obtain: 

i2_ r .288 Ew/150 - Ew/ 520 . (III-F-12) 

Since for all the experiments hei’ein described, w equals one 
millimeter, the formula becomes: 

= E/520 , (lII-F-13) 

viiere the mits of E and Ag are the same. 

4* Determination of 

The factor is impossible to obtain directly, and it can 

only be guessed at ’..'hen the resonance peaks' ^adths are analyzed « 
There is no reason to believe that it is even constant from one 
day to another under \jhat are presumed to be similar operating 
conditions. Since current variations cause greater beam voltage 
variations v/ith higher energies, \Je may guess tliat will general- 
ly increase \ath increasing proton energy, though tliis is by no 
means certain. 






G» Analysis of Unshifted Resonance Peaks and Detemination of 



Tarf^et Thictoesses, 

Table X gives the analysis of the statistical variations of 
nost of the resonance pcaJ:s studied ^dthout stopping nedia interposed. 
This analysis is based upon fonaula (II-D-5), and shows all computed 
values for . Observation of those results leads to the 

folloimig conclusions: 

(a) Thci-e are conparatively large uncertainties present in all of 

the data presented, so tliat the final probable error in the estima- 
tion of is on the order of 1 kev or mere. This is especially 

evident from the size and number of negative values obtained, since 
all values should Ite positive for variances, 

(b) Thei'e scens to be little reason to believe tliat is lai'ge 
enough to consider, except in a few cases, V7e shall, therefore ignore 
it in our computations bat mast consider its possible influence vrficn 
the matter of errors in results is discussed. 

(c) For Icn; energy values is at a maxiranm, \;hereas is at 

a minimum. Therefore in determination of target thickness, greater 
weij^t shotild be given to data, obtained from lov/or energy peaks. 

In ofdcr to ej:preES the target thicloicss in terns of stopping 

power, ire must standardize e’ch values to a specific value of proton 

energy (1 I'ev is ueial). This means that after having determined 

target thickness at a certain proton energy we must convert this 

value to thiclcness for a 1 Mev y^roton beam. Tliis can be done if we 

know the relative stopping poircr of the target substance at 1 Kev 

compared to the stopping power of the target at the energy of the 

~ I 'J ~ 



Line 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 



Table X. {Tvjo Pages) 

Aaalysls of Contributions to Variances in Unshifted 

Resonance Peaks, 



Date 


Target 


Peak 

Energy, 

(kev) 


I'agnet 
Current , 
(anp) 


F 

(Table DC) 


8-29-52 


LiF-1 


340.4 


.50 


1.562 


ft 


ff 


441.1 


.566 


1.533 


ft 


If 


669 


.70 


1.870 


9-24-52 


Al-3 


503 


.60 


1.620 


It 


. If 


630 


.68 


1.820 


If 


II 


986 


.85 


2.170 


fl 


n 


1112 


.91 


2.255 


9-26-52 


Al-3 


630 


.68 


1.820 


10-16-52 


LiF-1 


669 


.70 


1.870 


ft 


It 


873.5 


.80 


2.086 


n 

> 


II 


935.3 


.83 


2.135 


10-17-52 


Al-3 


630 


.68 


1.820 


n 


n 


986 


.86 


2.185 


ft 


If 


1112 


.91 


2.255 


10-23-52 


LiF-1 


441.1 


.566 


1.533 


12-2-52 


LlF-2 


340.4 


.50 


1.362 


1! 


It 


669 


.70 


1.870 


ft 


II 


873.5 


.81 


2.103 


ft 


ff 


935.3 


.84 


2.153 


12-18-52 


LiP-4 


441.1 


.57 


1.540 


ft 


II 


669 


.70 


1.870 


If 


If 


873.5 


.805 


2.095 


If 


ff 


935.3 


.835 


2.144 
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Table X. (Concl.) 



Line 


(ma) 


■ 

(kev) 


■Oo 

(kev^) 


■‘^nat 
(kev^) 
Table VIII 


*-s 

(kev^) 

(III-F-13) 


SlF + 

-*^T M 

(kev^) 
{II-E-5) 


1 


1.90 


2.60 


6.75 


0.94 


0.43 


5.38 


2 


4.95 


7.60 


57.7 


26.0 


0.72 


31.0 


3 


3.34 


6.25 


39.1 


9.9 


1.66 


27.5 


4 


0.80 


1.30 


1.69 


Negl. 


0.94 


0.75 


5 


0.90 


1.64 


2.69 


V 


1.47 


1.22 


6 


0.84 


1.82 


3.31 


0.01 


3.60 


- 0.30 


7 


1.35 


3.05 


9.30 


0.01 


4.58 


4.71 


8 


1.1 


2.00 


4.00 


Negl. 


1.47 


2.53 


9 


1.65 


3.10 


9.61 


9.9 


1.66 


- 1.95 


10 


1.2 


2.50 


6.25 


4.84 


2.82 


- 1.41 


11 


1.9 


4.05 


16.40 


11.23 

% 


3.24 


1.93 


12 


1.15 


2.10 


4.41 


Negl. 


1.47 


2.94 


13 


0.93 


2.03 


4.12 


0.01 


3.60 


0.51 


14 


1.03 


2.52 


5.38 


0.01 


4.58 


0.79 


15 


3.45 


5.30 


28.1 


26.0 


0.72 


1.40 


16 


3.58 


4.88 


23.8 


0.94 


0.43 


22.4 


17 


2.82 


5.27 


27.8 


9.9 


1.66 


16.2 


18 


2.97 


6.25 


39.0 


4.84 


2.82 


31.3 


19 


3.02 


6.50 


42.2 


11.23 


3.24 


27.7 


20 


3.2 


4,93 


24.3 


26.0 


0.72 


- 0.24 


21 


1.43 


2.67 


7.13 


9.9 


1.66 


- 4.43 


22 


1.2 


2.51 


6.50 


4.84 


2.82 


- 1.36 


23 


1.85 


3.97 


15.76 


11.23 


3.24 


1.29 
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Table XI 



Conputations of Target Thickness (kev) 



Target 


Peak 

Energy 

(Kev) 


Conv. 
Factor 
s YT 


T 

(ker^) 


ill 

Ocev) 


T ■ 

-^.p/.288 


T at 
1 Mev 
(kev) 


LlF-1 


,5404 


.583 


5.38 


2.32 


8.05 


4.7 


TT 


.4411 


.664 


31.0 


5,57 


19.3 


12.8 


Tt 


,44U 


.664 


1.4 


1.18 


4.1 


2.7 


Tf 


.669 


.818 


27.5 


5.24 


18.2 


14.9 


fl 


,669 


.818 


nega. 


- 


- 


- 


n 


.8735 


,935 


n 


- 


- 


- 


n 


.9353 


.970 


1.93 


1.39 


4.82 


4.7 


LlF-2 


.3404 


.583 


22.4 


4.73 


16.4 


9.56 ■ 


fl 


.669 


.818 


16.2 


4.02 


13.9 


11.4 


If 


,8735 


.935 


31.3 


5.60 


19.4 


18.1 


n 

♦. 


.9353 


.967 


27.7 


5.26 


18.2 


17.6 


LiF-4 


.4411 


.664 


nega. 








n 


.669 


.818 


Tt 


- 


— 


- 


ft 


.8735 


,935 


n 


- 


- 


- 


n 


,9353 


.967 


1,27 


1.135 


3.94 


3.8 


Al-3 


.503 


.710 


0.75 


0.866 


3,00 


2.13 


ff 


,630 


.794 


1.22 


1.105 


3,83 


3.04 


ft 


.630 


.794 


2.53 


1.59 


5.51 


4.37 


fl 


.630 


.794 


2.94 


1.714 


5.95 


4.72 


n 


.986 


.993 


nega. 


- 


- 


- 


n 


.986 


.993 


0,51 


0.714 


2.47 


2.45 


ff 


1.112 


1.055 


4.71 


2.17 


7.53 


7.94 


ff 


1.112 


1.055 


0,79 


0,889 


3.08 


3,25 



particular pea]c analyzed. 

Bohr ( 1948 ) indicates that vdtliin this energy region which v;e 
are concerned with, the stopping power is inversely proportional to 
the square root of proton energy. To the extent that data is 
available in tills region, such dependence is roughly confinacd 
(ISadsen and Venicateswarlu, 1948aj Kahn, 1953j and Figure 11 herein). 
In viei7 of the approxinate nature of the data given in Table X, it 
is sufficiently accurate to convert target thickness to the 1 I'ev 
standard energy ty nultiplying the target thicloiess by the square 
root of the proton energy at which the specific laeasurcnent was made. 

Table XI gives the target thickness analysis, with the aid of 
the data derived in Table X. In view of the approxinate nature of 
the final results, it is desirable to use other means to aid in the 
selection of definite answers. If the targets are thin the naxinun 
value of the yield at a certain resonance pealc vn.ll be a rou^ 
measure of the relative tliickness of the target. Thus a comparison 
of the maxima of the various peaks for the various litliium fluoride 
targets will add additional information on viiich to base final judg- 
ment. It is expecially desirable to do this for determination of the 
thickness of target LLF— 4» since it is evidently so thin as to make 
detection from an experimental pealc width practically Impossible, 
Table XII gives this data, as talcen from ejperimental results. The 
table indicates that in a rou{^ sort of v;ay, the relative yield 
froa the various lithium fluoride targets are in the ratio: 

LiF-1 : UF-2 : UF-4 = 0.3 : 1 : 0.1 



-y/ 
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Table XII 



Resoxiance Peak Ikudna for Lithiun Fluoride Targets 




* Evidently the target LiF— 1 has a sanex/hat variable 



thickness. We can speak only of its average value, 
n 

^ Counting geometry someiiiat variable. 



Table XIII. 

Target Thicknesses (kev) for Various Proton Energies 



Litliium Fluoride Targets 


Aluminum Target 


Proton 








Proton 




Energy 


LLF-1 


UF-2 


LiF-4 


Energy 


Al-3 


(Mev) 








(llev) 




1.000 


4.0 


12.0 


1.0 


1,000 


3.0 


• 3404' 


6.8 


20.6 


1.8 


.503 


4.2 


.4411 


6.0 


18.0 


1.6 


.630 


3.8 


.598 


5.2 


15.6 


1.2 


.771 


3.4 




4.8 


14.6 


1.2 


.986 


3,0 


.8735 


4.2 


12.8 


1.0 


1.112 


2.8 


.9353 


4.2 


12.4 


1.0 







On the basis of all information available, ue can nov; nalco an 
estimate of target thicknesses for all targets. Table HII gives 



i 



the thickness estinated at 1 Mev, and also gives the coii^juted values 
of target thicloioss at various resonant energies. In the cocqjutation 
the target thidoiess is assumed, ab before, to be inversely propor- 
tional to tlio square root of the energy. 

H. Enorfr;~-Currcnt Calibration. 

Since vre can no\/ find the value of for each peak measured, 
ly use of equation (lI-B-4)» we may find the calibration curves for 
accelerator beam energy (E^) versus analyzing magnet current (I). 

All data from eo^jeriment and previous analysis is listed in Table 
HV in such a fashion that the calibration curves can be put on a 
logarithmic scale. The oirves themselves are depicted in Figure 9. 

Table XV is a list of data similar to Table XIV, except that 
the data are for single resonance peaks and are not taJeen as part of 
a coEplete calibration scries. Their use is indicated below. 

In using and interpreting the calibration cmrves, several 
significant points have been discovered: 

(a) If the data is plotted to a logaritlmic scale, vertically and 
horizontally, the calibration line dra'wn is almost a strai^it line, 
^'his indicates that an crspirical formula of the type E r K I^ 
is approximately valid over the \iiole energy range, and is quite 
accurate ^jithin smaller limits, (This justifies equation (III-F-I).) 
This is a change from the previa >s custom (Grove, 1950; Taylor, 1952; 
Russell, 1952), by \jhich data :-jas plotted as yr versus I and the 
best strai^t line dravm. This older typo of plot presumes that an 
empirical formula of the type E = K (l+b) is valid. Plotting 
of a typical set of calibration data from Table XIV according to the 



Table XIV» (Ty:o Pages) 



Calibration Data 



Date; 
Target ; 
Slit 
Setting, 
mm. 


^nat 

(Hev) 


(I.:ev) 


Eo 

(Dev) 


log Eq 


I 

(amp) 


log I 


8-29-52 


.3404 


.0034 


.3438 


T.5363 


.4973 


T.6966 


LiF-1 


.4411 


.0030 


.4441 


1.6475 


.5665 


1.7532 


4.25- 


.669 


.0024 


.6714 


1.8270 


.6996 


1.8448 




.8735 


.0021 


.8756 


1.9423 


.8020 


1.9042 


9-26-52 


.503 


.0021 


.5051 


1.7034 


.6032 


1.7805 


A 1-3 


.630 


.0019 


.6319 


1.8006 


.6763 


1.8301 


4.75 


.986 


.0015 


.9875 


1.9945 


.8563 


1.9326 




1.112 


.0014 


1.1134 


0.0467 


.9095 


1.9588 


10-16-52 


.669 


.0024 


.6714 


1.8270 


.7023 


1.8465 


LiF-1 


.8735 


.0021 


.8756 


1.9423 


.8030 


1.9047 


4.75 


.9353 


.0021 


.9374 


1.9719 


.8315 


1.9199 


10-17-52 


.630 


.0019 


.6319 


1.8006 


.6789 


1.8318 


Al-3 


.986 


.0015 


.9875 


1.9945 


.8573 


1.9331 


4.75 


1.112 


.0014 


1.1134 


0.0467 


.9103 


1.9592 


12-2-52 


.3404 


.0103 


.3507 


1.5449 


.5036 


1.7021 


LiF-2 


.669 


.0073 


.6763 


1.8301 


.7036 


1.8473 


3.75 


.8735 


.0064 


.8799 


1.9444 


.8074 


1.9071 




.9353 


.0062 


.9415 


1.9738 


.aS77 


1.9231 


12-3-52 


.503 


.0021 


.5051 


1.7034 


.6115 


1.7864 


Al-3 


.630 


.0019 


.6319 


1.8006 


.6838 


1.8349 


3.75 


.986 


.0015 


.9875 


1.9945 


.8640 


1.9365 



Slit setting ecuals reading of nopition ricroineter minus one- 
balf the reading of the slit v/idtl ricrometer. It essentially 
indicates the nosition of the middle of the slit. 
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Table XIV. (Concl. ) 



Date; 
Target ; 

Slit 

Setting, 

ram. 


^nat 

(I'ev) 


It 

(liev) 


Eo 

(Fev) 


log Eq 


I 

(amp) 


log I 


lD-18-52 


.4411 


.0008 


.4419 


1,6453. 


.5690 


1.7551 


LiF-4 


,669 


.0006 


.6696 


1.8258 


.6995 


1.8448 


3.75 


.8735 


.0005 


.8740 


1.9415 


.8052 


1.9059 




.9353 


.0005 


.9358 


1.9712 


.8334 


1.9209 


12-19-52 


.503 


.0021 


.5051 


1.7034 


.6095 


1.7850 


^-3 


.630 


.0019 


.6519 


1.8006 


.6817 


1.8336 


4.25 


.771 


.0017 


.7727 


1.8880 


.7562 


1.8786 




.986 


.0015 


,9875 


1.9945 


.8620 


1.9355 




1.112 


.0014 


1.1134 


0.0467 


.9176 


1.9627 


1-9-53 


.4411 


.0008 


.4419 


I. 6453 


.5700 


1.7559 


LiF-4 


.669 


.0006 


.6696 


1.8258 


.7028 


1,3468 


3.50 


.8735 


.0005 


.8740 


1.9415 


.8074 


1.9071 




.9353 


.0005 


.9358 


1.9712 


.8367 


1.9226 


1-19-53 


.669 


,0006 


.6696 


1.8258 


.7014 


1.8460 


Li?-4 


.8735 


.0005 


.8740 


1.9415 


.8080 


1.9074 


3.50 


.9353 


.0005 


.9358 


1.9712 


.8367 


1.9226 


1-22-53 (an) 


.659 


.0006 


.6696 


1.8258 


,7040 


1.8476 


LiF-4 


.8735 


.0005 


.8740 


1.9415 


.8072 


1.9070 


3.50 














1-22-53 (pm) 


.669 


.0006 


.6696 


1.8258 


.7014 


1.8460 


LiF-4 


.8735 


.0005 


.8740 


1.9415 


.8056 


1.9061 


3,50 


.9353 


.0005 


.9358 


1.9712 


.8335 


1.9209 
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Table XV» (Tt’fo Pages) 



Data for Single Peaks (Unshifted) 



Date; 
Target ; 

Slit 

Setting, 

mm. 


E 4 . 

nat 

(Mev) 


(rev) 


Eo 

(rev) 


log 


I 

(amp) 


log I 


10-23-52 

LlF-1 

4.75 


.4411 


.0030 


.4441 


1.6475 


.5670 


1.7536 


10-28-52 

LU-1 

4.75 


.9353 


.0021 


.9374 


1.9719 


.8320 


1.9201 


10-31-52 

LiF-1 

4.75 


.9353 


.0021 


.9374 


1.9719 


.8303 


1.9192 


1-12-53 (am) 

LiF-4 

3.50 


.669 


.0006 


.6696 


1.8258 


.7009 


1.8457 


1-12-53 (pm) 

LlF-4 

3.50 


.669 


.0006 


.6696 


1.8258 


.7001 


1.8452 


1-13-53 

LiF-4 

3.50 


.669 


.0006 


. 6696 


1.8258 


.6998 


1.8450 


1-14-53 

LiF-4 

3.50 


.669 


.0006 


.6696 


1.8258 


.7002 


1.8452 


1-15-53 

LiF-4 

3.50 


.669 


.0006 


.6696 


1.8258 


.7023 


1.8465 


1-19-53 

LiF-4 

3.50 


.669 


.0006 


.6696 


1.8258 


.7024 


1.8466 
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Table XV, (Conel«) 



Date; 
Target ; 
Slit 
Setting, 
mm. 


E * 

nat 

(I'ev) 


(I.iev) 


(I.!ev) 


log 


I 

(amp) 


log I 


1-20-53 

LiF-4 

3.50 


.669 


.0006 


.6696 


1.8258 


.6996 


1.8449 


1-21-53 

LiF-4 

3.50 


.669 


.0006 


.6696 


1.8258 


.7030 


1.8470 


1-23-53 

LiF-4 

3.50 


.669 


.0006 


.6696 


1.8258 


.7018 


1.8462 
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previous custoB and ccoaparing such a plot vdth the double-loj^ith- 
nic plot has indicated that the latter plot has less curvature than 
tlie coiTe^onding plot the older nethod* In viev; of the fact, 
hovever, tliat a slif^t curvatia*e does exist for the calibration 
curves, even on the logarithnic plots, it vas considered more 
accurate to use the curves directly rather than to use a straight 
line fomula giving the closest fit over the uhole energy range, 

I 

The cc^jirical fonaula (lII-F-l) and Table IX are to be used only for 
coETOutations at localized ejiergy regions, 

(b) Different calibration series give curves which are generally 
sli£^tly displaced frcB one another, although they are essentially 
parallel. This indicates that personal errors during a series of 
runs are snail compared to those discrepancies betvjeen 
different daily series. Every effort \jas made in doing the experi- 
mental vork to use a standardised technique in operation of the 
accelerator and auxiliary apparatus; but scane variation from day 
to day seemed imavoidable. Such variations may be due to: (l) lade 

of precise control as to the horizontal and vertical position of the 
beam aids (as indicated ty the neutral beam, imdoviated ty the 
magnetic field); (2) lack of control over the tenperatwe of the 
magnet itself, Tlic first effect means that inhomogeneities in the 
magnetic field nay be sufficient to be noticeable. The latter effect 
is related to the fact that magnetic elements in an electromagnet 
are more easily aligned when the iron is at a higher ten^orature, 
so tliat the same current gives a sli^^tly greater deviating flux. 

The fact that different curves arc closely parallel neans that once 



a fo\t typical calibration curves have been plotted to give the 

general trend, farther curves nay bo indicated by having data for 

only a single calibration point on the cui’ve. This fact has been 

utilised in the latter part of the research, and Table XV gives 

the data for these single calibration points* 

(c) It is of some interest to note the actual variation of the 

calibration data for one specific value of energy, to give an idea 

as to the spread in the calibi'ation curves over a sis nonths* period, 

A chedc of Tables XIV and XV shows that for the 669 kev fluorine 

resonance peak, by use of the sane target and the saiae slit position, 

the average value of the analysing current was ,7014 anperes, with 

extremes of ,6696 and ,704D arroeres. The standard deviation of all 

readings fron tlie moan is ,0014 amperes. If on the other hand, one 

were given the current at about this point, the estrerie energy variation 

corrc£ponding to this calibration carve variation would be 6,2 kev 

(determined by using factoi'* F, Table IX), The standard deviation 

in energy cones out to be 2,6 kev from the mean calibration line, or 

a probable error, if one used the laean line, of about 1,7 kev. This 

fact indicates that if vc have some additional infoz’nal.ion as to the 

proper position of the true calibration curve, such as by a single 

calibration point, a much better choice can be made for energy 

correeponding to a certain value of magnet current; and the probable 

error in such a case would undoubtedly bo substaJitially less than 1 

kev. It shooild be noted, that this conclusion is justified only if 

the technique is used for data talcen on the sane roagnet cycle as the 

calibration point. A recycling of the magnet may require a new curve, 
(For cxarroles, sec Table XIV, data for 1-22-53, a,m, and p,n, ) 
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IV» Experipental Results on Stopping Poorer 



A, Possible Sovirces of Systematic Error and Corrections Regaired , 

When the experimental data on stopping are analyzed, the possi- 
bility of certai-n consistent factors requiring correction mst be 
considered, 

1, Nuclear Ck)l]ision Losses , 

It can be visualized that proton-nuclear col3.isions, thou^ 
elastic, vriU result in a certain small amotant of energy loss by 
the proton in the ■laboratory system. The usual formulas for 
stopping po\;cr do not include this effect, and therefore in theory 
the nuclear stopping effect should be subtracted from experimental 
data in order to j^ve a result vahidi can bo compared vd.th theory. 

Bohr ( 194 s) has derived a formula for nuclear stopping, givaa 
as follo^;s: 



In order to determine the magnitude of this effect, an example 
v;as calculated for laypothetical circumstances conducive to a large 
ansxier. For energr loss due to nuclear collisions of protons of 



thro'ugh a 0.05 mil nickel foil, the calculated valiie of E is 

V 

0.171 kev. It can be seen hy comparison \ri.th electronic stopping 
energy losses for the foils (Figure 10) that the nuclear effect is 




y 



(iV-Ar-l) 




(IV-A-2) 
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only on the order of 0.1^ of the total stopping effect. The nuclear 
stopping my therefore be neglected. 

2. Increase of Path Lcaifrbh ty ScatterinfT . 

Another effect, also due prinarily to nuclear collisions, is the 
deviation from the beam axis of many protons, especially for the case 
of slov; protons through large thidcnesses of matter. This effect is 
observable visually during stopping e^qjeriments with gases. The ex- 
citation of the gas atoms tgr the proton beam brings about the emission 
of light xiiich marks the proton paths. It can be seen that the beam 
does indeed diverge in a very naiToi; cone from the entrance foil to 
the gas chamber. Such effect increases the average amotint of mtter 
traversed ty the beam. The geometry of our diamber (Figure 4 ) is 
such that the target subtends at the entrance foil an mgle of 
5.6° to each side of the axis. Any proton deviated so much as to 
pass outside this ccne does not hit the target and cannot contribute 
to the gamma ray yield and thus has no effect on e:q5crimental values 
of stopping po\;er. Only at very low proton energies (as noted visu- 
ally) does there exist an appreciable deviation of the protons avray 
from the axis. . Even in such circumstances the correct Lon for longer 
path length is only a fraction of cos 5.6° , or less than one-lialf 
of one percent. At hi^er energies the correction is still less. 

Since this correction is much less than the over-all experimental 
error expected, i/e shall iguore it. 

3 . Carbonization Foils . 

There is a tejndency for all parts exposed to the beam to accu- 
mulate a dark deposit of liiai. is called "carbon”, though it 



probably includes a nodcrate percenta[;e of hydrogen. This deposit 
is a well-knoxjn problen for nany accelerators. It is presumably a 
result of deposition of volatile hydrocarbons on surfaces plus sub- 
sequent molecular brealcdovai by proton bombardment and by localized 
heating. It is even possible that if a foil or target is heated to 
a high tejjperature , hydrocarbon molecules in gaseous form my be 
decomposed upon random collisicJiid-th the spot heated by the bean, 
even if there is no tendency othervase to adhere as a semi-pemanent 
or permanent film. The sources of the volatile hydrocarbons are seve- 
ral, the most probable ones being: 

(a) vapors from the oil diffusion p'ornp of the vacuiun system; 

(b) volatile vapiors from organic substances such as vacuum rings and 
gaskets in the systoa; 

(c) vaporized solvents from freshly applied sealing paint (such as 
^lyptal) used to stop small leaks; 

(d) vapors from small amounts of acetone or similar penetrating 
fluids used for leak testing; 

(e) vapors from residual araounts of organic solvents used for 
cleaning parts of the cysten. 

These factors ere listed in such an order that the first ones are 
most difficult to avoid, and the latter ones the easiest. Fortu- 
nately, our experience seems to indicate that the very ones hardest 
to avoid are actually the ones having the least noticeable effect, 
and vice versa. 

Under normal operating conditions it vra.s our practice to discard 
a foil or target if the blackening became appreciable. This t;as 



es^eciaLly necessary for the foilj evidently tlie greater diffuseness 
of the bean on the targets and their greater ability to conduct the 
heat avrey prevented rapid “carbonisation” of the targets. 

In connection xri.th the foils, this carbonizing added sraall 
anounts to the stopping. The effect v:as checked by detemining 
energy loss through the foils at intervals throughout their use. 
Figure 10, i;hidi gives tlie energy loss in the various foils as 
function of the enteniig proton energy, dioxra this effect clearly. 
Note specially the progressive increase in stopping pover by the 
,05 nil copper foil for subsequent days (on the order of 2 kev per 
operating day for 1 llev protons). Also note the radical shift in 
energy loss through the ,05 nil nickel foil (#2), vihich occurred 
after a rather severe carbonizing due to release of gaseous liydro- 
carbons into the gas dianber from incojiqjlctcly dried applications 
of glyptal. 

Foil stopping pover data were discarded when it vas evident that 
such carbonizing had contributed to abnormally hig)i values of energy 
loss. This data is needed, however, when the carbonized foils are 
used as windox-ra for the esrocriniGnts on gas stopping power. It shoiild 
be nentioned tliat the greatest problem did ocenr in the gas stopping 
power experiments. In the foil stopping power experiments, the liiole 
system vra.s continually pumped by the vacum pui^sj bat in the gas 
experiments, the gas chamber x;as closed off and any hydrocarbon mo3.e- 
cules trapiaed therein had nai^ chances to strike the hot foil and be 
decon^josed and deposited there. The examples of severe carbonizing 

-~/^T - 



nentioned in the preceding paragraph occurred during cjq>erinients on 
gases. 

L* AccT-TTatc Conpiitation of Mean EnerCT' of Protons in Iledim . 

It is u Si-tally custonary in giving experiiiental results for 
stopping poi/cr to consider that AE/^x equals dE/dx for a value 
of E liiich is tlie arittoetic average of the initial and exit ener- 
gies of the bean in the storming mediun ( ~ Ep i- Eg 

2 

esanple, in case the stopping naterial is a foil). This approxi- 
mation is aijpropriate when AE is snail corroared to E. Uarshaw 
(1949) lias commented on this natter and proposes a fornula to give 
a correction to AE/ Ax so that the correct value for stopping 
can be obtained at the aritlmetic mean value of energy. His formula 
assumes that the stopping pouer curve over the region from Eq to 
Eq* is a straight line. He finds that the corrections thus con- 
put cd are less than his probable error but vrorth including to insure 
a correct shape for Ms curves, 

CX-ir values of both AE and E are in most cases larger than 
Warshav/'s, and we may also find a slight correction advantageous. 
Rather than Warshau’s fomula correcting stopping power ue prefer one 
irhich provides a correction to the aritlmetic mean energj'' to give the 
true nean (designated as E^) • In deriving the fomula for this 
correction let us nalce an assui:: 5 )tion as to the shape of the stopping 
poi-7er curve i/hich is more accurate than the linear shape assumed by 
Warshair, Assume that 

^ = k E“"^ . (lV-A-3) 

dx 
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Theory Indicates, as vie have noted before, that the exponent V" 
should have a value of 0,5 (See Section III-G) , Experimental re- 
sults (Madsen and Venkatesv/arlu, 1948a; Kahn, 1953) show that Y* is 
slightly dependent on Z and on E , but has in the cases wherein 
we are required to use it (E '^600 kev) an average value of about 



0.45 . 

The thickness of the medium traversed, ^ x , is easily seen 
to be given by: 

I E^- dE (IV-A-4) 

A 



f ° 

dE 

I dE/dx 



k (r+l) 



(eA 



r-t-l 



- E, 



r+i 



(IV-A-5) 



Now the correct value of mean energy, E , to use is that for which 

m ’ 






Ae/Ax 



Ax 



(IV-A-6) 



Solving for Ax and substituting equation (IV-A-3) gives: 

E* - E 

o o • 

Ax = 3y* (IV-A-7) 



k E 



m 



If we equate the tv;o expressions for Ax and solve for , we obtain: 



E 



m 



z f— i_ 
I 1 



1 Vr 

- E^ ' * 



E’ - E^ 
o o 



( IV-A-8) 



If we replace E» by E„4* -^AE , and E. by E„ - -^AE , where 

V/ Q, O G 

/r 



Eg is the arithmetic average, we find 






y-ti 



r-h 1 



(E^ +• ^AE) - (E^- if 

AE 






(IV-A -9) 






If we expand th.e parenthetical expressions according to the 



bincnnial fomnila, we find: 

\ "*■ ^ + other terms | ^ 

(IV-A -10) 

where the other terms are successive terms in (^E) /e“” , 

for n ■ 2, 4, 6, etc., with fractional coefficients for the terms. 

For appreciably less than E , this series converges rapidly, so 

that the other terms in the above equation may be neglected. 

Then let us expand the bracketed expression in equation (IV-A-10) 
once again by the binomial formula, to obtain: 



+ additional terms , (IV-A-ll) 
m a 24 E —— - ’ ' 

a 

where the additional terms can once more be dropped x<jhen Ae is 
appreciably less than E • 

Then we see, substituting 0,45 for V , that 



E„ * E^ - i E 
in a 

where 



(IV-A-12) 



<fE 



,55 (AE )^ 
24 E 



(IV-A-13) 



a 

V/e Tise this correction and apply it below in cases where it 
exceeds 1 kev . 

5, Corrections Required for Impurities in the Materials Stud led. 

One possible source of consistent error vrtiich should be care- 
fully checked is the impvirities in the elements studied. 

The metals are quickly disposed of. A semi-quantitative 



- 



1 



;ri.th sarnie foils fron 



spectro^jraphic analysis uas cari'ied out 
the sane set as those used, vdth the foilovdng resiilts: 

(a) For copper, the ir/purities present are as follo\7s: calci’jn 

\d.th the order of mgnitude of 0«l/a, nagncsiun less than 0 , 1 %, and 
a trace of silicon. 

(b) For nickel, the inpurities present are as foHoi;s: cobalt to 

an cL3.x)imt betv/een 0.1 to plus traces of nagnesiun, nolybdenun, 

berylliun, and iron. 

The error in stopping po;;er is not as great as the percentage of 
iE 5 )uritles. Each of the inpurities tends to produce an error only 
to the ortent that its stopping poucr varies from the stopping poxrer 
of tile prine constituent, so that the closer the iir^urity is to the 
prinary elonent in the table of elements, the less the effect of the 
ir^urity on the stopping power measured. With these considerations 
in mind, it becomes clear that the impurities in foils of either 
nickel or copper are such in quantity and type tliat no appreciable 
error can result T::!y considering the foils as pure, 

2 

The gases however rcqid.re more careful study. The gases as 
they arc supplied bj'- the company are guaranteed to be pure to the 
folloidng extent : 



The analysis iras carried out hy llr, San V7ohlfort at the Spectro- 
graphic Laboratoiy, Department of Chonistry, Ohio State Universi- 
ty, tlirough the courtesy of Professor W.ll. liacNevino 



The nitrogen and argon i/ere obtained in pressurised K cylinders 
from the Linde Air Products Company; i/hllc the rarer gases came 
from the cane company in glass flasks at approximately atmo- 
spheric pressure. 






(a) CoEirjercial nitroa'ai used is stated to be 99 , 7 % pure, 

(b) The arson sx 5 >plied is stated to be 99 , 6 % pure. 

(c) Tile neon irpurities are negligible, 

(d) The Icrypton nay have as mcli as 1 % xenon present, 

(e) The xenon nay have as nuch as 7 .% Icrypton present. 

Since the above amounts are naximun values, the exact extent 
of the irpurities for our samples uas still sonev/hat in doubt. 

There existed a still more serious factor in connection with the 
latter tliree gases, which were supplied in one liter or one-half 
liter glass flaslcs. Difficulties were encountered in the attach- 
ment, of the flasl^s to the system in such a way as to be absolutely 
air-tight idien the duct system is under partial or cooqjlete vacuum. 

The results were that a certain amount of air contamination was 
probable for the gases in eacli of the three flasks, before final 
gas-tightness was attained. It imis considered necessary therefore 
tliat an analysis be made on the gases in the flasks after the pos- 
sible contamination occurred. 

The analysis was carried out by mass-spectrograph^ , and results 
are as given in Table XVI, 

In order to determine the effects of these iD5>urities, it is 
necessary to note that the amount of matter traversed is indicated 
ly the measured values of path length, gas pressure, and gas tenpera- 
ttirc. Since the pressures used are quite lovf (a few centimeters of 

1 

The anal^’-ois was carried out by lir. Leonard I'auk, at the Kass- 
Spectrographic Laboratary, Department of Chemistry, Ohio State 
University, through the coiurtesy of Professor II, L, Johnston and 
Dr, David Vihite, 






raercury) the ideal gas relationships are valid, and the variables 
cited above are sufficient to give tho nunber of nolecules per square 
centimeter of cross-section perpendicular to the proton path. For 
the rare gases, the atonic stopping povrer is the sarae as the molecu- 
lar stopping po\;erj and since the irpuritics are given as volvme 
percent (equivalent to mole percent), the corrective factor to be 
applied to the e^qperinental results is the relative moleculqr stopping 
power of the iiqjure gas as con^ared to that for the p'ure gas. 

Table XVII giTOs tho values assumed for atomic and nolecular 
stopping povrer for the various elements and compounds concerned. It 
has been necessary in some cases to make rough preliminary estimates 
of stopping povrer from the data taJeen in our experiments, so that 
the values talcen can be justified only ty the final ansi/er. Fortu- 
nately, hi^ly accurate values are not necessary in this particular 
set of computations. The sources of data arc indicated in the table. 

Table XVIII tabulates the calculations of the corrective factors 
for tho gases neon and xnnon. It is quite evident from the analysis 
of the gases in Table XVI that no appreciable leaJiage into the system 
occurred for the Inypton experiments, so that no correction is r*e- 
quired for this gas. The correction factors computed must divide the 
ejqDerimental rcsulfcs to obtain final corected results. 

We have neglected oxiy effect of possible impurities in nitrogen 
and argon. This is justified for t\;o reasons. First, the gas cham- 
ber was supplied v;ith nitrogen or argon directly from the pressurized 
tank, so that the filling system up to the inlet valve of the chamber 
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Table 'XVI 



Impurities in Volume Percent in Gases, As 
Determined by Lass Spectrograph 



impurity 


Primary Gas 


Leon 


Krypton 


Xenon 


HgO 


2,72^. 




0.02f 




8.15 


o.W. 


3.73 




1.11 


.011 


.72 




.102 


- 


.017 


COp 


.59 


.005 


- 




1,19 


- 


,172 


Ne 


— 


.026 


.07 



Table xVII. 

Approximate Values of 7 \tCBiic and l olecular Stopping 
Fo^ver of Varioxis Gaseous Constituents for lh*otons at 
Ji’nergies of 500 and 1000 kev 



Atom or 
f olecule 


Atonic 

Stopping } Oiver 
(ev-cm^ X 10^‘^) 


I'olecular 
Stopping Pov/er 

(ev-cm^ X 10^^) 


Ref. 


500 kev 


1000 kev 


500 kev 


1000 kev 


ir 


1.1 


0.6 


2.2 


1.2 


(2) 


c 


7.25 


4.35 


- 


- 


(2) 


N 


8.0 


4,9 


16,0 


9,8 


(1) 


0 


8.8 


5.5 


17.6 


11.0 


(3) 


I'Je 


10.1 


6.45 


10.1 


6.45 


(3) 


A, 


14.4 


9.8 


14.4 


9.8 


(1) 


Kr 


21.1 


14.2 


21.1 


14.2 


(1) 


Xe 


24.2 


16.25 


24.2 


1C. 25 


(1) 


IlgO 


- 


- 


11.0 


6.7 


(4) 


COg 


— 


— 


24.85 


15.35 


(4) 



(1) Estimated from "reliminary rcs\ilts, nr^nent research. 

(2) iixtrapolation from data given by Livingston and Bethe (1937), 

(3) Interpolated value, according to atonic number, from values 
obtained by references (1) and (2), 

(4) Obtained by addition of atonic stopping pov;ers. 
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Table X\riII. 

Computation of Inpurity Correction Factors, to Convert Erpertmentally 
Detcrnined j.tonic Stopping Fov/er to the Correct Value - Valid for 
Protons of 500 and 1000 ker . 



L’olecule 


^mol 


(xlO^^) 


Relative 


1 

mol 




Vol. 

55 


(Vol.^) 


X (Rel.v-) 




500 kev 


1000 kev 


500 kev 


1000 kev 


500 kev 


1000 kev 


(a) iJeon 


















I'Jeon 


10.1 


6.45 


1.00 


1.00 


86.15 


86.15 


86.15 


IlgO 


11.0 


6.7 


1.09 


1.04 




2.7 


2.94 


2.81 


% 


16.0 


9.8 


1.58 


1.52 




8.15 


12.9 


12.4 


^2 


17.6 


11.0 


1.74 


1.70 




1.1 


1.91 


1.87 


A 


14.4 


9.8 


1.43 


1.52 




.1 


.14 


.15 


CO 2 


24.85 


15.35 


2.46 


2.38 




.6 


1.48 


1.43 


Kr 


21.1 


14.2 


2.09 


2.20 


1.2 

100.0 


2.51 

103.03 


2.64 

107.45 


(b) Xenon 


Xenon 


24.2 


16.25 


1.00 


1.00 


95.3 


95.3 


95.3 


Ne 


10.1 


6.45 


.42 


.40 




.1 


.04 


.04 


N2 


15.0 


9.8 


.66 


.60 




3.7 


2.44 


2.22 


O2 


17.6 


11.0 


.73 


. 68 




.7 


.51 


.48 


ICr 


21.1 


14.2 


.87 


.87 


.2 

100.0 


.17 

98.46 


.17 

98.21 



Suranary; To correct experimental data, divide by factors: 



Prime 

Stopping 

Substance 

neon 

Xenon 



Proton Energy 
500 kev 1000 kev 



1.080 

.985 



1.075 

.982 



i;as under a positive pressi’-re relative to atnospheric. This pre- 
vented inuard leaking of air, so that most of the In^urlties pre- 
sent could be attributed to those orisinally present in the gas 
containers purchased. In re<;;ard to these original irajurities, v;e 
note tliat the percenta;^-e allowed is snail. Purthemore such im- 
purities are priroarily air constituents. Table XVII shows that 
the molecular stopping ix>wer of nitrogen, oj^ygen, and argon are all 
very close together, so that small amounts of one gas in the other 
as irpurity would laavc an extremely small effect on stopping po^;er. 

6. Possible Errors D’ue to Air Leakage During Experiment . 

It ins discovered during most of the cjncrinents on the gases 
that the gas chamber pressiu'e increased slov;ly, by varying amotmts 
at different times. At the worst, this change amounted to about 25 % 
of the initial pressure. Tliis indicated a leak of seme kind into 
the gas system during operation. Continual monitoring of the pres- 
sui'O indicated that such leaking occurred at a fairly constant rate 
for a single fil.ling, although the rate varied from one day to the 
ne:di. The errors due to such lealcing could be alloxred for therefore, 
provided the type of gas leaking vrere Icnoiai. It night be possible 
that the lealcages were either of air, or of the gas being studied, 
throu^ the filling system. The possibility also existed that a 
combination of both effects in unJeno^m proportions occuri'ed. 

Great precautions v/cre ta3:cn to reduce this effect to a minimum, 
and frequent observation ira.s made on the rate of leaJeage into the 
chamber ;ath no gas present, Ebccept for a couple of times ■vilien 
definite air-leaicgge t:as sliown (data taicen on these days shovfcd radi- 
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cal depart' ’.res from other data on the same gas, and \;ere discarded), 
no appreciable pressi-ire cliange vras discovered in the e'vac'aatcd and 
closed gas chamber over periods of several honrs duration. 

Despite these checks for air leaks, the effect persisted during 
periods that gases bfeing steadied uere in the system. After some dis- 
cussion v/ith other persons in the departments of physics and ctaenistry 
vdio vrorked -idth high vacuums, it \ras concluded that the effect v/as 
due to the entrjyDiiient of minute bubbles of gas in the stop-cock 
grease at the stop-cock which occiuved during 'the filling process. 

Such bubbles gradually traveled to the low pressure port and dis- 
charged the gas into the chamber. The prominence of such effect de- 
pends upon rijny factors - type of stop-cock used, type of vacuum 
grease used, ameunt of -vacuum gi'ease applied, pressure \/ith which the 
valve of the stop— cock is forced onto its seat. It is concluded that 
the main difficulty v;as due to use of a type of grease used for oi'di- 
nary chemical applications, but not esj-^ecially s.iited for high vacuum 
v;ork. For futiTe ;ror3cers, Apiezon, Type II, grease, applied strictly 
in accordance xdth instructions on the label and in sparing quantities, 
is recommended. 

lIo't\d.thstanding the failure to reduce the effect of pressiare 
change to zero, the c^gjerinents made and precautions talren insured 
that a negligible proportion \ias due to air leakage. All that is 
necessary for o'ur analysis is to loiov; the precise value of gas pres- 
sure in the chamber at the tine that a resonance peak is neas'ured. 

Since the pressure in the system at varlo'as times during the se- 
quence of runs v:as recorded, and follov;ed very closely a linear 
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change \vith tine, precsares and variables dependent upon tine can be 
detemined through a single linear fomula for each gas filling. 
Tliese formulas will be given in the conpilation of data \diere they 
apply. 

B. Stopping Power Results - Foils . 

1. Determination of Ener^ of Displaced Peahs, Foils Inserted . 
Tables XIX and XX shov; the tabulated data and results of deter- 

nlnation of displaced peak energy values (E^ ). The values of l’ , 
the displaced peaJc position on a nagnet current scale, are given ly 
the yield plots (Figure 8(b), for The values of E^' are 

then detemined fron the calibration corves (Figure 9). The proper 
calibration curve to use on any particular day is determined by one 
or none calibration points taken in the sane sequence of runs. This 
is discussed in detail in Section III.H, above. 

2. Comm tat ion of Foil Stopping Poorer . 

Tables XXI and XXII indicate in tabulated fom the cor^utations 
for stopping poxjer of nickel and copper, respectively. The defini- 
tions and methods of conputing the variables specified in the column 
headings are already completely set forth. Stopping po\/er is com- 
puted both in the "experimental" units and also, aa atcsnic stopping 
power, in those tinits closer to basic theory, A coliinn is provided 



stopping of slov7 particles in heavy nedia provides that stopping 



for tabulation of values of 




Since Bolir's theory for 



power is proportional to Z^/2 (Equation I-B-74)» such a tabulation 



gives a better means of comparing the results of stopping power ex- 



Table XIX 



Detemination of Shifted Resonance Peak Positions in the Energy 
Spectrun, with Various Uickel Foils Inserted into the Proton Beam. 



Foil 


Date* 


Eo 

(I'ev) 


Shifted Peak Position 


I* 

(amp) 


log I» 


^'log E» 


(Mev) 


.05nll Hi-1 


8-29-5E 


.6714 


.7848 


1.8948 


1.9241 


.8397 


n 


tf 


.8756 


.8696 


1,9393 


0,0115 


1,0268 


ft 


9-26-52 


.9875 


.9150 


1.9614 


0.0517 


1.1264 


,05mil Ni-2 


1-19-53 


.6696 


.7908 


1.8981 


r.9245 


.3404 


n 


tt 


.8740 


.8776 


"1.9433 


0.0107 


1.0249 


n 


tt 


.9358 


.9015 


1.9550 


0.0330 


1,0790 


#' " 


1-21-53 


.6696 


.7955 


1,9006 


1.9274 


.8461 




ft 


.8740 


.8823 


1.9456 


0,0130 


1.0304 




tt 


.9358 


.9075 


1.9578 


0,0360 


1.0864 


.04 mil Hi 


10—16—52 


.6714 


.7720 


1.8876 


r.9084 


.8098 


ft 


ft 


.8756 


.8620 


1,9355 


0,0027 


1.0062 




tt 


.9374 


.8860 


1.9474 


0.0261 


1.0619 


n 


10-17-52 


.6319 


.7560 


1,8785 


1.8905 


.7771 


ft 


tt 


.9875 


,9069 


1.9576 


0.0435 


1.1053 


ft 


10-23-52 


.4441 


,6681 


1.8248 


1,7870 


.6123 




10-28-52 


.9374 


.8868 


I. 9478 


0,0270 


1.0641 


H4 " 


10-31-52 


.9374 


.8870 


1.9479 


0.0283 


1.0673 


.04 mil Hi 


12-19-52 


.9875 


.9145 


1.9612 


0.0440 


1.1066 


ft 


1-9-53 


.6696 


.7765 


1.8901 


1,9080 


.8091 


n 


It 


.8740 


.8650 


1.9370 


1,9986 


.9968 


.03 rail N1 


1-22-53 


.6696 


.7519 


1.8762 


1.8840 


.7656 


It 


It 


.8740 


.8427 


1.9257 


1.9805 


.9561 


n 


It 


.9358 


.8707 


1.9399 


0.0079 


1.0184 



Rotes: * The dates not only serve to identify the various sets of 

data extracted from plotted yield curves, but also 
indicate the proper calibration curve to use (Fig. 9), 

^ Obtained from calibration curves (Fig. 9) , with use of 
Tables XIV and XV. 

iHr Appreciable carbon accumulation indicated. 






Table 7X 



Detemination of Shifted Resonance Peelc Positions in the Energy 
Spectrum, with Various Cooper Foils Inserted in the Proton Beam. 









Shifted Peak Position 


Foil 


Date 


So 

(Hev) 


I' 

(amp) 


log I* 


log E5 


SS 

(Mev) 



mil Cu 


12 - 2-52 


.3507 


.6310 


1.8000 


1.7372 


.5460 


ft 


II 


.6763 


.7860 


1.8954 


1.9215 


.8346 


tl 


II 


.8799 


.8735 


1.9413 


0.0080 


1.0186 


fl 


II 


.9415 


.8970 


' 1,9528 


0.0298 


1.0710 


?f 


12 - 3-52 


.6319 


.7705 


1.8868 


1.8996 


.7936 


ft 


fl 


.9875 


.9218 


1.9644 


0.0485 


1.1182 


fl 


18 - 18-52 


.4419 


.6760 


1.8299 


1.7956 


.6246 


If 


fl 


. 6696 


.7834 


1.8940 


T .9188 


.8295 


IT 


If 


.8740 


.8710 


1.9400 


0.0073 


1.0170 


II 


fl 


.9558 


.8960 


1.9523 


0.0307 


1.0732 


n 


1 - 12-53 


.6696 


.7860 


1.8954 


1.9203 


.8323 


n 


(PEI) 


.6696 


.7855 


1.8951 


1.9210 


.8337 


mil Cu 




.4419 


.7580 


1.8797 


1.8877 


.7721 


ft 


fl 


.6696 


,8485 


1.9287 


1.9819 


.9592 


fl 


ft 


.8735 


.9214 


1.9644 


0.0515 


1.1259 



# Appreciable carbon accumulation indicated. 
See also the notes in Table XIX, 
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Toble ( Ti'.’O Pages) 



Stopping Power Conputation for Nickel 



Line 


Foil 

Number 


Foil' 

Thickness 

(mg/cm2) 


Eo 

(kev) 


E6 

(kev) 


A E 
(kev) 


1 


.05mil Ni-] 


1.150 


671.4 


839.7 


168.3 


2 


ff 


ff 


875.6 


1026.8 


151.2 


3 


Tt 


ff 


987,5 


1126.4 


138.9 


4 


.05mil Ni-2 1,175 


669,6 


840.4 


170.8 


5 


ft 


ff 


874,0 


1024.9 


150.9 


6 


ff 


ff 


935,8 


1079.0 


143.2 


7 


tf 


ff 


669.6 


846.1 


176.5 


^ 8 


ff 


ff 


874.0 


1030.4 


156.4 


# 9 


ff 


ff 


935.8 


1086.4 


150.6 


10 


,04 mil Ni 


0.977 


671.4 


809.8 


138.4 


11 


ff 


ff 


875.5 


1006.2 


130.6 


12 


ff 


ff 


957.4 


1061,9 


124.5 


13 


ff 


ff 


631.9 


777.1 


145.2 


14 


ff 


tf 


987.5 


1105.3 


117.8 


15 


tf 


ff 


444.1 


612,3 


168.2 


ne 


ff 


ff 


937,4 


1064,1 


126.7 


#17 


ff 


ft 


937.4 


1067.3 


129.9 


18 


,04 mil Ni 


0.977 


987,5 


1106.6 


119.1 


19 


ff 


ff 


669.6 


809.1 


139.5 


20 


ff ' 


ff 


874.0 


996,8 


122.8 


21 


.03 mil Ni 


0,640 


669,6 


765.6 


96.0 


22 


tf 


ft 


874.0 


956.1 


82.1 


23 


ff 


ft 


935.8 


1018.4 


82.6 



(kev) 



755.6 

951.8 

1057.0 

755.0 

949.4 

1007.4 

757.8 

952.2 

1011.1 

740.6 

940.9 

999.6 

704.5 

1046.4 

528.2 

1000.8 

1002.4 

1047.0 

739.4 

935.4 

717.6 

915.0 

977.1 



Appreciable carbon accumulation indicated 



Table XXI. (Concl,) 



Stopping Power Computation for Nickel 



Line 


6e 

(kev) 


(kev) 


Stopping 

Pov/er 

(kev-cm^/mg) 


^ <7“ , At, 

Stpg.Pov;er 
{ ev-cm^ ) 


^ / 

^zl/3 


Prob. 

Error, 

p 


1 


negl. 


757 


146.3 


14.27 


4.70 


1.5 


2 


Tl 


951 


131.5 


12.51 


4.12 


2.0 


3 


n 


1057 


120.7 


11.77 


3.88 


2.0 


4 


T» 


755 


149.3 


14.18 


4.67 


1.5 


5 


ft 


949 


128,3 


12.51 


4.12 


2.0 


6 


n 


1007 


121,9 


11.88 


3.91 


2.0 


7 


Carbon 


accumulation renders invalid. 






8 


fT 


TT 


II 


II 






9 


tt 


II 


If 


II 






10 


negl. 


741 


141,7 


13.80 


4.55 


2.0 


11 


n 


941 


133.7 


13.02 


4.29 


2.0 


12 


IT 


1000 


127.4 


12.42 


4.09 


2.0 


13 


II 


704 


148.7 


14.48 


4.77 


2.0 


14 


If 


1046 


120.4 


11.74 


3.87 


2.5 


15 


1.2 


527 


172.0 


16.78 


5.52 


1.5 


16 


Carbon 


accumulation renders invalid. 






17 


If 


II 


II 


11 






18 


ngl. 


1047 


122.0 


11.89 


3.92 


2.5 


19 


I? 


739 


142.8 


13.91 


4.58 


2.0 


20 


It 


935 


125.7 


12.25 


4.04 


2.0 


21 


IT 


718 


150,0 


14.62 


4.81 


5.0 


22 


II 


915 


128.2 


12.50 


4.12 


3.0 


23 


tl 


977 


129.0 


12.57 


4.14 


3.0 



1 ^ 

Values tabulated in these columns to be multiplied by 10 . 
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Table XXII. ('IVfO Pages) 



Stopping Power Conputation for Copper 



Line 


Foil 

Nxmber 


Foil 

Thickness 

(mg/cm^) 


®o 

(kov) 


(kev) 


AE 

(kev) 


(kev) 


1 


,05 mil Cu 


i,155 


350,7 


546.0 


195.3 


448.4 


2 


ff 


Tf 


676,3 


834.6 


158.3 


755.5 


3 


ff 


n 


879,9 


1018,6 


138.7 


949.2 


4 


TT 


ff 


941,5 


1071.0 


129.5 


1006.2 


5 


TT 


fi 


631,9 


793,6 


161.7 


712.8 


6 


TT 


ff 


987,5 


1118.2 


130.7 


1049.6 


7 


Tf 


ff 


441.9 


624.6 


182.7 


533.2 


# 8 


ff 


n 


669,6 


829.5 


159.9 


749.6 


# 9 


TT 


ff 


874,0 


1017.0 


143.0 


945.5 


#10 


T? 


ff 


935,8 


1073.2 


137.4 


1004.5 


#11 


n 


ff 


669.6 


832.3 


162.7 


751.0 


#12 


ff 


ff 


669.6 


833.7 


164.1 


751.6 


13 


0,1 mil Cu 


2,20 


441,9 


772.1 


330.2 


607.0 


14 


ff 


ff 


669,6 


959,2 


289.6 


814.4 


15 


TT 


ft 


‘ 873,5 


1125.9 


252.4 


997.8 



# Appr^iable carbon accunulation indicated. 
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Table XXII, (Concl«) 



Stopping Pov/er Computation for Cooper 



Line 


iE 

(kev) 


Em 

(kev) 


Stopping 

Power 

(kev-cm^/mg) 


^ <r , At. 
St pg. Power 
(ev-cm^) 


tf , 


Prob. 

Error, 


1 


1.9 


446 


169.0 


17.86 


5.815 


1.5 


2 


negl. 


755 


137.0 


14.48 


4.715 


2.0 


3 


n 


949 


120.0 


12.68 


4.13 


2.0 


4 


fi 


1006 


112.0 


11.83 


3.85 


2.0 


5 


ft 


713 


139.9 


14.78 


4.81 


1.5 


6 


It 


1050 


113.1 


11.95 


3.89 


2.0 


7 


1.4 


532 


158.0 


16.70 


5.44 


1.5 


8 


Carbon accumulation renders 


invalid. 






9 


n 


ff 


ff 


tf 






10 


fi 


ff 


ff 


ft 






11 


ff 


ff 


ff 


ff 






12 


fi 


Tf 


ff 


ff 






13 


4.1 


603 


150.1 


15.87 


5.165 


1.5 


14 


2.4 


812 


131.6 


13.90 


4.525 


1.5 


15 


1.4 


996 


114.7 


12.11 


3.94 


1.5 



-15 



jf Values tabulated in these columns to b e multiplied by 10 
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perincnts on different elements. 

3« Discussion of Errors in Exporiraental Data and Results . 

Probable errors in the individual values of stopping power ob- 
tained are difficult to judge accurately, but even an approximate 
idea of such deviations would be useful. 

Values obtained from the literature for the resonance energy 
values used are probably correct to about i 1 kev on the average. 
Some are laiovn. vdth better accuracy, since they are quoted to tenths 
of a kev. liany other peak values are determined in relation to 
those few peaks x/hose energy have been accurately obtained ty abso- 
lute methods, and thus their error may be much larger, fhe value of 
1 kev is considered a rough average for all those used herein. 

Errors in target thickness are lilcely to be of the same order 
of magnitude. Ho^;evcr, although such errors affect the values taken 
for Eq , they cause at the same time an erroneous shift in the cali- 
bration exarve used vrhich introduces an error of the same order of 
magnitude in the value of Eq'. Thus for AE , these two errors 
largely cancel. V/c xrf.ll therefore neglect such errors. 

The irqjortant errors to be considered in Eq' are due to error 
in obtaining the proper value of l' (-value of shifted peak position 
in the magnet current spcctr-um) and to an error in choice of the cali- 
bration curve used. In regard to the foincr it must be noted that 
the stopping phenmenon not only shifts the peak but also broadens 
it, so that the position of the maximum point is more difficult to 
estimate. Whereas the xmde-viated peak can be judged to an accuracj' 
of about 0.1 milliaiqDores, or about 0.2 kev, the de-viated peak lias a 



greater probable error - probably on the order of 1 kev. 

In regard to the error from calibration inaccuracies, let us 
refer to Section III.H v/herein a typical value of probable error 
is about 1,7 kev by use of the mean curve. We can improve on this, 
hoxjever, since \;e determine one or more calibration points during a 
given day's sequence of operations, A value of 1 kev should be 
about the proper value for this error. 

Summarizing the above, ve see that the probable error in Eq' 
is about 1,4 kev and about 1 kev for Eq, so that the probable 
error from the diffei'cnce of these tvra is about 1,7 kev. We feel 
it necessary to add to this an estimate for surface impurities, 
primaiily carbon, of about 1,4 kev (see Section IV,A,3.), This 
gives a final probable error in A E of about 2,2 kev. 

The percentage error i/hich tliis figure provides depends on the 
value of AE, Reference to Tables 3QCI and XXII indicates that the 
probable error varies from 0,7^ for thidc foils and lo\; proton ener- 
gies to 2,7% for thin foils and high proton energies, ^o take into 
account the possible increase in path length at lou energies, let 
us raise the loiter figure to 1 kev. 

It has already been shorn in Section III.B, tlmt the probable 
ein*or in foil thicknesses is about 1,1%, This means that the 
probable error in values for foil stopping pov;er varies between 1,5% 
and 3%« Tables XXI and XXII tabulate the estimates of probable 
error for individual values of stopping povrer. 

It must be emphasized that the probable errors given are based 
on the assimiption of all errors conforming to the Gaussian theory of 
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errors, \/ith positive and negative errors equally likely. Further- 
more this analysis does not guarantee that unknown errors of either 
a consistent or random nature have not caused the tabulated values 
of probable error to be less than they should be, 

C. Stopping Po\/er Results - Gases . 

1« Energy Loss in Windov; Foils . 

It is evident that in obtaining the energy loss tlirou^ the gas 
in the chamber, it is necessary not only to subtract Eq from Eq* , 
but from this difference to subtract the energy loss in the \/indow 
foil (Equation II-B-9). Since for vri.ndov;s v;e used the some foils 
on which the stopping power measurements were made, we can determine 
the value of vrindov; foil energy loss, \jj making use of a 

chart whereon is dravm for each foil used a curve shoxdng foil 
energy loss versus proton energy at the foil entrance, E^’. Using 
the data provided in Table XXI and XXII, this can be readily accom- 
plished, Figure 10 gives these curves. It shoTild be noted that 
certain points are plotted vAiich indicate clearly the effect of 
carbon accumulation. Such data enable one to malce a better estimate 
of Bp so as to iuiclude carbonizing when appropriate. 

2« Energy of Displaced Peaks, Gases and Window; Foils Inserted . 

These values are determined precisely as they are for the case 
of foils alone (Section IV.D.l). Tables XXIII-XXVII, inclusive, 
show this determination, 

3« Determination of the Formula for Atomic Stopping Power of Gases . 



The fundamental conversion formulas for stopping power have al- 



ready been evaluated. It is necessary hov/ever to develop a formula 
which idll give atonic stopping povjer as a function of gas tempera- 
ture (T), pressure (p), and bc^ energy loss, I’he other factor, 

t 

path distance, is alreacfe^ knoxna from measurenent of the length of 
the gas chamber (Fig.re 4)» This distance, Ax, is 10,56 cn. From 
Section I, A,, vie see tliat 



Cr (ev-cm2) - . AE , (lV-C-1) 

vAiere A E is measured in kev. 

The gas pressures used are so lov; that the ideal gas relations 
are sufficiently accurate for our purposes. Then, under standard 
conditions, 

IIcTP “ 6«023 X 10^^ • n , (IT— 0-2) 

22,414 



vAiere n is the n.nber of atoros per molecule of the gas considered, 
Furthemore , 



~ ^ ^ ‘ ^^STP » (IV-C-3) 



vrfaere p is measured in centimeters of mercury, and T is measured 
in degrees Kelvin, 

Combination of the above formila gives the relation sought : 



<r = 



G - AE 






vAiere 



(IV-C-4) 



J 



1,021 X -p-n 

Since p and T are functions of time (Section IV, A, 6), it is possible 



to determine G as a function of time for each gas filling. Table 



Table XXIII 



Detemlnation of Shifted Resonance Peak Positions in the 
Energy Spectrum, with Nitrogen Gas and VJindov# Foil Inserted 

in the Beam. 



W Indow 
Foil 


Date 


®0 

(Mev) 


Shifted Peak Position 


I* 

(amp) 


log I* 


"log 


(Kev) 


,04 mil Ni 


. 10-27-52 


.3438 


.6820 


1.8338 


1.8040 


.6368 


tf 


ff 


.4441 


.7237 


1.8596 


i.8543 


.7150 


ft 


ff 


.6006 


.7846 


1.8946 


1.9225 


.8366 


If 


ff 


.6714 


.8132 


i.9102 


1.9525 


.8964 


ft 


tf 


.8756 


.8936 


1.9511 


0.0318 


1.0760 


ff 


ff 


.9374 


.9164 


i.9621 


0.0533 


1.1306 


.04 nil Ni 


10-28-52 


.3438 


.6852 


1.8358 


1.8085 


.6434 


ft 


fi 


.4441 


,7255 


1.8606 


1.8563 


.7183 


ft 


ft 


.6006 


.7856 


1.8952 


1.9234 


.8383 


ft 


ft 


.6714 


.8150 


1.9112 


1.9544 


.9003 


M 


ft 


.8756 


.8940 


1.9513 


0.0321 


1.0767 


ft 


ff 


.9374 


.9168 


1.9623 


0.0535 


1.1311 


.03 mil Ni 


1.23-53 


,3413 


.7050 


1.8482 


1.8295 


.6753 


ff 


ft 


.4419 


.7413 


1.8700 


1.8716 


.7440 


ff 


ff 


.6696 


.8269 


i.9175 


1.9622 


.9166 


ft 


ff 


.8740 


.9045 


1.9544 


0.0368 


1.0884 


ft 


If 


.9358 


,9279 


1.9675 


0.0580 


1.1429 



Obtained from calibration curves (Fig. 9) , xvith the use of 
Tables XIV and XV. 
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Table XXIV 



Determination of Shifted Resonance Peak Positions in the 
Energy Spectrum, v/ith Argon Gas and V/indow Foil Inserted 

in the Beam, 



IVindovf 

Foil 


Date 


O <D 


Shifted Peak Position 


I* 

(amp) 


log I» 


*log Efc 


% 

(Mev) 


.04 mil Ni 


10-29-52 


.3438 


,7492 


1.8746 


1.8837 


.7651 


ft 


ff 


.4441 


,7852 


1.8950 


i.9238 


.8391 


ff 


n 


.6006 


.8427 


1.9257 


1.9840 


.9638 


n 


ff 


.6714 


.8668 


1.9379 


0.0077 


1.0179 


ft 


ff 


.8756 


.9370 


1.9717 


0.0740 


1.1858 


.04 mil Ni 


10-30-52 


.3433 


,6978 


1.8437 


1.8235 


.6660 


If 


ff 


.4441 


.7415 


1.8701 


1.8750 


,7499 


If 


If 


,6006 


.7990 


1.9025 


1.9385 


.8680 


If 


ti 


.8756 


.9054 


1.9568 


0,0450 


1.1092 


ft 


n 


.9374 


.9280 


1.9675 


0,0658 


1.1636 


.04 mil Ni 


10-30-52 


,8756 


,9080 


1.9581 


0.0474 


1.1153 




p.m. 













Obtained from calibration curves (Fig, 9) , with use of Tables 
XIV and XV* 



— /"K 









I 



Table XXV 



Determination of Shifted Resonance Peak Positions in the 
Energy Spectr\nn, with Neon Gas and Window Foil Inserted in 

the Beam. 



- 


S 


hifted Peak Position 


Window 

Foil 


Date 


(Mev) 


I* 

(amp) 


log I* 


*log E» 


(l!ev) 


# .05nil Ni-2 


1-20-53 


.3413 


.7164 


i.8552 


1.8453 


.7003 


n ^ 


11 


.4419 


.7572 


1.8792 


1.8910 


.7780 


ft 


If 


.5986 


.8180 


1.9128 


1.9553 


.9022 


IT 


IT 


.6696 


.8440 


1.9263 


1.9810 


.9572 


If 


Tl 


.8740 


.9210 


1.9643 


0.0535 


1.1311 


If 


TT 


.9358 


,9452 


1.9755 


0.0746 


1.1874 


# ,05mll Ni-2 


1-21-53 


.3413 


.7220 


1.8585 


1.8475 


.7039 


n 


If 


,4419 


.7617 


1.8818 


1.8921 


.7800 


ft 


If 


.6696 


.8469 


1.9278 


1.9802 


OKRA 


n % 


t If 


,8740 


.9250 


1.9661 


0.0535 


1.1311 


If 


If 


.9358 


.9466 


1.9762 


0.0725 


1.1817 



Obtained from calibration curves (Fig, 9) , with the use of 
Tables XIV and XV. 

^ Foil was blackened by heavy carbon deposit during runs on 

20 January. Target thickness ivos checked before use on 

21 January (Figure 10), 






Table XXVI 



Determination of Shifted Resonance Peak Positions in the 
Energy Spectrum, with Krypton Gas and V/indow Foil Inserted 

in the Beam. 



V/indov/ 

Foil 


Date 


Eo 

(Mev) 


Shifted Peak Position 


I’ 

(amp) 


log I* 


"log Ei 


Ei - 

(Kev) 


.04 mil Ni 


1-9-53 


' .4419 


.7518 


1.8761 


1.8818 


.7617 


If 


ft 


,5986 


.8130 


1,9101 


1.9472 


.8855 


M 


If 


.6696 


.8422 


1,9254 


1.9765 


.9473 


If 


tl 


.8740 


.9206 


1.9641 


0.0515 


1.1259 


.05 mil Cu 


1-12-53 


,3413 


.7070 


1,8494 


1.8332 


.6811 


If 


IT 


.4419 


.7492 


1.8746 


1.8814 


.7610 


If 


II 


.5986 


.8110 


1.9090 


1.9478 


.8867 


If 


II 


.6696 


.8404 


1.9245 


1.9775 


.9495 


rt 


II 


.8740 


.9170 


1.9624 


0.0503 


1.1228 


If 


If 


.9358 


.9394 


1.9729 


0.0705 


1.1762 



* Obtained from calibration curves (Fig, 9) , with use of Tables 
x'rv and XV, 
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Table XXVII 



Determination of Shifted Resonance Peak Positions in the 
Energy Spectrum, with Xenon Das and Window Foil Inserted 

in the Beam. 



Window 

Foil 


Date 


^0 

(Kev) 


Shifted Peak Position 


I' 

(arap) 


log I» 


*10g 


E» 

(Kev) 


.05 rail Cu 


1-14-53 


.3413 


.7270 


1.8615 


1.8572 


.7198 




ft 


.4419 


.7728 


1.8881 


1.9080 


.8091 


If 


ft 


.5986 


.8268 


1.9174 


1.9645 


.9215 


n 


ft 


.6696 


.8561 


1.9325 


1.9932 


.9845 


If 


Tf 


.8740 


.9292 


1.9681 


0.0613 


1.1516 


• 05 nil Ni-2 


1-19-53 


.3413 


.7351 


1.8663 


1.8635 


.7303 


n 


ft 


.4419 


.7747 


1.8891 


1.9073 


.8078 


It 


ft 


.5986 


.8300 


1.9191 


1.9647 


.9219 


tt 


ft 


.6696 


.8581 


1.9335 


1.9923 


.9824 




ft 


.8740 


.9517 


1.9785 


0.0791 


1.1998 



Obtained from calibration curves (Fig. 9), with use of Tables 
XIV and XV. 
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Table XXVIII. (TV/o Pages) 

Computation of G(t) , a Factor Required in Com- 
puting Stopping Power of Gases (Equation rV-C-4) 



Date 


T®Abs, , 
Initial: 
Final : 


p (cm). 
Initial: 
Final : 


* 

Initial: 
Final : 


Hour,t, 
Initial; 
Final : 


* 

G(t) 


(a) Fitrogen 










10-27-52 


298.5 


2.02 


72.3 


- 


72.3 




299,5 


2,03 


72.3 






10-28-52 


297.4 


2.076 


70.1 


9:45 


70.1 - 0.847 (t-9. 25) 




298,9 


2.20 


66,5 


14:00 




1-23-53 


'298.1 


3,915 


37.25 


- 


37.25 


(b) Argon 












10-29-52 


297.5 


5.29 


55,1 


10:15 


55.1 -0.392(t-10.25) 




298.8 


5.60 


53.3 


14:50 




10-30-52 


297.8 


3.10 


94.0 


12:30 


94.0 + 0.20(t-12.5) 


a.m. 


299.2 


3.09 


94,8 


16:30 




10-30-52 


299.2 


3.25 


90.1 


— 


90.1 


p.m. 












(o) Neon 












1-20-53 


298.0 


4.685 


62.25 


13:27 


62.25-0.925(t-13,45) 




298.2 


4.815 


60,6 


15:14 




1-21-53 


298.0 


4.28 


68,2 


11:05 


68.2 




298.0 


4.39 


66.4 


11:25 


66.4 




298.0 


4.515 


64.6 


13:10 


64.6 




298,5 


4.85 


60,25 


14:00 


60.25 




298,5 


4.86 


60.1 


14:10 


60.1 



♦ 



Values tabulated are to be rmltiplied by 10 
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Table XKYIII. (Concl«) 



Date 


T®Abs., 
Initial; 
Final ; 


P (cm), 
Initial; 
Final : 


♦ 

G* 

Initial; 
Final ; 


Hour,t, 
Initial: 
Final ; 


* G(t) 


(d) Krypton 










1-9-53 


297.5 


2.183 


133.4 


13:55 


133.4-0.85(t-13.916) 




298.5 


2.675 


109.3 


16;45 




1-12-53 


298.0 


2.097 


■'142.0 


13:00 


142.0 - 4.16(t-13.0) 




298.5 


2.27 


131.6 


15:30 




(e) Xenon 










■ 


1-14-53 


298.5 


2.183 


133.9 


10:30 


133.9 - 12(t-10.5) 




299.0 


2.482 


117.9 


11:50 




1-19-53 


298.5 


2.204 


132.7 


13:33 


132.7-12.84(t-13.55) 




298.0 


2.482 


117.5 


14:44 




1-19-53 


298.4 


3.36 


87.0 


16:00 


87.0 



* 
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Values tabulated are to be multiplied by 10 



XXVIII shOTJs the values of G(t) for each sequence of runs in 
studying gas stopping power, 

4» Conputation of Gas Stoupinr; Poi;er . 

Tables XXIX through XXXIII, inclusive, give the corqjutations in 
tabulated form for stopping pox/er of the gases nitrogen, argon, 
neon, Icrypton, and xenon. The various factors involved have al- 
ready'’ been fully defined and e:qDlained, The results are e^q^ressed 
in the sane units as the foil stopping power, and vdth the sane end 
in view, (See Section IV, B, 2) 

5, Discussion of Errors in Gas Stopping Power Results , 

As in the case with the foils, ue should expect the probable 
error in position of the undeviated peak to be on the order of 1 kev. 
Errors in the value of the siiifted pcalc are, as before, about 1,4 kev. 
We have to subtract from the energy shift the energy loss in the 
\dndou foil. The values of the curves plotted in Fig re 10 are 
correct to rouglily 2 kev, and an additional 1,5 kev variation in the 
stopping power of the carbon fito is estimtcd. These probable 
errors conbine statistically to give an over-all probable error in 
AE of about 3 kev. This anovmts to about 1^ for large values of 
AE and up to IS for the smallest values of energy loss in the gas. 

The probable error in pressure readings is estimated to be IJS, 
and the uncertainty in the quality and quantity of ir^Durities, 
though doubtful, is taken to cause probable errors in stopping 
pox/cr of about £ov nitrogen and argon, and about 2^ for the 
rarer gases. 






Table XXIX. (Three Pages) 



Computation of Stopping Power of Nitrogen 



Line 


Foil 


Date 


Time 


(kev) 


* 

AEp 

(kev) 


(kev) 


1 


.04 mil N1 


10-27-52 


11:00 


636.8 


164.5 


472.3 


2 


ft 


ft 


12:30 


715.0 


154.5 


560.5 


3 


ft 


It 


13:45 


836.6 


141.0 


695.6 


4 


n 


If 


15:00 


896.4 


136.5 


759.9 


5 


ft 


tf 


16:00 


1076.0 


124.0 


952.0 


6 


ti 


t? 


17:00 


1130.6 


120.0 


1010.6 


7 


.04 mil N1 


10-^8-52 


10:30 


643.4 


165.5 


477.9 


8 


n 


ft 


11:15 


718.3 


155.5 


562.8 


9 


ti 


ft 


12:15 


838.3 


142.5 


695.8 


10 


ft 


1? 


13:00 


900.3 


137.5 


762.8 


11 


ft 


ft 


13:30 


1076.7 


125.0 


951.7 


12 


ff 


ft 


14:00 


1131.1 


121.5 


1009.6 


13 


.03 mil Ni 


1-23-53 


13:50 


675.3 


102.5 


57B.8 


14 


ft 


ft 


14:05 


744.0 


97.5 


646.5 


15 


ft 


If 


14:20 


916.6 


86.5 


830.1 


16 


n 


n 


14:40 


1088.4 


79.0 


1009.4 


17 


n 


tf 


14:55 


1142.9 


77.0 


1065.9 



Allopjance made for slight carbonization of foil on 27 and 
28 October (Figure 10) , 
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Table XXIX. (Cont,) 



Line 


^0 

(key) 


AE ■ 

“^0 

(kev) 


^(Eo + E») 

(kev)^ 


(kev) 


(kev) 


Correction 
Factor for 
Impurities 
(Divide) 


1 


343.8 


128.5 


408.0 


negl. 


406 


1.0 


2 


444.1 


116.4 


502.3 


II 


501 


ft 


3 


600.6 


95.0 


648.1 


If 


648 


fl 


4 


671.4 


88.5 


715.6 


If 


716 


If 


5 


875.6 


76.4 


913.8 


If 


914 


IT 


6 


937.4 


73.2 


'974.0 


If 


974 


n 


7 


343.8 


134.1 


410.8 


1.0 


409 


fl 


8 


444.1 


118.7 


503.4 


negl. 


502 


If 


9 


600.6 


95.2 


648.2 


If 


648 


ft 


10 


671.4 


91.5 


717.1 


ft 


717 


ft 


11 


875.6 


76.1 


913.6 


If 


914 


ft 


12 


937.4 


72.2 


973.5 


If 


973 


ft 


13 


341.3 


231.5 


457.0 


2.7 


452 


If 


14 


441.9 


204.6 


544.2 


1.8 


541 


ft 


15 


669.6 


160.5 


749.8 


negl. 


748 


ft 


16 


874.0 


135.4 


941.7 


fl * 


942 


ft 


17 


935.8 


130.1 


1000.8 


If 


1001 


ft 
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Table XXIX. (Concl,) 



Line 


Corrected 
Zi E 

(kev) 


♦ 

G 

(Kiiltip. ) 


^ At. 

Stopping 


* yzi/s 


Stopping 

Power 

(kev-cm^/mg) 


Prob, 

Error, 


1 


128.5 


72*3 


9.30 


4.86 


400 


4.0 


2 


116.4 


tf 


8.42 


4.40 


362 


4.0 


3 


95.0 


ti 


6.87 


3.59 • 


296 


4.5 


4 


88.5 


n 


6.40 


3.34 


275 


4.5 


5 


76.4 


If 


5.53 


2.89 


238 


4.5 


6 


73.2 


ti 


5.30 


2.77 


228 


4.5 


7 


134.1 


69.5 


9.32 


4.87 


401 


4.0 


8 


118.7 


68.8 


8.16 


4.26 


351 


4.0 


9 


95.2 


68.0 


6.47 


3.38 


278 


4.0 


10 


91.4 


67.3 


6.15 


3.21 


264 


4.5 


11 


76,1 


66.9 


5.09 


2,66 


219 


4.5 


12 


72.2 


66.5 


4.80 


2.51 


207 


4.5 


13 


231.5 


37.25 


8.62 


4.50 


371 


2,5 


14 


204.6 


Tt 


7.62 


3.98 


328 


3.0 


15 


160.5 


fl 


5.98 


3.12 


257 


3.5 


16 


135.4 


If 


5.05 


2.64 


217 


3.5 


17 


130.1 


It 


4.85 


2.53 


209 


4.0 



* Daia in this colurin to be multiplied by 10“^® , 

^ Data in these colinans to be multiplied by 10“^ • 






Table XXX, (Two Pages) 



Computation of Stopping Pov/er of Argon 



Line 


Foil 


Date 


Time 


^0 

(kev) 


(kev) 


®3 

(kev) 


1 


.04 mil Ni 


10-29-52 


10:30 


765.1 


152.0 


613.1 


2 


ft 


ff 


11:15 


839.1 


144.5 


694.6 


3 


ft 


fl 


12:30 


963.8 


134.5 


829.3 


4 


fi 


fl 


13:00 


1017.9 


130.5 


887.4 


5 


ff 


ff 


14:00 


1185.8 


120.0 


1065.8 


6 


ft 


10-30-52 


13:15 


666.0 


164.5 


501.5 


7 


tr 


ft 


14:00 


749.9 


154.5 


595.4 


8 


ft 


ff 


14:45 


860.0 


144.0 


724.0 


,9 


ft 


ft 


15:30 


1109.2 


126.0 


983.2 


10 


ft 


ft 


16:00 


1163.6 


122.5 


1041.1 


11 


ft 


” (P*n) 


17:30 


1115.3 


126.0 


989.3 



Allowance is made for moderate carbonization of the window foil 
(Figxire 10) . 



Line 


(kev) 


As = 
(kev) 


*(^+' 15 ) 

(kev) 


(kev) 


(kev) 


Correction 
Factor for 
Dapurities 
(Divide) 


1 


343.8 


269.3 


478.4 


2.8 


476 


1.0 


2 


444.1 


250.5 


569.3 


2.5 


567 


ff 


3 


600.6 


228.7 


714.9 


1.7 


713 


ft 


4 


671.4 


216.0 


779.4 


1.4 


778 


ff 


5 


875.6 


190.2 


970,7 


negl. 


971 


ft 


6 


343.8 


157.7 


422.6 


1.3 


421 


ft 


7 


444.1 


151.3 


519.7 


1.0 


519 


ff 


8 


600.6 


123.4 


662,3 


negl. 


662 


ff 


9 


875.6 


107.6 


929.4 


ff 


929 


ft 


10 


937.4 


103.7 


989,2 


ft 


989 


ff 


11 


875.6 


113.7 


932.4 


negl. 


932 


ff 



■^/3S- 




I 



I 



Table XXX, (Concl.) 



Line 


Corrected 

(kev) 


♦ 

G 

(Multip. ) 


At. 

Stopping 

Power 


^ ‘V , 
/zi/3 


Stopping 

Power 

(k 0 Y-cn^/mg) 


Prob. 

Error, 

io 


1 


269.3 


55.0 


14.81 


5.66 


224 


2.0 


2 


250.5 


54.7 


13.71 


5.24 


207 


2.5 


3 


228.7 


54.2 


12.39 


4.73 


187 


2.5 


4 


216.0 


54.0 


11.66 


4.45 


176 


2.5 


5 


190.2 


53.6 


10.20 


3.89 


154 


3.0. 


6 


157.7 


94.2 


14.87 


5.68 


224.5 


3.5 


7 


151.3 


94.3 


14.28 


5.45 


215 


3.5 


8 


123.4 


94.5 


11.67 


4.45 


176 


4.0 


9 


107.6 


94.6 


10.19 


3.89 


154 


4.0 


10 


103.7 


94.7 


9.82 


3.75 


148.5 


4.0 


11 


113.7 


90.1 


10.24 


3.91 


155 


4.0 



% 

Data in this coltmm to be multiplied by 10 



^ Data in these columns to be multiplied by 10 






Table XXXI, (Tf/o Pages) 



Conputation of Stopping pov/er of Neon 



Line 


Foil 


Date 


Time 


(kev) 


(kev) 


(kev) 


1 


•05 mil Ni-2 


1-20-53 


13:35 


700.3 


190.0 


510.3 


2 


ft 


ff 


13:55 


778.0 


179.0 


599.0 


3 


It 


tf 


14:15 


902,2 


165.5 


736.7 


4 


tf 


ft 


14:30 


957.2 


160.0 


797.2 


5 


If 


n 


14:50 


1131.1 


144.5 


986.6 


6 


ft 


If 


15:05 


1187.4 


140.0 


1047.4 


7 


ft 


1-21-53 


11:05 


703.9 


198.0 


505.9 


8 


n 


ff 


11:25 


780.0 


185.5 


594.5 


9 


ft 


ff 


13:10 


955.4 


164.5 


790.9 


10 


t? 


ff 


14:00 


1131,1 


146.0 


985.1 


11 


ff 


ff 


14:10 


1181.7 


141.0 


1040.7 



The foil was blackened by heavy carbon deposit during run on 
20 January. Carbonized target checked for stopping power and 
used again 21 January, This effect has been taken into account. 



Line 


(kev) 


= 

H - Eo 

(kev) 


(kev) 


Se 

(kev) 


(kev) 


Correction 
Factor for 
Impurities 

(Divide) 


1 


341.3 


169.0 


425.8 


1.5 


424 


1.081 


2 


441.9 


157.1 


520.4 


1.1 


519 


1.080 


3 


598.6 


138.1 


667.6 


negl. 


668 


1.078 


4 


669.6 


127.6 


733.4 


ff 


733 


1.078 


5 


874.0 


112.6 


930.3 


ff 


930 


1.076 


6 


935.8 


111.6 


991.6 


ft 


992 


1.075 


7 


341.3 


164.6 


423.6 


1.5 


422 


1.081 


8 


441.9 


152.6 


518.2 


1.0 


517 


1.080 


9 


669.6 


121.3 


730.2 


negl. 


730 


1.078 


10 


874,0 


111.1 


929,5 


ft 


929 


1.076 


11 


935,8 


104.9 


988.2 


ft 


988 


1.075 






Table XXXI. (Conel.) 



Line 


Corrected 

AE 

(key) 


* 

G 

(Mult.) 


^ <7",At. 
Stopping 
Power 


* / 


Stopping 

Power 

{key-cmVng) 


Prob, 

Error, 


1 


156.2 


62.1 


9.70 


4.50 


290 


3.5 


2 


145.4 


61.8 


8.99 


4.17 


269 


3.5 


3 


128.1 


61.5 


7.88 


3.66 


235.5 


4.0 


4 


118.4 


61.3 


7.26 


3.37 


217 


4,0 


5 


104.7 


61.0 


6.38 


2,97 


191 


4.0 


6 


103.8 


60.7 


6.30 


2.93 


188 


4.0 


7 


152.1 


68.2 


10,37 


4.82 


310 


3.5 


8 


141.3 


66.4 


9.38 


4.35 


280 


3.5 


9 


112.6 


64.6 


7.27 


3.38 


217 


4.0 


10 


103.3 


60,25 


6.22 


2.89 


186 


4.0 


11 


97.6 


60.1 


5.86 


2,73 


175 


4.0 



* IQ 

Data in this column to be multiplied by 10“-^° . 

^ Data in these colxmms to be multiplied by 10“^^ . 



-/V-/ 




Table XXXII, (TVro Pages) 



Computation of stopping Power of Krypton 



Line 


Foil 


Date 


Time 


% 

(ker) 


(kST) 


(kev) 


1 


.04 mil Ni 


1-9-53 pm 


15:10 


761.7 


146.0 


615,7 


2 


ft 


ff 


15:45 


885.5 


134.5 


751.0 


3 


« 


If 


16:05 


947.3 


130.0 


817.3 


4 


Tf 


ff 


16:35 


1125.9 


118.5 


1007.4 


5 


,05 mil Cu 


l-12-53pm 


13:25 


681.1 


178.5 


502.6 


6 


If 


ff 


13:55 


761.0 


169.0 


592.0 


7 


If 


ff 


14:15 


886.7 


156.0 


730.7 


8 


fl 


ff 


14:45 


949,5 


152.0 


797.5 


9 


fl 


ff 


15:00 


1122.8 


138.0 


984.8 


10 


ft 


ff 


15:20 


1176.2 


133.0 


1043.2 



Line 


^ 

(kev) 


ZiE = 
(kev) 


(kev) 


(kev) 


(kev) 


Correction 
' Factor for 
Impurities 

(Divide) 


1 


441.9 


173.8 


528.8 


1.3 


527 


1.0 


2 


598,6 


152.4 


674,8 


negl. 


675 


ff 


3 


669.6 


147.7 


743.4 


If 


743 


ff 


4 


874.0 


133.4 


940.7 


ff 


941 


tf 


5 


341.3 


161.3 


421.9 


1.4 


420 


tf 


6 


441.9 


150.1 


516.9 


1.0 


516 


ff 


7 


598.6 


132,1 


664.6 


negl. 


665 


tf 


8 


669.6 


127.9 


733.5 


ff 


733 


fl 


9 


874,0 


110.8 


929.4 


ff 


929 


If 


10 


935.8 


107.4 


989.5 


ff 


989 


tf 



/V*. 



Table XHCII, (Concl.) 



Line 


Corrected 

Ae 

(kev) 


♦ 

G 

(LIult. ) 


^ <r,At. 
stopping 
Power 


* / 
7^1/3 


Stopping 

Power 

(kev-cm^/mg) 


Prob. 

Error, 


1 


173.8 


122.7 


21.3 


6.46 


153,5 


3.5 


2 


152.4 


117.8 


17.95 


5.44 


129 


3.5 


3 


147.7 


114.9 


16.98 


5.15 


122 


3.5 


4 


133.4 


110.7 


14.77 


4.48 


106,5 


4.0 


5 


161.3 


140.3 


22.65 


6.87 


163 


5.5 


6 


150.1 


138.2 


20.75 


6.29 


149.5 


3.5 


7 


132.1 


136.8 


18.07 


5.48 


130 


4.0 


8 


127.9 


134.7 


17.21 


5.22 


124 


4.0 


9 


110.8 


133.7 


14.81 


4.49 


107 


4.0 


10 


107,4 


132.3 


14.22 


4.31 


102.5 


4.0 



Data In this column to be multiplied by 10” 

^ Data in these columns to be multiplied by 10“^® 



Table KiCXIII. ('IWo Pages) 
Computation of Stopping Power of Xenon 



Line 


Foil 


Date 


Time 


-^o 

(kev) 


* 

AEj, 

(kev) 


(kev) 


1 


,05 mil Cu 


1-14-53 


10:35 


719.8 


183.5 


536.3 


2 


n 


If 


11:05 


809.1 


172.5 


636.6 


3 


tf 


Tf 


11:10 


921.5 


161.0 


760.5 


4 


If 


Tf 


11:20 


984.5 


154.5 


830.0 


5 


Tf 


Tf 


11:35 


1151.6 


139.5 


1012.1 


6 


,05mll Ni-2 


1-19-53 


14:00 


730.3 


184.5 


545.8 


7 


TT 


fl 


14:15 


807.8 


174.0 


633.8 


8 


Tf 


tf 


14:30 


921.9 


161.5 


760.4 


9 


rt 


Tf 


14:40 


982.4 


155.5 


826.9 


10 


f1 


TT 


16:00 


1199.8 


136.0 


1063.8 



Allov^ance made for carbonization of window foil (Fig. 10) , 



Line 


(kev) 


s 

-Eo 

(kev) 


*(Eo+EJ) 

(kev) 


(kev) 


Em 

(kev) 


Correction 
Factor for 
Dnpurities 

(Divide) 


1 


341.3 


195.0 


438.8 


2.0 


437 


.985 


2 


441.9 


194.7 


539.2 


1.6 


538 


.985 


3 


598.6 


161.9 


679.5 


negl. 


679 


.984 


4 


669.6 


160.4 


749.8 


fl 


750 


.983 


5 


874.0 


138.1 


943.0 


IT 


943 


.982 


6 


341.3 


204.5 


443.5 


2.2 


441 


.985 


7 


441.9 


191.9 


537.8 


1.6 


536 


.985 


8 


598.6 


161.8 


679.5 


negl. 


679 


.984 


9 


669.6 


157.3 


748.2 


TT 


748 


.983 


10 


874.0 


189.8 


968.9 


Tf 


969 


.982 







r* 



■Kpn 

V 



Table XXXIII. (Concl,) 



Line 


Corrected 

(kev) 


♦ 

G 

(Mult.) 


<r,At. 

stopping 

Power 


* , 

<ry.^lz 


Stopping 
Pov;er _ 
(kev-cm^/mg) 


Prob. 

Error, 

% 


1 


198.0 


132.9 


26.3 


6.96 


120.5 


3.5 


2 


197.8 


126.9 


25.1 


6.64 


115 


3.5 


3 


164.7 


125.9 


20.72 


5.48 


95 


3.5 


4 


163.3 


123.9 


20.22 


5.35 


92.7 


3.5 


5 


140.6 


120.9 


16.99 


4.49 


77.8 


4.0 


6 


208 


126.9 


26.4 


6.98 


121 


3.0 


7 


194.9 


123.7 


24.1 


6.38 


110.5 


3.5 


8 


164.5 


120.5 


19.85 


5.25 


91 


3.5 


9 


160.1 


118.3 


18.95 


5.01 


86.9 


3.5 


10 


193.3 


87.0 


16.8 


4.44 


77 


3.5 



* —IB 

Data in this column to be multiplied by 10 



n —15 

S" Data in these columns to be multiplied by 10 



/y~s~’ 



Thus, the pi'obahle er^’or anounts to 2% to 4«5f<> for nitrogen and 
argon, and 2,5/S to h.*5% for the other gases studied. The higher 
values are valid for staaller values of energy loss, and vice versa. 
Tables XXIX throu(^ XXXIII include colunns for the probable errors 
coEg)uted in this fashion, 

D, Straggling in Stopping Power . 

1, Cornutation , 

We sh al l consider il?/ Ax as the ncasare of stopping poi/er 
straggling, and assmae it to bo valid for Ejj , the mem value of 
energy for the protons passing tlirough the stopping nediun. Results 
xvere considered only for the foils. Results for the gases are not 
easily obtained Igr using the technique erployed for our stopping 
power experinents. There are tv;o reasons fca' this: 

(1) Straggling results for vrindou foils must be knovm v/ith great 
aceiracy in order to subtract such effects from the overall results 
provided for gas plus vrindow foil. Such accuracy for the foils has 
not been obtained. 

(2) The statistical spread of the resonance peaks, idth both gas 
and windoir foil inserted, is so Icrge that adjacent pealcs interfere 
somewhat vath one another, Tliis introduces uncertainty in how much 
of the yield data are ascribable to each individual peak. 

Table XXXIV gives the data and computations for straggling effect 
in proton stopping power. The colintin headings have already been de- 
fined and discussed in previous sections. The basic formula for de- 
termination of is equation (II-B-IO), 



— /V-M ^ 



For reasons uhich have boon oxtjlained in Section III.G. , the 
values of St Ij are not usually large enough to consider. Thus 
the values of ^ considered in the final 

detemination of . 

In obtaining the values of n froa the yield 

O o 

curves the follomng rules of procedure were used (see Figures 8 

(a) and (b) for illustration): 

(a) The non-resonance portion of -the yield \jas estimted and sub- 
tracted fron the yield values to give only the resonance portion 
of the yield for analysis. TMs is estimted as the level vrhich 
the yield coi’ve seens to approach on either side. Since displace- 
nent of the non— resonance portion of the curve should not affect 
its hoi^t, the sane non-resonance datun level should be talcen for 
corresponding shifted and unshifted pealcs. In case the non-reso- 
nance base lines do not appear equal, some conpronise has to be 
reached. 

(b) If the pealcs appear Gaussian in sliar>e, \ic can ectii.-B.te the 
standard dev at ion quite acc'urately by takil^g a ha]! of the pealc 
\/idth at 60,6% of its height above the datun level. This is easily 
proved by reference to any table of the error function, or nornal 
distribution. 

(c) If the peak is not Gaus'ian in shppe, the standard deviation 
nust be conputed by nunerical analysis in accordance with basic 
definition. 

(d) The in;?luoncec of adjacent pcalcs nust be subtracted if they 

— 



Table XXXIV. (O\iro Pages) 



Computation of Straggling in Stopping 
Povsrer of Copper and Nickel Foils 



Line 


Date 


Foil 


/I I = 

Thickness 

(mg/cm^) 


Resonance 
Peak Ener, 
(kev) 


I» 

(amp) 


F» 


1 


8-29-52 


.05mil Ni-1 


. 1.150 


669 


.7848 


2.06 


2 


10-16-52 


,04 mil Ni 


0,977 


669 


.7720 


2.03 


3 


ff 


tl 


ff 


873.5 


.8620 


2.21 


4 


10-17-52 


ff 


tf 


630 


.7560 


2.00 


5 


tr 


ff 


ff 


986 


.9069 


2.29 


6 


10-23-52 


ft 


fl 


441.1 


.6681 


1.79 


7 


1- 9-53 


ff 


ft 


669 


.7765 


2.04 


8 


12- 2-52 


.05 nil Cu 


1.155 


340.4 


.6310 


1.70 


9 


rt 


n 


ff 


669 


.7860 


2.06 


10 


n 


n 


ff 

* 


073.5 


.8735 


2.14 


11 


ft 


ff 


1.18 


935.3 


.8970 


2.27 


12 


12-18-52 


ff 


»i 


669 


.7834 


2.06 


13 


n 


ff 


n 


873.5 


.8710 


2.13 


14 


ff 


ft 


1.155 


935.3 


.8960 


2.27 



35^ added to thickness because of carbonization. 



Line 


(kev) 


E 

m 

(kev) 


(ma) 


(kev) 




< 


1 


839.7 


757 


4.9 


10.10 


102,0 


39.1 


2 


809,8 


741 


4.13 


8.39 


70.4 


9.61 


3 


1006.2 


941 


4.5 


9.95 


99.0 


6.25 


4 


777.1 


704 


4.5 


9.00 


81.0 


4.41 


5 


1105.3 


1046 


4.75 


10.80 


118.4 


4.12 


6 


612.3 


527 


6.0 


10.74 


115.5 


28.1 


7 


009.1 


739 


4.35 


8.87 


78.7 


12.04 


8 


546,0 


446 


4.71 


8.00 


64.0 


23.8 


9 


834.6 


755 


4.71 


9.70 


94.1 


27.8 


10 


1018.6 


949 


5.10 


11.10 


123.2 


39.0 


11 


1071.0 


1006 


5.5 


12.49 


156.0 


42.2 


12 


829.5 


749 


4.4 


9.06 


82.1 


7.13 


13 


1017.0 


945. 


4.3 


9.16 


83.9 


6.30 


14 


1073.2 


1004 


4.6 


10.21 


104.3 


15.76 



— — 



Table XXXIV. (Concl.) 



Line 








il ^ 




iiV 

(kev2) 


-Of Ax 

(keT^-cffl^/mg) 


1 


62.9 


1.61 


2.59 


1.66 


0.93 


62.0 


53.9 


2 


60,8 


1.56 


2.43 


1.66 


0.77 


60.0 


61.4 


3 


92.75 


1.93 


3.72 


2.82 


0.90 


91.85 


94.0 


4 


76.6 


1.49 


2.22 


1.47 


0.75 


75.85 


77.6 


5 


114.3 


2.13 


4.54 


3.60 


0.94 


113.4 


116,2 


6 


87.4 


1.18 


1.39 


0.72 


0.67 


86.7 


88.8 


7 


66.7 


1.55 


2.40 


1.66 


0.74 


66.0 


67.6 


8 


40.2 


1.05 


1.10 


0.43 


0.67 


39.5 


34.2 


9 


66.3 


1.60 


2.56 


1.66 


0.90 


65.4 


56.5 


10 


84.2 


1.96 


3.84 


2.82 


1.02 


83.2 


72.0 


11 


113.8 


2.06 


4.24 


3.24 


1.00 


112.8 


95.5 


12 


75.0 


1.59 


2.53 


1.66 


0.87 


74.1 


62.8 


13 


77.6 


1.96 


3.84 


2.82 


1,02 


76.6 


64.9 


14 


88.5 


2.06 


4.24 


3.24 


1.00 


87.5 


75.6 







seen to be present to any extent. 

For the saJee of cor 5 )arison xjith theory, values of strangling are 
plotted in Fig^a*e 15. 

2, Discussion of Iii»rors . 

It is difficult to nalcc good estinates of probable errors of 
straggling results, since nany factors of indoteminable size are 
present. The best esticate of the size of the randon errors can be 
gained by noting the si^read in the results as tabulated and as 
plotted in Fig;irc 15. It vdll be noted that the values for copper 
arc ijuch none self-consistent tliat those for nickel. This my be 
due to the fact that the data for copper were in general talcen at a 
later period tlian the data for nickel, and the inprovenent in tech- 
nique and procedures obtained mth ’ ractice tended to mke later 
results norc precise. On the other hand, the anaint of data is not 
sufficient to insure that such self -consistency for the copper re- 
sults is norc than pitre cliance. 

Factors leading to large errors in the straggling results include 
such natters as are listed: 

(a) There is reason to think that the target LiF-1, uliich v;as used 
for sone of the earlier data is senevrhat uneven. Since the beam 
v/ithout the foil interposed is concentrated b.\t is norc diffuse v/ith 
the foil interposed, the portion of spread in corresponding pealcs due 
to target thiclsacss nay not cancel in the subtraction of the peak 
variances. 

(b) As indicated in Section III.B, very fine, -undetectable thick- 
ness -variations nay cause a stra :gling effdet in aildition to the 
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theoretical stopping pov/er variations. This error is always positive. 
If we arc nulling to accept theoretical predictions in stopping as 
roughly con'ect, the experinental results may indicate the lade of 
appreciable contribution due to this factor. (See Figure 15). 






V« Sirnnar?/' of Results . 



The results on stopping power are sumoarized in Table XXXV, 
giving the values for stopping pa;cr in kev-cm /ag, and in 
Table XXXVI vxhidi gives the results for atonic stopping pov;er ( O') 

p 

in ev-cn , Figure 11 gives in graphical fom the resul-fc for 

stopping pOfirer, \/hile Figures 12, 13, and 14 give in graphical 

1/3 

fom the atonic stopping poi/er results ejqpressed as ®/Z • 

The results in stopping poorer straggling, eroTressed as 

(A^) are given in Table XXXVII and Figure 15. 
dx 







I 



Table XXXV« (IV^o Pages) 



SUI.P'JfflY OP DATA ON STOPPING POVER FCH PROTOi^ 

(kev-on^/ing) 



(a) Nickel 






(b) Copner 






Proton 


Stopping 


Prob. 


Proton 


Stonping 


Prob. 


Energy 


Pov/er 


Error , 


Energy 


Pov;er 


Error , 


(key) 




c" 

- P 


(key) 




oK 

_ _ P 


527 


172.0 


1.5 


446 


169.0 


1.5 


704 


148.7 


2.0 


532 


158.0 


1.5 


718 


150,0 


3.0 


603 


150.1 


1.5 


739 


142.8 


2,0 


713 


139.9 


1.5 


741 


141.7 


2.0 


755 


137.0 


2.0 


755 


145.3 


1.5 


812 


131.6 


1.5 


757 


146.3 


1.5 


949 


120.0 


2.0 


915 


128.2 


3.0 


996 


114.7 


1.5 


935 


125.7 


2.0 


1006 


112.0 


2.0 


941 


133.7 


2,0 


1050 


113.1 


2.0 


949 


128.3 


2,0 








951 


131,5 


2.0 








977 


129.0 


3.0 








1000 


127.4 


2.0 








1007 


121.9 


2.0 


(d) Nitrogen (^Np) 




1046 


120.4 


2.5 








1047 


122.0 


2.5 


Proton 


Stopping 


■^rob. 


1057 


120.7 


2,0 


Energy 


Power 


Error , 



(key) 









406 


400 


4.0 


(c) /iTgon 






409 


401 


4.0 








452 


371 


2.5 


Proton 


Stopping 


Prob, 


501 


362 


4.0 


Energy 


Power 


Error, 


502 


351 


4.0 


(key) 






541 


328 


3.0 








648 


278 


4,0 


421 


224.5 


3.5 


648 


296 


4.5 


476 


224 


2.0 


716 


275 


4.5 


519 


215 


3.5 


717 


264 


4.5 


567 


207 


2.5 


748 


257 


3.5 


662 


176 


4.0 


914 


219 


4.5 


713 


187 


2.5 


914 


238 


4.5 


778 


176 


2.5 


942 


217 


3.5 


929 


154 


4.0 


973 


207 


4.5 


932 


155 


4.0 


974 


228 


4.5 


971 


154 


3.0 


1001 


209 


4.0 


989 


148,5 


4.0 









— 



Table XXXV. (Concl, ) 



(e) Heon 



Proton 

Energy 

(kev) 


Stopping 

Power 


Prob. 

Error 


422 


310 


3,5 


424 


290 


8.5 


517 


280 


3.5 


519 


269 


3,5 


668 


235.-5 


4,0 


730 


217 


4.0 


733 


217 


4.0 


929 


186 


4.0 


930 


191 


4.0 


988 


175 


4.0 


992 


188 


4.0 



(f) Krypton 



Proton 

Energy 

(kev) 


Stopping 

Pov;er 


Prob. 

Error 

el 

! 


420 


163 


3.5 


516 


149.5 


3.5 


527 


153.5 


3.5 


665 


130 


4.0 


675 


129 


3.5 


733 


124 


4.0 


743 


122 


3.5 


929 


107 


4.0 


941 


106.5 


4.0 


989 


102.5 


4.0 



L&L Xenon 



Proton 

Energy 

(ker) 


Stopping 

Pov/er 


Prob, 

Error 

f 


437 


120.5 


3.5 


441 


121 


3.0 


536 


110.5 


3.5 


538 


115 


3.5 


679 


91 


3.5 


679 


95 


3.5 


748 


86.9 


3.5 


750 


92.7 


3.5 


943 


77.8 


4.0 


969 


77 


3,5 






STMIARY OF DATA ON ATOUIC 
STOEPIIG POVfflR FOR ^OTONS (10“^ OT-cn^) 



(a) Nickel (b) Copper 



Proton 


Atonic 


Prob. 


Proton 


Atotoic 


Prob. 


Energy 


Stopping 


Error. 


Energy 


Stopping 


Error, 


(KeT) 


Power 




(Kev) 


Power 


%... 


527 


16.78 


1.5 


446 


17.86 


1.5 


704 


14.48 


2.0 


532 


16.70 


1.5 


718 


14.62 


3.0 


603 


15.87 


1.5 


739 


13.91 


2.0 


713 


14,78 


1.5 


741 


13.80 


2.0 


755 


14.48 


2.0 


755 


14.18 


1.5 


812 


13.90 


1.5 


757 


14.27 


1.5 


949 


12,68 


2.0 


915 


12.50 


3.0 


996 


12.11 


1.5 


935 


12.25 


2.0 


1006 


11.83 


2.0 


941 


13.02 


2.0 


1050 


11.95 


2.0 


949 


12.51 


2.0 








951 


12.51 


2.0 








977 


12.57 


3.0 








1000 


12.42 


2.0 








1007 

1046 


11.88 

Uc74 


2.0 

2.5 


(d) Nitrogen 


JlSd 




1047 


11.89 


2c5 








1057 


U.77 


2.0 


Proton 


Atomic 


Prob, 








Energy 

(Kot) 


Stopping 


Error, 








Power 


. 


(c) Argon 






406 


9.30 


4.0 








409 


9.32 


4.0 


Proton 


Atomic 


Prob. 


452 


8.62 


2,5 


Energy 


Stopping 


Error, 


501 


8.42 


4.0 


iKevl 


Power 




502 


8.16 


4.0 








541 


7,62 


3.0 


421 


14.87 


3.5 


648 


6.47 


4.0 


476 


14.81 


2.0 


648 


6.87 


4.5 


519 


14.28 


3.5 


716 


6.40 


4,5 


567 


13.71 


2.5 


717 


6.15 


4.5 


662 


11.67 


4.0 


748 


5.98 


3.5 


713 


12.39 


2.5 


914 


5,09 


4.5 


778 


11.66 


2.5 


914 


5,53 


4.5 


929 


10.19 


4.0 


942 


5.05 


3.5 


932 


10.24 


4.0 


973 


4.80 


4.5 


971 


10.20 


3.0 


974 


5.30 


4.5 


989 


9.82 


4.0 


1001 


4.85 


4.0 









(e) Neon 



Proton 

Energy 

(Key) 


Atomic 

Stopping 

Pov/er 


Prob, 

Error 

% 


422 


10,37 


3,5 


424 


9.70 


3.5 


517 


9.38 


3,5 


519 


8.99 


3.5 


668 


7.88 


4.0 


730 


7,27 


4,0 


733 


7.26 


■ 4.0 


929 


6.22 


4.0 


930 


6,38 


4,0 


988 


5,86 


4.0 


992 


6.30 


4.0 



(f) Krypton 



Proton 

Energy 

(Key) 


Atonic 

Stopping 

Power 


Prob, 

Error, 


420 


22.65 


3.5 


516 


20.75 


3.5 


527 


21.30 


3.5 


665 


18.07 


4.0 


675 


17.95 


3.5 


733 


17.21 


4.0 


743 


16.98 


3.5 


929 


14.81 


4.0 


941 


14.77 


4.0 


989 


14,22 


4.0 



(g) Xenon 



Proton 

Energy 

(Key) 


Atomic 

Stopping 

Power 


Prob, 

Error 

$ 


437 


26.30 


3.5 


441 


26.40 


3.0 


536 


24.10 


3,5 


538 


25.10 


3.5 


679 


19.85 


3.5 


679 


20.72 


3,5 


748 


18.95 


3.5 


750 


20.22 


3.5 


943 


16.99 


4.0 


969 


16,80 


3.5 






Table XXXVII 



SUWIARY OF DATA ON STOFPHKJ Pa.'ER STR/sGGLIKG OF 

PROTONS ( * (il2)) 

dz 

( kev^-cm^/mg) 



(a) Nickel 



Proton Energy 
(key) 



Straggling In Stooping ^ower 
(Variance per Unit Thickness) 



527 


88,8 


704 


77.6 


739 


67,6 


741 


61.4 


757 


53.9 


941 


94.0 


1046 


116,2 



Copper 



Proton Energy 
(key) 


Straggling in Stopping Power 
(Variance per Unit Thickness) 


446 


34.2 


749 


62.8 


755 


56.5 


945 


64.9 


949 


72.0 


1004 


75.6 


1006 


95,5 
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VI« Conclusion - Discussion of Results . 

A. Conparison of Prescait Results vdth Other Published Results * 

1« Stopping Pot.’er - CoTTparison v/ith Theory . 

In Ficnre 11, v;e lia.ve plotted mean curves obtained from our 
results on atanic stopping power. Superin^iosed on our resulting 
curves are curves obtained fron theoretical predictions. The set 
of curves on the right side of the figure are obtained use of 
equation (I-B-69), as modified inclusion of the correction Cj^ 
in the manner indicated try equation (I-B-83). The values of 
were obtained from llano (1934) > aiid are listed in Table II. ^he 
values of vrere obtained frcaa the graphs of Walske (1952), As 
this theoretical computation cannot be esq^ected to be very accurate 
in the energy region of interest to us , (primarily due to the 
failure to include corrections for deficiency in stopping pov/er of 
L and higher shell electrons), it is not suiprising that the theo- 
retical curves tend toxrard too high a value as the proton energy 
decreases. It is seen that a smooth curve of transition could easily 
be dra\n from ovir eijqjerinental curves to approach the theoretical 
curves at higher energies. 

The coup -ted o-irves of Hirschf elder and liagee (1948), which 
traverse the energy region of interest to \is, are plotted for the 
gases argon and xenon. Since their approach has included the stopping 
pov/cr corrections for all shells, one may cspect their curves to 
agree better idth experimental data than those based on the older 
formulations. In both cases the experimental curves lie belox; Jrhe 






theoretical c^irves ly about 5 to 10%, It would be unreasonable to 
espect better agreement (except ly accident), since the theoretical 
considerations on v/hich Hirschf elder and ^kgee based their analyses 
were admittedly approximate. 

Livingston and Bethe’s (1937) curve for proton atomic stopping 
po\/er in air is likewise plotted, as a rou;^ conparison ^d.th our 
results for nitrogen. Here, too, e^qjeriiaent seeras to be about 10^ 
below the predictions in this energy region. 

No eiTfort is made to include the predictions by Bohr’s (1948) 

formla for slow particles in heavy media (Equation I-B-24), since 

1 

it is not expected to give quantitatively correct ansirers. As a 

natter of fact, the fornula gives predictions for of 

2 2 
8.5 ev-cm at 500 kev proton energy and 6.1 ev-cm at IflX) kev. 

Ccanparison vdth Figure 14 shov;s that these values are consistently 
hif^er than experiment ly an appreciable amount. The theory can 
be checked hoirever for^correctness of dependence on Z and v. 

The plots of O" all seem to show that for heavy media the 

value of this variable is almost independent of Z. This conforms 
to the predicted variation of atomic stopping power v/ith atomic 
number - at least for elements of atonic number ^30. The in- 
verse dependence on v predicted by Bohr's theory can be checked 
by detcnnining from b^igure 11 wiicthcr the dependence on E is as 
the inverse square root of the latter variable. Since Figure 11 
is plotted on a log-log scale, the slope of the curve should deter- 
mine this e:p)onent. It is seen that all the curves have a slope of 
approximately minus one-half, in accordance vdth the theory. The 






heavier gases ore more nearly in accord, as is e^q^ected. The lij^ter 
gases shov; a steeper slope, and the netals a soiae^iiat less steep 
sloi>e« The anomalous results from the netals will be discussed 
in the concluding section, 

2. Stopping Poorer - Connarison \ath Other Experimental Information , 
The sources of other experimental infon.iation have already been 
indicated in Section I,C, Figures 12, 13 > and 14 give data on 

1/3 

0-/Z in such a v/ay that the comparison of our results with that 
reported from other laboratories is readily seen. 

In regard to the metals. Figure 12 includes results, attributed 
to Kahn, Our curve for copper is about 4^ belo;; that of Kahn, 
v;hich discrepancy is not alarming, Kahn has matched his curves to 
V7arshav;'s (1949) work at the 350-400 Icev region, but our curve for 
copper can be matched with V/arshav/’s just as easily, Ihe data of 
1-ladsen (1953) on copper is much different from other experimental 
data and is not included. The earlier \;’ork of Madsen and Venkates- 

I 

uarlu (1948a) on berylliian is essentially in agreement v/ith the 
curve given fcy Kahn, so that a separate curve for this contribution 
is unnecessary in Figme 12, althou{^ credit is certainly due to 
these earlier v;or]:ers. 

Our data for gases is compared in Fig;.Ti’e 13 v;itli the data from 
the California Institute of Technology up to 6OO kev, reported by 
Dunbar, et al, (1952) . The work of Ueyl (1952) on the lighter gases, 
especially on argon, up to the maxlnum cner.y he utilized of 4OO kev, 
agrees vri.th the data from the California Institute of Technolo^r, so 
that his rcaolts aro not ejSDlicitly included. It can be seen that 



our data a^ree v/ell vdth other results for the lighter gases. The 
data for krypton differ by a little more than the probable error, 
bat the results for xenon differ ly nuch more than the probable 
error. 

3« Conparison of Results of Straggling Measiu*cinents xdth Theory . 

There are no experinental measurenients of stopping poi/er 
straggling lath irfiich vre can corrpare oar results on nickel and 
copper. Since Madsen and Venkatesv;arlu>s (194Sb) neasurenents 
for beryllium checked so well with the theory of Bohr (194S), 
it is of special interest to check our results with theory. 

Equations (l-B^5) and (I-B-96) provide the relation, accord- 
ing to Bohr (1948), for theoretical prediction of straggling. 

They are readily evaluated in the cases of the metals nickel 
and copper. We should note that according to (I-B-96) the number 
of effective electrons, Z ’ , becomes about 27 for 500 kev protons 
and 39 for 1000 Icev protons. Since the latter figure is considerably 
more than the total number of electrons in the copper or niclccl atom, 
it is obvious that the estimate is too hij^. V/e have seen in sub- 
section 1, hovrever, that Bohr’s use of the Thomas-Femi model gives 
too high an ans\^er for stopping pov:er, iMch leads us to the suspi 
cion that the formula obtained from the model (l-B-72) gives too 
large a value. 

Evaluation of the theoretical predictions according to Livingston 
and Bethe (1937) is a little more difficult. (See equation I-B-97). 
The values of are not known precisely, nor is fr easy to evalu- 



ate. A calculation vao made indicating the contribution made ly the 
tem in which these constants appear, using 4/3 for each gnd 

the ionization energy in place of the average excitation ener- 
gy • This gives a value to this expression of about 8*5 for 
500 kev protons and about 11.5 for 1000 kev protons, valid both for 
copper and nickel. For Z’ we nust subtract froa t!"© atonic numbers 
of the element sbidid the oscillator strength of the K and L electrons, 
in accordance v;ith the criterion of equation (I-B-71). Using the 
values cenputed ty Honl (1933) for oscillator strength we get for 
Z’ values of 19.8 for nickel and 20.8 for copper. 

The theoretical curves and experimental data are plotted in 
Figure 15. It can be seen tliat our reafts are too scattered to per- 
nit us to nalce eny positive statements. The foUov/ing conclusions seem 
to be indicated, however: The theory of Livingston and Bethe is im- 

perfect in that it does not permit gradual changes in the value of 
Z» , the effective number of electrons; that of Bohr takes account 
of this change with proton energy, but does not include the extra 
term stemming from the kinetic energy effect of thecrbital electrons. 

On the other hand, the equation of Bohr over-estimates the number of 
effective electrons, and in this particular region this error nearly 
cancels the failizre to include the ld.net ic energy effects. The ex- 
perimental results roughly follow the trends v;ith energy as predicted 
ty Bohr (1948) though they seem to indicate that in absolute value, 
theory may give too high a prediction in this region. 






B. Concl;idinf: ConDonts , 

This dissertation ^ould not be concluded vdthout sone oorrnent 
on the general signifionce of the experinental results on stopping 
power in this energy region. There are a fev discrepancies between 
the results reported by the various ejsperimenters, including the 
writer, but nevertheless there is sufficient agreement to detect the 
pattern of the stopping pen/er effect at those moderately low proton 
energies. Figures 12, 13, and 14 - especially the latter - bring 
out a very interesting fact. Above 700 kev there is little differ- 
ence betv/een tlie stopping poi/er of gases and the solids studied ex- 
cept the predicted variation va.th atonic number, Belov: 700 loev 
however the results indiste a consistency anong the gases or anong 
the netals, b.vt a definite variation betv/een the tv:o groups of elo- 
nents. 

Theory loas not developed to the extent of permitting a definite 
explanation of this effect. In a qualitative v;ay, hov:ever, v;e nay 
speculate upon the possibility ihat those special, effects in con- 
densed nedia (oondufcting and non-conducting) discussed in Section 
I,B,6 nay provide the explanation desired. Whereas the polariza- 
tin effects in ccndaised nedia of hij^i atonic number are considered 
to have a negligible effect in viev-; of the snal!’ proportion of 
loosely bound electrons contributing to polarization, this presump- 
tion my no longer be valid at lov/ proton energies at which only 
the outer electrons contribate to the stopping power phenomenon. 
Since theoreticians have admitted that deficiency in stopping nay 



e 



exist for both conducting and non-conducting types of condensed nedia, 
v;e cannot state whether our resulting difference bet\;een the gases 
and nctals in stopping of3ou energy protons is due to the special 
conducting properties of the netals, or simply to their condensed 
state. It is proper and desirable, of covirse, to majee no final judg- 
ment on such a matter until experiments are forthcoming on the 
stopping pa;er of non-conducting elements in the solid state for pro- 
tons in this energy region. 

One conent nay be made on the stopping power of neon, v^ch 
according to Figure I 4 apjjcars to have an anomalously low value of 
stopping poTrer, This is not too surprising however upon some con- 
sideration of theory. The theoretical prediction of smooth vari- 
ation with atonic number is a i*esult of Bohr*s use of the Thonas- 
Fcrmi statistical model of the atom, which is v;ell-knov/n to be 
least valid for the light elements. Since neon contains only closed 
shells, its outer decbrons lie deeper in the Coulombian potential 
well than those of the atoms adjacent in the atonic table. Its 
electrons are tlierefore likely to be loss effective than those of 
the neighboring atoms, even aside from the general trend idth atomic 
number. 
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Basic Thooif/- and Deed,.'-?! Details of Volta'':e 
RorrixLator for Ib/giet Circait 




Basic Circuit of Regulator 
Fipure Apt>~1 



Load 

Rl 




Theoretical Triode Characteristic Curves 
Fifyre App-2 



Figure App-1 shox/s the basic circuit of the voltage regulator 
for the raagnet circuit. Figure App-2 gives the characteristic curves 
for an ideal triode. The tube used, 6AS7, has characteristic curves 
quite close to the ideal case — at least in region in viiich we expect 
to operate. The characteristics inay be theoretically expressed 




■fay the folloiriJig formula r 



_ Eg - (App-I-1) 

\^ere^‘is the tube arplification factor, Rp is the plate resis- 
tance of the tube, and the other ^Tabols are obvious from the 
above figures* 

From Figure App-1 ue con see tliat 

n I , (App— I— 2) 

and 

Eg = Ea - (I +Ip)(R3^ + Rg) . (App-I-3) 

Also we note that 

Eg - Ej_ - (I+Ip) Rg - Vg . (Api>.I-4) 



The above four independent equations can be reduced to a single 
equation l:^'’ the elinination of the three variables Eg , Eg, and * 
^ doing so ve my get an equation for I , which can be reduced to 
the folloiang form: 



I 



A + // V- -hr 




+ _S_Sfc 

Bl+Rj 



(App-I-5) 



with 

r s Rp - • (App-I-6) 

From these last tvro equations we note that I is independent 
of xjhen r z 0 , or t;hcn R^ r; V/j ■ Since tliis last 
factor is the reciprocal of tlio transconductance of the tube, the 
criterion that outiDut curreiit (and output voltage) be constant is tlaat 
conductance of the plate resistor equal the transcondLictance of the 
tube - whence the name of the regulator type: "transconductanco”. 
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Unfortunately no actiial tube characteristic curves correspond 
to the ideal curves, so that different operating points on the 
actual cliart have different values of transconductance* The proper 
value for is taken then as the reciprocal of triode transcon- 

ductance at some mean operating point located in the most linear 
portion of the actual set of curves* 

Once is selected then r may vary from zero only if the 

characteristic point at iriiich ire are operating is one ;d.th a dif- 
ferent transconductance value than the one used to select * If 
the actual characteristic curves are not too different from the 
ideal (this is one reason for selection of the 6AS7 tube - the 
others arc its high current capacity and its lo\/ transconductance, 
permitting a reasonably loa-; value for R^^), then r will not deviate 
markedly from zero. We see thus from equation (App-I-5) tliat I 
is affected very little by variations in r of the order of mag- 
nitude expected* 

We note also that if r is approximately zero, I is not only 
practically independent of but also of valuations in I ?2 • Thi 

resistance is not needed theoretically but has a very good use in 
actual practice. Assuming that r is small enough to be considered 
zero, we may use equations (App-I-2), (App-I-3), and (Apip-I-5) to 
deteivaino that: 

I - [k, . 1 

|_ ^ Rp + Rj. (1 +^) J 

(App-I-7) 

We see that ip does depend siffiificantly on R 2 . This enables 
us to vary I appreciably vdthout markedly changing I or . 

p - 2 



The value of tills will be seen in a monent 



The E^.'-steci nust be flexible enough so tliat any dosirea value of 
regulated current may bo obtained from about 0,4 ac^ieres to about 1*2 
ar^icres. (During the process of cycling the magnet, a naxiraum 
current of not less than 3 ainperes is desirable, bat precise regu- 
lation is not required in this process, and the design is made so that 
the regulator may be removed from the system by a single switching 
operation,) Tliis means that must be allowed to vary at vill 

irithout affecting the regulation effect. It happens that lines of 
constant transconductance can be drawn on the characteristic chart 
for the tube selected, v/liich lines proceed in general fi*om the 
lower left side of the chart to the upper right side. If the 
system is designed so that the operating point, even if it changes, 
stays on or near the line selected to natch , then the regular 
tion property remains for all operating currents. In case the sys- 
tem characteristics curmot be made to follow the proper transconduc- 
tance line (■v;iiich ue nay caill the operating line) exactly, changes 
in the value of R 2 made to bring one exactly to said opera- 

ting line. In practice it is quite easy to determine the proper value 
of to keep the system operating at the proper transconductance. 
With r sliglitly different from zero, any change in v/ill make 
a just perceptible change in current I . ^2 siiiqjly clianged until 

the position is reached for wliich slight changes make no appreciable 
clianges in magnet current. 

Gross changes in magnet crxrrent are effected tr/ cl^magcs in E^j^, 
wiiich is caused in turn by changing the setting on the field current 



- 



variac, as explained on Section III.C. In order to prevent the 
regulator fron opposing this desired current change, it is necessary 
to cliange sinultaneously some other factor in the equation (.^jp-I-S) 
on viiiich the current depends. Obviously V is the only one v;hich 
it is practicable to change. We provide this bias ty means of a 
bank of batteries across a potentiometer, the middle connection of 
which, being variable, can provide a variable value of V . This 
potentiometer is then "ganged" with the variac, so tliat they arc 
changed together, ^ proper selection of the potentiometer resistance 
and batteig’’ voltage, the value of grid bias will change along \riLth 
tlie plate voltage (E 2 ) so that the operating ppint moves approri- 
mately along the line of proper transconductance on the characteris- 
tic chart. 

The complete diagram of the magnet control panel, including 
all the control and regulator elements, is given in Figures App-3(a) 
and App-3(b), separated into the upper and lo't/er sections. This 
separation can easily be carried out in practice if necessary, 
permitted by the use of multiply-ganged connectors between the two 
sections. 
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